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INVESTIGATIONAL WORK CARRIED OUT IN THE USSR 
WITH REGARD TO ICE ENGINEERING PROBLEMS IN 

HYDRO-ELECTRIC CONSTRUCTION 

by 

B. V. Proskuryakov 
Professor, Doctor of Science 


Tho severe climatic conditions of the Soviet Union left their special 
mark on the design of hydraulic structures and on the construction and 
operation methods. 

The Soviet design-engineers, constructors, operating personnel of 
hydro-electric power plants and scientific workers of research institut- 
ions had to work hard in order to overcome numerous difficultiee which 
merciless nature created on Uie way of building and operation of 
hydraulic structures. This paper presents a brief description of the 
results achieved by the Soviet ice -engineers. 

Main cona tr uctio® methods in winter time . Excavation. 


Peculiar are the escavation methods used in the USSR, in view 
excavation of frosen soils being ImmeaQurably harder than that of that»y 
ones, the conetructors use as the principal means of warming soils ei&e? 
heat insolation in pre-winter period or special heating ©£ frosan ooils. 
Hero, in this country, a number methods and procedures ©f heating uiese 
developed and neceosary equipment was designed for this purpose. 

Wtotor cofficretii^ fee potpforfflfted in special housings of different sise and 
forms. Wide tise is m&de in the Soviet Union of tiie "^ermoe mediod" 
which c^sists in pre-heating off ingredients of concrete mixture before 
its pouring in blocks of structures with subsequent warming of concrete by 
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heat-ins dating cover. While concreting in cold weather use is also 
made oi electric heating o£ placed concrete, either of its exposed 
surfaces or of the whole mass of interior concrete. In building of 
industrial structures wide use is made of the method of steam cur^ 
of concrete. Within the last few years investigations have been 
carried out with regard to using in hydro-electric construction the 
method of cold concrete with chemicals which lower the temperature 
of water crystallization and provide by this hardening- of concrete at 
low temperatures. 

Methods of design of water temperatures in reservoirs and stream flow^. 


In the 20-iQ8 th© Russian researcher Rudovea used the equation o£ 
thermal equilibrium for the design of evaporation from the si^face of 
the Caspian Sea. It was the first attempt to establish interrelations 
between water temperature variations in a reservoir 

meteorological factors. This method of temperature design on the basis 
of thermal equilibrium was further developed by a well known physicist. 
Professor Academician V. V. Shuleykin. as a result of his investigations 
on the Kara Sea. A considerable stride forward on this way was made by 
engineer N.M. Bernadsky who in his works presents the eq^tion of 
thermal equilibrium in differential form and proposes methods of Its 
Integration. Moreover. N.M. Bernadsky’s water 
accounted for the design of water flows in reservoirs . for which 
he worked up a method of plotting the stream flow lines. The 
design of water temperature distribution over the whole depth stream 
flows and reservoirs were studied by the B.M. Makkaveev’s school by 
th® VNIIG (the All-Union Scientific Research Institute of Hydrotechmca) 
and by the Hydro-project. The combined efforts of these ins titutimo 
resulted in solving a number of practically important problems of ice- 
englsieeringa 

Ice processes in rivers, reservoirs and channel^ 

In the 20-ie8 large scale development of hydro-electric projects started 
in the USSR and in connection with this there were carried out 
investigations of winter regime in rivers and ^ 

eu«c©8S of hydro-electric power plants operation. We 
^ investigations carried out at the construction^ the 
Dnieper ^ Chirchlck end other power projects. The rosulta a£ 
lit^elugatlons are aummariaed la the works on i 

A.M. Istifeev, F.S. Bydin. V.E. Timonov a®i o^ers 
30-iea an Ice -Engineering Laboratory was organised In the 
Scientific Research Institute of Hydrotechnics in Leningr^. ^whlch the 
engineering staff made a considerable contribution to 

pr^lemt of Ice-englneering. Worthy of note are also the wurfe. of th® 
State Hydrological Institute (GGL) ^d the Hydroproject. 
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Aa the reoult of the investigational wos-k carried out by the above 
ino titutiono , w® have available at present the design prediction methods 
which permit uo to solve with an acceptable degree of accuracy th© 
problemo pertaining to fraasll ice regime, ice Jamming and clogging 
phenomena in rivers. Beeideo, the above mentioned inatitutioM w@?rksd 
o out the methods of combatting ice troubles during the construction and 
operation of power plants. 

Measure® were accepted for? 

a) controlling of ice processes upstream and downstream of 
hydro-plants during spring and autumn ice-drift. 

b) ice passing through hydraulic structures during the 
cosis tz'uction and operation periods© 

c) preventing metal structures of power house from freezing. 

The problems concerning the design of static and dynamic ice pressure 
against hydraulic structures are also thoroughly studied and recommend 
ations are given. Artificial frost is also widely used in hydro-electric 
construction. 

This large-scale Investigational work in the field o£ ice-engineering 
permitted also to study in detail the j^yslcal constant® of natural ic®. 

The afore-said Investigations were widely discussed in th© Soviet 
technical literature - in periodicals as well as in special publications. 
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ICE AFFECTING ENGINEERING STRUCTURES ON 
THE SIBERIAN RIVERS DURING ICE-RUN 
by 

K. N, Korzhavin 
Doctor of Technical Sciences 


1) Introduction 

At the present time in the U«S«S. R, widespread development of hydro- 
electric power projects and transportation facilities io in progress,. 
Particularly in Siberia, the effects of ice on these structures, during 
construction and operation, require special consideration. In the more 
severe climates, a considerable thickness of ice usually forms on the 
rivers, followed by a violent spring break-up, heavy ice runs, and 
frequent formation of ice jams and packs. The associated problems 
have been studied extensively, and the results of long-term investi- 
gation by the author are contained in this paper. 

2) Behaviour and Performance of 
Bridge Piers during Ice-Run 

Field observations of the action of ice on hydraulic structures and 
bridge piers in Siberia were made between 1934 and 1958, and the 
pro€8ss of ice breaking by piers during the 1945-1958 ice runs was 
filmed. Some of the results noted were as follows: 

W^n an ice cake meets a pier, the vertical 
the pier cuts it for a length of about 
0.5 meters causing local compression, 
aad ^rces are set up which either break the 
cake or bring it to a stop. 

An inclined ice-cutting pier face facilitates 
break-up of an ice sheet by cutting from 
below, and reduces the possibility of ice 
jamming. 
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The width of a atrip cut by the pier in an 
ice sheet does not exceed the maximum 
width of the pier® 

3) Mechanical Properties 

of Ice During lce~Run; 

In order to determine ice strength during an ice-run, the author has 
tested more than 1.500 ice samples, and has drawn the following 
conclusions; 

a) - At the beginning of the spring ice-run, compressive 
strength and resistance to bending of ice may reduce by 

as much as three times® ■“ 5 • 

b) - Local compression of ice due to pier penetration 
may increase interaction forces 2 to 2«5 times® 

c) - As the rate of deformation of ice increases, the 
compressive strength decreases® The following relation- 
ship was obtained experimentallyt 

R = R - - - (1) 1 

com 1 » 

which is true when > 0®! m/ sec® 

Design values for the ultimate strength of ice during an ice-run are 
given in Table I© 

Table I 


2 

Design values of ice ultimate strength, t/m._ 


Character 

of 

Vo 

Interaction m/ 

Rivers 

of the North and 
Siberia 


Direction 
of force 

® 0® 5 

sec 

Vo * 1®0 
m/ sec 

Vo ® 1®5 
m/sec 

toward the 
axis of 
crystal- 
lization 

1, Compression 

65 

50 

45 

Normal 

Zo Local 

bearing 

160 

125 

110 

86- 

3a Bending 

60 

45 

40 

Ql 
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Character 


of 

Interaction 


Vo 2 0,5 
m/ sec 


Vo 2 loO 
m/ Bee 


Vo = L5 
m/ sec 


Direction 
of force 
toward the 
axis of 
crystal- 
lise tioii_ 


4, Shearing 


30-40 


parallel 


5» Tension 


70-90 


tt 


Tests by the author on 73 ice samples showed that ice-run resistance 
is influenced greatly by the shape of a penetrating objects The form 
coefficient (m) for a pier triangular in plan is given by the following 
equations confirmed cxperiemcntallys 

m « 0.39 (1 + ) - - - (2) 

where oL is 1/2 of the pier nose angles 


In the author's opinion, existing standards for determining ice loading 
upon structures require refinementSs 

4) Proposed Design Methods? 


Assumptions are based on N, A. Rynin's supposition that the maximum 
possible ice pressure against a structure would be that necessary to 
break the ice or to bring it to a stop, 

a) - Piers with Vertical Cutting Face ? 

Ice pressure (P) developed at the moment of cutting a large 
ice sheet by a pier with vertical cutting face is the greatest 
possible and is given by the equation 

P « Kj, b^ h - - - (3) 

where h is ice sheet thickness 
bo is width of pier 

K- is a coefficient for local conditions 
” given by the expression 


2.5 K Rcom 
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where 

m 

k 

o 

om 

Vo 


is local bearing coefficient 

depends upon the form of the pier plan 

is the coefficient of contact 

is the compressive strength of ice 

is the velocity of the ice sheet 


b) “ Piers with Inclined Cutting Face ? 


When an ice sheet meets a pier with inclined cutting face, 
the ice may break due to shear bending* An analysis of 
observations made during break-up of an ice cover by 
piers having an inclined cutting fall produced the following 
relationships for determining ice pressure against piersi 


Breaking by Shear 

P o bo h K - - - (4) 

sino6 

where Pg is ultimate strength of ice in shear 


Breaking by Bending 


P 


0.66 Rfeh* 


^-|a_ + l2C0Sd ) 

12 sinoi - 


(5) 


where Rb is ultimate strength of ice in'bending 

is angle of inclination of cutting face 
to the horizontal. 

The equations for shear and bending were compared with the results 
of 13 field observations, and in 11 cases showed good agreement, 

5) The Method of Determining Ice 

Pressure in Field Conditions; 

When a large ice field is cut by a pier, the pressure (H) against the 
pier is given by the proposed equations 

H . c + £A (Vj - - - - (6) 

where c is given by the expression 
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In (6) and (7) above; 

£ is the coefficient of overall resistance 

to movement of the ice field 
A is the area of the ice field 

Vj is the velocity of the ice field 

Vo is the velocity of the flowing water 

(tj - t^) is the time of cutting 
m is the mass of the ice field. 


TABLE U 

Comparison of the results obtained by means of equations (6) and (7) 
with actual measurements are given in Table 11, 

Comparison of Actual and Design Ice-Pressure 

(in tons) 


ICE FIELDS 



Noa 

No. 2 

No, 3 

No. 4 

No, 5 

No. 6 

No. 7 

Determined 
by the 
Author's 
method 

72.5 

25,5 

54.5 

56.0 

66.5 

92.0 

88.5 

Actual 

Pressure 

45.8 

17.5 

46o 4 

45.3 

21.2 

62. 6 

55.0 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 



Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 

ZOA I 



Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 



Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


PROTECTION DES OUVRAGES HYDROTECHNIOUES 
CONTRE LA SOUS-PR ESSION ENGENDREE PAR 
L'ECOULEMENT DE L'EAU 


5ie 


PROF. R. CEBERTOWICZ PROF. AGREGE W.^ JAROCKI 
Laboratoire d'hydraulique et d'hydrologie k I'Ecole 
polytechnique de Gdansk 


sic sjc :{c 


SOMMAIRE 


L'article traite des evacuateurs de fond d'un barrage, ainsi 
que d'un clapet de vanne d'un barrage mobile oh, dans certaines 
conditions^ se produit une aspiration d'air, engendr6e par I'ecoule 
ment de I'eau, 

Les conditions dans lesquelles I'aspiration de I'air s'etait pro- 
duite ont 6t6 r^alisees sur des modules hydrauliques des ou- 
vrages precites. On y a observe la marche de ce ph^nombne, 
ainsi que des ph6nomenes associ^s qui 6taient accompagnes de 
bruits et d'6branlements de la construction. 

Afin de proteger 1' installation contre I'aspiration de Pair, divers 
dispositifs ont 6t6 am6nages sur modules hydrauliques; ces 
derniers 6taient destines h interrompre la couche d'eau afin de 
faciliter I'acces de Pair h I'espace par en-dessous de la dite 
couche. L'usage de ces dispositifs n'a pas donn6 lieu h. des 
solutions g6n6rales, vu qu'ils remplissaient leur fonction 
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uniquement dans des conditions fayorables, en d6pendance de 
l*epaisseur de la couche d*eau retotnbante, de la hauteur de la 
chute, de la position des dispositifs de fermeture, etc, 

Les etudes realisees sur modules ,r6duits ont donn4 lieu St la 
conclusion essentielle suivante! dans le cas ou l*aspiration se 
produit il convient tout simplemenit d^amener l*air d’une manifere 
artificielle k I'espace ou a lieu la sous-pression; toutefois, les 
recherches ulterieures concernant l*aeration naturelle des jets 
d*eau ne sont pas k negliger. 


INTRODUCTION 

Certains ouvrages hydrotechniques existants en Pologne accusent la 
naissance de chocs hydrauliques, une deformation des jets d*eau, des 
6branlements des constructions, etc, II a 6t6 constat6 que lesdits 
phenom^nes se produisent en raison de Uentrainement de l*air par les 
filets d*eau en ecoulement. 

D*article presente la description de ces ph6nom^nes en cours d*essais 
effectuees sur un module r^duit de pertuis evacuateurs et sur le module 
d*un clapet de vanne d’un barrage mobile, ainsi que Uetude des disposi“ 
tifs de defense centre les suites nuisibles desdits ph^nom^nes. 


ETUDE DES EVACUATEURS D»UN BARRAGE 

Le module d*un deversoir de barrage, ferme ^ Uaide de segments ainsi 
que des evacuateurs de fond, ferm6s ct I’aide de vannes, a 6t6 construit. 
La hauteur des Evacuateurs de fond en nature se chiffrait k 3, 5 m. k 
I’amenee et ^ 4, 75 m, h Uechappement (Fig, l). L*ecoulement maxi“ 
mum ^ travers le dEversoir se montait ^ 33 m, sec, par m^tre courant 
du deversoir, tandis que UEcoulement maximum ci travers un Evacuateur 
Etait de 180 m, sec. 

On a constatE qu*en dirigeant UEcoulement de l*eau h, travers Uouverture 
de 1* Evacuateur de fond seulement une partie de la section Elargie de 
Uechappement Etait remplie d’eau, la partie supErieure de la section 
etant occupEe par Fair, Le niveau de l*eau dans la section Elargie de 
1* evacuateur dEpendait de la hauteur de Uouverture de la vanne. 

l 

Dans le cas ou, au^dessus de 1* Evacuateur , le segment du dEversoir Etait 
ouvert, la couche d*eau s*Ecoulant cl travers le dEversoir recouvrait 
compl^tement I’ouverture de UEvacuateur , Dans ces conditions l*eau 
s*ecoulant ^ grande vitesse k travers le dEversoir et UEvacuateur en- 
trainait par aspiration Fair de UEvacuateur . 
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Si l*epaisseur de la couche d*eau au deversoir est plus grande que 0, 8 m, 
et la hauteur de Touverture de la vanne est inferieure 'k 20% de sa hauteur 
totale, alors, au fur et k mesur*- de 1’ entrainement de l*air par aspiration 
dans Tevacuateur, le niveau d*eau commence k monter, Au-dessus de la 
surface de l*eau dans I’evacuateur il se produit un rouleau hydraulique k 
axe horizontal (similaire au ressaut de Bidon) qui engendre des chocs 
hydrauliques et des ebranlements de la construction* La Fig* 2 montre 
le schema de la formation du rouleau hydraulique, et la Fig* 3, les 
diagrammes de la pression engendree alors dans l*evacuateur de fond* 

Si I’epaisseur de la couche d*eau au deversoir est inferieure k 0, 75 m* 
et Fouverture de F^vacuateur ne d6passe pas 2, 6 m*^ le rouleau hydrau 
lique ne peut se former au“desf'jus de l*eau s*ecoulant dans 1 evacuateur* 
Cependant k une difference suffisante de pressions d*air k I’interieur et 
k I’exterieur de 1’ evacuateur , I'air exterieur peut faire irruption k tra- 
vers la mince couche d’eau dans I’evacuateur en y engendrant des vibra 
tions de la construction et des bruits r etentissants* 


Divers dispositifs ont ete amenages sur le modblei aux fins deaeration 
de I’evacuateur. II s'est trouve que les interrupteur s des filets d’eau 
en forme de sauts en.ski ou de pyramides, instalies au-dessus de l*en- 
tree des evacuateur s, sont les plus avantageux (Figs* 4 et 5)* 

A I’aide d’ interrupteur s en forme de sauts en ski on a ete k meme ce- 
pendant d'obtenir I'aeration de I'evacuateur pour un ecoulement in- 
ferieur k 6 m* y seCa par in^tre courant du deversoir. Aux 6coule“ 
ments plus importants, sur les c6t4s du jet tombant du saut en ski, 
de minces couches d'eau se formaient en forme de pellicules qui in- 
terdissaient I'acc^s de I'air ci I'^vacuateur du c&te de l*echappement. 
Les inter rupteur s en forme de pyramides engendraient la rupture de la 
couche d'eau s'^coiilant cl travers le deversoir et I'aeration de I'^vacu- 
ateur aux 6coulements inf^rieurs ^ ZO m» /sec» par m^tre courant du 
deversoir (Fig. 6). Lors des 6coulements plus importants I'eau 
s'gcoulait au-dessus de I'argte sup6rieure de 1' interrupteur en couche 
d'eau continue, interdisant I'acc^s de I'air ^ 1' evacuateur. 


etude du clapet de la vanne au-dessus du barra ge mobile 

Un clapet ^ 20 m. de longueur et 1, 8 m. de largeur est install^ au- 
dessus de la vanne du barrage. La section du clapet transversal 
pivotant sur I'axe horizontal a 6t6 projet4e conform6ment ^ la courbe. 
Divers ecoulements d'eau k travers le clapet jusqu'k un volume de 
50 m, ^^sec. ont 4te pr6vus. 
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Lr*^tude effectuee sur le module d*un clapet de ce genre a montr6 que 
l*6coulement d*eau est accompagn^ d*un entrafneinent d*air par 
aspiration par en-dessous du clapet en y produisant une deformation 
des jets d*eau retombant et engendrant des ebranlements de la con- 
struction (Fig. 7). Aux ecoulements plus importants les jets d*eau 
viennent frapper la construction de la vanne et l*eau remplit Fespace 
^ air rarefie* 

Deux variantes d’un dispositif d'amenee d*air au-dessous du clapet 
ont ete projetees afin d’assurer I’aeration. Une des variantes 
pr^voyait 1* installation des interrupteur s ’k Uextremite du clapet, 

Uautre, Fusage d*un clapet k inclinadson variable de ses extremites 
(Fig. 8 et 9). 

L'etude sur module d’un clapet k interrupteur s a demontre leur 
efficacite lorsque Fepaisseur de la couche d*eau sur le clapet est 
inferieure k la hauteur des interrupteur s; dans le cas contraire, 
les interrupteur s donnent lieu cl une diminution de l*6paisseur de la 
couche d*eau; la rupture des filets d*eau ne se produit pas cependant 
et Fair n*a pas d*acc^s k Fespace par en-dessous du clapet# LVusage 
d’interrupteur s k hauteur sup4rieure> lors de la naissance d*6coule- 
ments plus importants, n*est gu^re possible pour des raisons pilre^ 
ment pratiques. 

L'usage de clapets k inclinaison variable de leurs extremit6s s*imposait, 
vu la supposition qu*en ce cas une disjonction des jets d*eau tombant du 
clapet aurait lieu, c*est-^-dire, que dans la partie du clapet k plus forte 
inclinaison Feau retomberait plus pr^s, tandis que dans la partie k 
inclinaison inferieure elle retomberait plus loin de Fouvrage; entre ces 
deux jets il se produirait un espace donnant lieu k Fa^ration par en- 
dessous du clapet. 

L* etude d*un clapet de ce type a cependant demontre que la disjonction 
de la couche de Feau retombante ne se produit pas immediatement, 
mais k une certaine distance du clapet, relativement k Fepaisseur de 
la couche d*eau. A une hauteur inferieure de la chute, la disjonction 
des jets d'eau ne se produit pas et Fespace par en-dessous du clapet 
n*est pas sujette ci Faeration. 

. CONCLUSIONS 

Les etudes realisees sur modules hydrauliques donnent lieu aux con- 
clusions suivantes: 

a) lors d*un ecoulement d*eau simultane k tr avers les evacuateurs 
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et les deversoirs d'un barrage, I'eau en s*ecoulant produit une 
aspiration d'air engendrant une sous-pression dans l'4vacuateur; 

b) ladite sous-pression engendre dans l'6vacuateur un rouleau 
hydraulique; si la couche d'eau sur le d6versoir est suffisamment 
mince une rupture de la meme couche peut se produire donnant 

^ acc^s el I'air dans I'^vacuateur; I'air qui fait irruption y restera 
jusqu'au moment d'une nouvelle aspiration; ces ph4nom^nes sont 
accompagn^s de bruits et d'4branlements nuisibles S. la construct- 
ion; 

c) Le rouleau hydraulique diminue la section utile de l'4vacuateur, 
ce qui donne lieu ^ un 4coulement ^ grandes vitesses dans 
I'^vacuateur ces vittesses sont quelquefois sup6rieures i celles 
qui se produisaient dans l'4vacuateur avant la naissance du 
rouleau hydraulique; 

d) el un ^coulement d'eau simultan6 ^ travers les 6vacuateurs et les 
deversoirs, le coefficient de I'^coulement dans I'evacuateur est 
ordinairement plus eieve, relativement au coefficient de I'^coule- 
ment calcuie pout le cas d'un deversoir ferm^ (inactif); 

e) I'usage des interrupteur s, resp, d'autres dispositifs d' interruption 
des jets d'eau et de I'a^ration, est impropre, car dans ce cas ils 
ne remplissent que partiellement leurs functions, k voir dans 
certaines conditions, relativement el I'^paisseur de la couche 
d'eau retombante, de la hauteur de chute, de la disposition des 
dispositifs de fermeture, etc; 

f) il apparaft done qu'il serait utile, dans le cas oh I'aspiration se 
produit, d'amener I'air par ,un dispositif special k I'espace oh la 
sous-pression a lieu; 

g) dans le cas oh I'air n'est pas amen6 au dit espace il conviendrait 
alors d'^viter I'^coulement simultan6 de I'eau el travers les 
4vacuateurs et les deversoirs du barrage. 
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Fig. 1 - COUPE TRANSVERSALS 
DU BARRAGE. 



L'EAU PAR LE^DEVERSOIR ET 
LES PERTUIS EVACUATEURS 
DU BARRAGE. 

a) d4but de la formation du rouleau 
hydraulique , 

b) rouleau hydraulique partiel, 

c) rouleau hydraulique total. 
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Fig. 3 - PRESSION HYDRAULIQUE 
SURGISSANT DANS LES 
EVACUATEURS DE FONDS. 

a) debut de la formation du 
rouleau hydraulique , 

b) rouleau hydraulique partiel , 

c) rouleau hydraulique total. 


9-D-7 


J If t 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 



Fig. 4 - DETAIL DU DEVERSOIR EN SAUT 

EN SKI AU-DESSUS DE L'EVACUATEUR. 

a) coupe en travers 

b) vue d'en face 



Fig. 5 - INTERRUPTEUR INSTALLE AU-DESSUS 

DE L'ECHAPPEMENT DE L'EVACUATEUR. 

a) coupe transversale 

b) vue d'en face 
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Fig. 6 - FONCTIONNEMENT DES INTERRUPTEURS A 

DIVERS VOLUMES D'ECOULEMENT PAR METRE 
COURANT DU DEVERSOIR. 
a) 6 m^/sec. b) 14m^/sec. 

c) 17 m^/sec. d) 21 m^/sec. 
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Fig. 7 - ECOULEMENT DE L'EAU A TRAVERS 

LE MODELE DU CLAPET DE LA VANNE 

DU barrage mobile. 

a) Q 22,8 m^/sec. b) Q 31 , 0 m^/sec. 
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Fig. 8 - MODELE^ DU CLAPET DE LA 
VANNE A INTERRUPTEURS. 



Fig. 9 - MODELE DU CLAPET DE LA 

VANNE A DIVERSES^ INCLINAISONS 
DE SES EXTREMITES. 

a) detail d*un clapet a forte inclinaison 

b) detail d'un clapet a inclinaison 
inf^rieure. 
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AIR ENTRAINMENT AND VERTICAL DOWNWARD MOTION 


25X1 


OF AERATED FLOWS 


^ ^ 


A. G. CHANISHVILI 
Master of Science, the Tbilisi Scientific 
Research Institute of Engineering Structures 
and Hydro-Electric Development 


^ ^ 


SUMMARY 


The experiments conducted on several models under the conditions 
of vertical downward motion of liquid with air entrainment yielded 
the following results. 

Within a wide range of the downfall jet velocity variation (I to 
ISm/sec) most of the bubbles formed in the flow are 5 -7mm*dia, 
when the normal settling velocity Wo can be considered as a 
constant value. 

An approximate theoretical analysis of the problem concerning the 
steady air bubble size in water agrees well with the experimental 
data. 

The main reasons for formation of unsteady aerated flows in vertical 
pipes are insufficient free aeration of downward stream, existence 
of local obstacles and reverse air flows. 

Minimum discharge of entrained air corresponds with uniform motion 
of the downward aerated flow. 
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Quantity of air entrained by water falling down the pipe depends 
upon water discharge, height of fall, etc. The velocity of downward 
flow causing entrainment of bubbles drops with increasing degree 
of aeration’ which is a result of their dynamic interaction under the 
conditions of decreasing settling velocity. 


A similar phenomenon has been observed in liquids with heavy 
content of solids. 


An analysis of the existing methods for taking account of the solids 
concentration effect on settling velocity shows that determination 
of the latter by suspending a granular layer in an upward flow yields 
the results closest to the actual conditions. 


Modification of this method for the case of an aerated flow enabled 
to investigate the effect of the degree of aeration upon settling 
velocity of air bubbles. 


The experimental results within the variation range of the aeration 
factors A from 0 to 0.30 can be expressed by the following empiric 


formula: 


where Vt = W.*Wo is relative settling velocity. 


SOMMAIRE 


Les essais effectues sur plusieurs installations experimentales en 
conditions du mouvement de liquide vertical descendant et I'entrai- 
nement d'air correlatif ont donn6 les r^sultats suivants. 


Dans une large 6tendue de variation des vitesses de jet (de 1 jusqu'^ 

15 m/s.) , la plupart des bulles qui se forment dans I'^coulement 
variant en dimensions entre 5 et 7 mm,de diametre lorsque la vitesse 
de chute Wo peut 8tre consid4r4e comme une valeur constante. 

L'analyse th^orique approximative du problfeme de stability des 
dimensions des bulles d'air dans I'eau est en conformite des essais 
sur modules. 

L'entrainement d'air dans l'4coulement en conditions insuffisament 
libres, I'existence des resistances locales et les mouvements d'air 
inverses sont les causes principales de I'origine du mouvement non- 
permanent des ecoulements a6r4s. 
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Le d^bit minimum de Pair entrain^ correspond au mouvement 
uniforme d’^coulement a6r6 descendant. 

La quantity d*air, entram^e par l*6coulement descendant dans 
un tube, depend du d6bit d*eau, de la hauteur de chute, etc* 

La vitesse de l*6coulement descendant, suffisante pour I’en- 
trainement des bulles diminue avec ^augmentation du degr6 
deaeration ce qui est caus6 par la diminution de la vitesse de 
chute au cours de leur interaction dynamique. 

Le phenomfene identique est observe dans les liquides avec 
des additions importantes de solides. 

L’analyse des m6thodes existantes de determination de I’effet 
de la concentration sur la vitesse de chute des additions, demontre 
que les r6sultats les plus rapproch6s des conditions r^elles sont 
obtenus par la m6thode de determination de la vitesse de chute 
au moyen de suspension de la couche granuiee dans I’ecoulement 
ascendant. 

La modification de cette methode pour les cas de l*ecoulement 
aere a permis d*etudier le problfeme de Leffet de degre d*aeration 
sur la vitesse de chute des bulles d*air, 

Les resultats des essais dans les limites de variation des coefficients 
d’a6ration A entre 0 jusqu'^0. 30 peuvent 6tre traduits en formule 
empirique: 

6.17A^-3.72AM 

oh W»Vs/o - la vitesse de chute relative. 


INTRODUCTION 

Air entrainment, motion of aerated flows and the process of deaeration 
are rather complicated and little known physical phenomena. 

There is no general theory developed for two-phase or multi-phase systems 
up to the present. The methods for analytical calculation of the phenomena 
involving aeration are not reliable and the methods of scale -modelling 
are little developed. For this reason it is almost impossible to generalize 
or extrapolate the few experimental results that have been obtained by 
various researchers at different times as well as to estimate quantitat- 
ively the results of theoretical analysis of individual problems. 
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The problem of aeration is of great importance in hydraulic power 
engineering. In most cases the aeration process results both for open 
and closed type hydraulic structures in negative effects causing consid- 
erable deterioration of service qualities of the structures and in some 
cases even leading to serious failures. In pressure type hydraulic 
structures the effects of the aeration phenomena are usually more 
pronounced than in pressureless structures. 

The problems of the aerated flow motion in hydraulic engineering, 
particularly the' air entrainment phenomena in water flowing down vert- 
ical shafts and pipes, were given much consideration at the 5th Congress 
of the International Hydraulic Research Association in 1953. 

^ 1) , 2) , 3) 

Several papers present the results of extensive research of 

some air entrainment phenomena and transport of air by water flowing 
down vertical pipes with flooded outlets. Some conclusions and results 
of the author's experiments are given in another paper. 

'The above-mentioned experimental research enabled to investigate the 
effect and significance of individual factors in the air entrainment and 
transport phenomena under specific conditions as well as to establish 
certain qualitative and quantitative relationships between them. 

The model experiments performed at the Tbilisi Scientific Research 
Institute of Engineering Structures and Hydro-Electric Development 
are generally similar to the above-mentioned experiments and provide 
a qualitative confirmation of their results. At the same time, these 
model experiments enable to define more exactly certain stream and 
structural peculiarities of a downward aerated flow which can be used 
to advantage for more profound investigations of the phase interaction 
phenomenon resulting in varying settling velocity of air bubbles, depend- 
irig upon the degree of aeration, as well as of other highly significant 
problems that are outside the limits of the particular case of motion 
under consideration. 

1) Marquenet, G. Entramement d'air par un ecoulement en conduite 
verticale. - Proceed. Minnesota Intern. Hydraulic Convention, 

Sept. 1953, p. 489. 

2) Kalinske A. A. Hydraulics of Vertical Drain and Overflow Pipes. 
University of Iowa Studies in Engineering, Bull. No. 26, p. 26. 

3) Laushey L.M. and Mavis F. T, Air Entrained by Water Flowing 
down Vertical Shafts. - Proceed. Minnesota Intern. Hydraulic 
Convention, 1953, p, 483. 

4) Viparelli M. Transporto di aria da parte di correnti idriche in 
condotti chiusi. -'T'Energia Elettrica", 1954, v. 31, No. 11. 
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CERTAIN PHENOMENA OBSERVED IN THE CONDUCTED 
EXPERIMENTS 


Detailed description of the experiments and experimental installations 
omitted, only those results are presented below that are of interest from 
the viewpoint of defining some aeration phenomena or are directly con- 
nected with the problem of settling velocity of air bubbles. 

The experiments aimed at the investigation of air entrainment by down- 
ward jet and the motion of aerated flow have been carried out under 
two different types of specific conditions. 

1) When the velocities of water flowing down a vertical pipe were con- 
siderably less than the settling velocity of the bubbles entrained by fall- 
ing stream (the experiments designed to investigate the penetration depth 
of air bubbles inside the liquid) ; 

2) When the liquid phase velocities were commensurate with the settling 
velocity. of separate bubbles (the experiments conducted on vertical pipes 
with the lower outlet flooded) . 

The experiments designed to investigate the penetration depth of air 
bubbles inside water have been conducted under the conditions described 
below. 

A jet of water was run into a vertical cylinder, 30 cm diameter and 3 m^ 
high, from two conoid nozzles, 11.1 and 13.2 mm*diameter respectively. 
When striking against water surface, the jet entrained some air. The 
height of free fall of the water jet varied from zero to 200 cm, rate of 
water discharge amounted from 0.19 to 2.25 1/sec. 

The experiments enabled to establish that under the conditions of air 
entrainment through jamming, at the atmospheric (gauge) air pressure 
in the point of contact, the overwhelming majority of air bubbles formed 
thereby within a wide range of jet velocities (1 to 15 m/sec.) are 5 to 7 mm, 
diameter. The aerated flow structure involves uniform distribution of 
the bubbles by volume. 

The penetration depth of bubbles entrained by the water jet, 

with the liquid phase moving through the vertical cylinder and fall remain- 
ing constant, increases with the jet velocity and rate of water discharge. 
With the fall increasing up to a certain value and rate of jet discharge 
remaining the same, the depth of air penetration suddenly drops, though 
it begins to increase slowly with further increase of fall. 

Use of product a *01 as a parameter gives a satisfactory scatter of 
points along the straight line (I, Fig. 1) expressed by the following equat- 
ion: 


H^= 0.052 u.d 
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In fig. 1 are also plotted the straight lines (3 and Z) obtained from the 
formulae of Abramovitch for the maximum and of Konovalov for 
mean velocity of a free turbulent jet, assuming that the bubbles penet- 
rate up to the depth where the jet velocity equals the settling velocity 
of separate particles and introducing some other simplifications. 

It is seen from the above figure that the intensity of the vertical flow 
velocity attenuation in case of air entrainment (non-free turbulent flow) 
is considerably higher than for a free turbulent stream under the same 
conditions . 

The experimental investigations of air entrainment by water flowing 
down vertical pipes having different diameters and working with flooded 
outlets confirm that pulsating motion of air and water mixtures with 
periodical entrainment of large air volumes is caused by insufficiently 
free aeration of suction zones as well as existence of local resistances, 
reverse air flows and so on. 

Motion of the mixture becomes steady under the conditions of free 
aeration and uniform air suction by the flow which is effected by supply- 
ing air to the suction zones, A structure formed in this case is char- 
acterized by uniform distriDution of bubbles which, for the most part, 
irrespective of the flow fall velocity or other factors, are 5 to 7 mm, 
diameter, as mentioned above. From the viewpoint of settling velocity 
such a mixture can be considered uniform (figs. 2 and 3) . 

Mention should also be made of another property of a uniform downward 
aerated flow confirmed by our experiments. 

It lies in that, with other factors being equal (water discharge, fall etc.) , 
minimum discharge of entrained air corresponds to the uniform motion. 

The experiments have shown that at a certain rate of the liquid phase, 
discharge velocity of liquid flowing down a pipe is such that the flow can 
ensure through entrainment of air bubbles, that is transport of the whole 
entrained air . 

1) Abramovitch G.H. Free turbulent liquid and gas flows. - Gosenergoizdat, 
Moscow-.Leningrad , 1948. 

2) Mackaveyev V.M. and Konovalov I.M. Hydraulics. - Techizdat, 1940. 


8-D-6 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 



Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


Air entrainment was observed also with downward flow velocities which 
were smaller than settling velocity of separate bubbles. Besides it 
proved to be that the greater was the degree of flow aeration the smaller 
were these velocities. 

This points out to the fact that, similar to the heavy solids content, 
the entrainment flow ability increases with the degree of aeration due 
to decreasing settling velocity of the bubbles or increase of actual 
velocities . 

The above-mentioned stream and structural peculiarities of the downward 
uniform motion of aerated flows somewhat simplify the general pattern 
of phase interaction, enable to reduce the number of factors to be taken 
into consideration and, similar to other two-phase flows, offer a possib- 
ility to investigate the problems of the settling velocity variation with 
concentration as well as of sliding between the phases and other phenomena 
becoming most apparent under specific conditions. 

APPROXIMATE THEORETICAL ANALYSIS OF THE PROBLEM CONCERNING 
STABLE SIZE OF AIR BUBBLES 


The problem of size of the air bubbles forming in the water flow aeration 
is of great significance for determination of the settling velocity, satur- 
ation, transport ability and calculation of power balance of the aerated flows. 

The attempts at theoretical determination of stable air inclusions known 
to us have yielded results which are in considerable disagreement with 
the author’s experiments. 

Consideration of the bubble stability subjected both to static ( P 5 ) 
and dynamic pressure heads enables us to work out the following 

equation: . 

p= p, Pd I- A p ,1, 

where ^ 

curvature. ( 2 ) 

If equation (1) is not observed that would mean transition of the bubble 
boundary surface to unsteady state. 

1 ) Levitch V. C. Physical and chemical hydrodynamics. - Akademia 
Nauk USSR, Moscow, 1952. 

2) Adam N. K. Surface physics and chemistry. - Gostechizdat , 
Moscow-Leningrad, 1947. 






is additional pressure due to the surface 
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For approximate determination of the bubble size and its ultimate strain 
causing damage of the separation surface it can be assumed that the 
latter phenomenon takes place when the radius of curvature in some point 
of the surface passes the infinity. In this case an infinitely small variation 
of the applied force will cause finite strain of the separation surface, that 
is, application of force will result in unsteady state of the surface. 


The point in question will be the lower point located on the vertical axis 
of symrnetry. 


If we do not take into consideration the dynamic pressure head and assume 
that the internal pressure does not depend upon the degree of the bubble 
strain, the following equation can be drawn from the equation of stability: 



iz 

1 -f 


2 

IRc 


Z<5’ 


26 " 


( 2 ) 


Assuming = 1 g/cm and ^ - 0.075 g/cm#positive values of Rb 
and Rd ^ can be obtained from this formula, with initial radii of 

an unstrained bubble Ro equalling 0 toQ. 38 cm. 

Taking into account the hydrodynamic pressure corresponding to a value 
of settling velocity equalling W© » 30 cm/sec.has an insignificant bear- 
ing upon finite velocity of a stable bubble. The tests designed to break 
the bubbles blown out from vertical pipes submerged with their lower 
ends in water have shown that the adopted hypothesis concerning breakage 
of the separation surface, when the sign of the curvature varies, is very 
close to reality. 


Therefore, both the approximate theoretical analysis and the experiments 
confirm that with the air and water mixture in uniform motion and the 
maximum interphase sliding velocity amounting up to 30 cm/sec. the stable 
bubble velocities are quite finite. 


Within the bubble velocities observed the settling velocity of their isolated 
motion is constant, i.e,, of the order 25-30 cm/sec, . 2) , 3) , 4) » which, 
from the viewpoint of settling velocity, permits us to consider the air and 
water mixture under the conditions Of our experiments as a homogeneous one. 

1) Levitch V. G. Physical and chemical hydrodynamics. - Akademia Nauk 
USSR, Moscow, 1952. 

2) Hasegawa T. USSR Abstracts Journal, Mechanics Series, 1957, No! 4, 
Abstr. 4212. 

3) Kozlov V. K. and Mologin M.A. USSR Academy of Sc. , Proceed. 

Techn. Sc. Dpt. , 1951, No. 8. 

4) Ladyzhensky R. M. , "Journal of Applied Chemistry" , 1954, v.27, No.l. 
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BERNOULLPS EQUATION FOR TWO-PHASE FLOW MOTION 

One of the most convenient and ready means used to determine the 
relationship between separate factors under the conditions of tworphase 
liquid flow motion is the Bernoulli's equation. However, correct 
determination of this relationship depends on the degree of approximation 
with which this equation was worked out, upon consideration of specific 
characteristics of the two-phase mixture motion as well as on the fact 
of how precisely each member of the equation reflects physical nature 
of the phjsnomena under study. 

Most precise definition of all basic members of the Bernouilli's equation 
for the two-phase flow motion has been given by M. A. Dementiev (1) ^ho^ 
took into account along with other kinds of power the solid phase pressure 
power absorbed by the flow, when deriving the equation. 

The member of Bernoulli's equation that defines the power of the dissip- 
ative forces is expressed by various authors by means of three methods 
differing in their physical nature. 

1. If the power losses are related only to the liquid phase unit weight, 
then the power of the dissipative forces can be expressed as follows: 

* Q-w * 

2. If the power losses are related to the liquid phase unit weight with 
fictitious discharge equalling the mixture discharge, then: 

|\J;= ' ClwS •hv7S2 

3. If the losses are related to the unit weight of the mixture, then: 

f\j ~ Q.WS KwSj 

The additional member of Bernoulli's equation will assume different 
physical meaning and various forms depending on what particular equation 
of the three given above is used to express the power of the dissipative 
forces . 

It seems to us that from the viewpoint of the phenomena under consider- 
ation the third equation is the best one since it allows for a possibility 
of utilization of the pressure energy and composition of mixtures. The 
above expression for the losses has also been used by A.N. Klimentov * K 

1) Dementiev, M.A. Proceed, of the All-Union Inst, for Hydraulic 
Engineer ing Research , v. 52. 

2) Klimentov A. N. "Hydraulic Power Engineering", 1953, No. 9. 
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If the rotai'y motion power, non-uniform distribution of velocities and 
pressures in both phases are not taken into consideration for uniform 
motion of the two-phase flow, Bernoulli's equation assumes the following 
form: 



The pressure drop Hws and specific loss Hw* related to the unit 
weight of the mixture are expressed in this formula by the fall of the 
liquid phase. 

If the mixtures are lighter than the medium, for example, air bubbles 
in water ( iJ'w ^a)» Bernoulli's equation for the downward ( Z< > Z21 
^ = Z< “Z2 ) uniform two-phase flow motion is written in the follow- 

ing form: 


or 


, taking into consideration that ta ^ ^ 


•wa n wa 


Hwa 


Qc 




z 


( 4 ) 


(5) 


The above equations make it possible to analyze the phenomena in 
question and establish certain relationships between individual factors 
peculiar for uniform motion of two-phase flows. 

ON SETTLING VELOCITY UNDER THE CONDITIONS OF INTERACTION 
BETWEEN PARTICLES 


The majority of researchers do not take sufficient account of variation 
of settling velocity of particles under the conditions of interaction between 
the particles or give erroneous interpretation of the physical and naechan- 
ical factors causing this phenomenon. 

The existing methods enabling to allow for the effect of concentration of 
the mixtures, that are heavier than the medium, upon settling velocity 
as well as physical and mechanical diagrams of the interaction of forces 
between the mixture phases can be divided into the groups described below. 
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According to Richards and Karaushev » variation in settling velocity 
is caused by increasing density of the mixture, the latter defined as 
some fictitious averaged density of a mechanical mixture. In fact, the 
medium surrounding a particle does not change with concentration, and 
introduction of the conception of averaged density within the range commen- 
surate with the volumes of separate particles in a mechanical mixture 
is not rightful. 

Finkey, Monroe, Velikanov and Paltchevsky (5) connect variation of 
settling velocity with increase of the liquid flow rate in the gaps. This 
diagrana is not sufficient to allow for structural peculiarities of a two- 
phase mixture, the dynamic phase interaction effect developing due to 
super-imposing of a velocity field forming around each article as well as 
for a correct estimate of the suspension effect. 

Thus, for example, the velocity of fall of suspended particles with respect 
to a fixed coordinate system under the conditions of dynamic equilibrium 
of these particles equals zero, though the settling velocity and energy 
needed for suspension are not equal to zero. 

Lyashchenko (6) , Mints and Shubert (7) substitute the conditions of the liquid 
motion through the granular layer in dynamic equilibrium for the conditions 
of contracted movement of particles and interaction between them. 

(1) Lyashenko P. V. Gravity Methods of Separation. - Gostoptechizdat, 
Moscow-Leningrad, 1940. 

(2) Karaushev A. V. "Hydraulics of storage lakes and rivers. - "River 
Transport", Leningrad, 1955. 

(3) Lyashenko P. V. Gravity Methods of Separation. - Gostoptechizdat, 
Moscow-Leningrad, 1940. 

(4) Velikanov M. A. Bed flow dynamics , vol. 11. - Gostechizdat, Moscow, 
1955. 

(5) Paltchevsky A. P. "Hydraulic Power Engineering", 1953, No. 2. 

(6) Idem see p. 15, nl 

(7) Mints D.M. and Shubert S.A. Hydraulics of Granular Materials. - 
Moscow, 1955. 
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The way this problem is formulated, except for a certain simplification 
concerning hydromechanical aspects, ensures maximum approximation 
to actual conditions of the phase interaction and provides extensive ex- 
perimental possibilities to obtain reliable and physically well-grounded 
values of settling velocity within a wide range of varying concentrations. 

To compare settling velocity in conditions of particle interaction accord- 
ing to the formulae recommended by various authors, in fig. 4 are plotted 
the curves iP(s) for sand possessing different grain size and a volume 

weight of 2.65 t/ra^. 

In case of the contracted fall-out of particles conception of settling velocity 
assumes somewhat different nature due to lack of coincidence between 
the absolute and relative velocities of fall-out. 

Mints and Shubert determined settling velocity \a/ 
as the averaged velocity of an upward flow W= iT = 
equilibrium of a suspended granular layer. 

From the viewpoint of correct allowance for the suspension work in dynamic 
equilibrium of particles this expression does not give rise to doubt. How- 
ever, to extend the conception of settling velocity equally well to the uniform 
upward or downward motion of the mixture, on condition of observing by 
the above principle, the average velocity should imply the average mix- 
ture flow rate: 

Hf- === Q. S 

^ws CO 

Then the settling velocity is determined from an expression: 

W *= ^ws — ^ 
vv = m ( \3w — 

where rn is the volume factor of water .content coinciding in meaning 
with volume porosity which is used with success when considering the 
problems of flow motion in suspended granular layers or porous mediums 
to give geometrical characteristics of the particular structure .( 1) 

In this case the aeration factor is defined by an expression A — 1~IT\ , 

while the reduced velocity and actual velocity of a liquid (u) are con- 
nected by the relationship: 

\)'= m-u 

In a particular case of dynamic equilibrium ( U's= 0 ) these express- 

ions give us W«= tn - while for suspending without sliding 
( "Vl'w*** ) i.e. 0 . The above expressions can be used equally 

well to determine settling velocity under the intermediate conditions of 
uniform vertical motion of two-phase liquids. 

(1) Idem see p. 16, n. 2. 


in their experiments 

ensuring dynamic 
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EXPERIMENTAL INVESTIGATIONS OF SETT LING VELOCITY OF AIR 
BUEBLES UNDER CONDITIONS OF THEIR INTERACTION 


Bei'noulli's equation for the uniform vertical motion of an aerated flow 
(4) or (5) gives an opportunity for experimental determination of 
settling velocity of air bubbles possessing various degrees of aeration. 


If the discharge in individual phases is expressed by appropriate flow 
velocities: — rn • U> ' and 


CL a A CO \Ja 


then: 


hwa “ Hwa ~ 


m +{< - rn) 


da 

Vw 




( 6 ) 


The condition of dynamic equilibrium for air bubbles in a flow implies 
that 


'wa 


H 


wa 


(7) 


If we take for settling velocity the value then formula (6) 

expressing the relationship between the loss, concentration and flow 
velocities of individual phases will assume the following form: 

hwa= Hwa- )Z 

The conditions most convenient for experimental determination of the 
degree of aeration are created iti case of dynamic equilibrium of the 
bubbles ( ) in a downward flow. Below is given an expression 

obtained for this case, if we do not consider friction losses between the 
flow and piping walls as compared to the interphase sliding resistance, 
and making an assumption that has been confirmed by recent experiments ^ 
according to which the pressure drop is always equal to the weight of the 
admixture column, the latter having a cross-section equal to unity. 

hwa- Hwa - ( f - If') Z (9) 


(1) Idem see p. 16, n2 
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Using the above expression, the air content, that is, the aeration factor 
i — rn ^ can be determined if the values of the pressure 

drop HwCL for the particular region Z are known. 

This method can be used equally well with uniform motion of the air and 
water mixture. 


To determine settling velocity of the bubbles in dynamic equilibrium 
using the formulae given above, experimental determination of the liquid 
discharge and drop is necessary, and in case of uniform motion air 
discharge through the particular cross-section should also be known. 

Experimental determination of the air bubble settling velocity was 
carried out using an experimental installation (fig. 5) consisting of a 
3 m^high and 10 cm.dia. transparent vertical cylinder. Water was 
supplied to the upward pipe end while its lower end was lowered in an 
air trap from which the air evolved was passed through a measuring 
nozzle (air gauge) calibrated in terms of pressure and provided with 
a micromanometer. 

Every 50 cm.of the cylinder length were mounted piezometer s enabling 
to determine pressure drops on the regions of equal length. The only 
measurements used to determine settling velocity were those obtained 
for the lowest region having a depth ^ — 50 cm. and characterized 

by uniform flow conditions of the mixture motion to which Bernoulli’s 
equation can be applied without introducing appreciable errors. 

Water discharge, as measured for the model, varied from I to 2,85 1/sec, 
while air discharge varied from zero to 650 cm^/sec. The aeration 
factor varied with pressure drop within the range of 0.06 + 0,30. 


The following parameters have been determined using the measured 
value’s of Q.vv , 0^1 and Hv 

a) aeration factor, 

A= 1 - m- 




Hwq 


b) mixture discharge and flow velocity, 
Q-w'*' Oa and V. 


wa 


CO 


c) 

d) 


velocity of air bubbles, 
finally, settling velocity, 




Q-g 

A CO 


W — — 




The experimental results (fig. 6) show that with smaller concentrations 
settling velocity tends to a value of the order 27 cm/sec.which is in good 
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agreement with the data presented by other authors. With the degree 
of aeration increasing, the settling velocity gradually drops to 12 cm/ sec. 
when A ==^ 0.30. 

Thus, within the range of the volume aeration values covered by the 
above-described experiments settling velocity decreases more than 
twice, as compared to its value for a separate bubble having the same 
size . 

The reason for such a drop of settling velocity is, as it is well-known, 
the dynamic interaction between the bubbles caused by the velocity 
field interference and variation of the flow factor. 

In fig. 6 is also plotted an experimental curve expressing the relation- 
ship between relative settling velocity and the aeration factor. Settling 
velocity of separate bubbles is assumed to be Wq 27 cm/sec. The 
curve is in satisfactory agreement with the following equation: 

H - 6.17 a“- 3.717 A + ) (10) 


The settling velocity drop which, as found out by means of the above 
experiments, takes place with the degree of aeration increasing, 
supplies an explanation to a number of phenomena observed in two-phase 
flow motion, including increase in the transport ability up to a cert^iin 
limit of saturation. 
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H, - 

Hg - 

H - 
K - 
t - 

a . 

CO 

'O' 

u 


d . 
D - 
N . 
P . 

H - 

PdL . 
aP - 

Z = 2rZa. 
m 


LIST OF SYMBOLS 
jet fall, cm, 

depth of bubble penetration, cm. 
pressure drop, cm, 
specific pressure head loss, cm. 
volume weight, g/cm? 

volume discharge, cm^/sec. or litre/sec. 
sectional area, cm? 
reduced velocity, cm/sec, 
actual velocity, cm/sec. 

The indexes W , $ and Cl denote the 

liquid, solid and gas (air) phases respectively. 

The indexes WS and CL denote a corresponding 
mixture 

jet diameter, cm. 

pipe diameter, cm. 

dissipated power 

internal bubble pressure 

external static pressure 

external dynamic pressure 

pressure due to the surface curvature 

difference in the sectional centre of gravity marking, cm, 

volume porosity or water content factor 

volume concentration of mixtures or aeration factor 
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W 


X 


w 

Wo 

cr 


settling velocity under contracted conditions, 
cm/ sec. 

normal value of settling velocity (for a separate 
particle) , cm/sec. 

relative settling velocity 

surface tension factor 


Ri )R 


2 " 


main radii of bubble surface curvature in a given 
point, cm. 


Ho 


unstrained bubble radius, cm. 



curvature radii in the upper and lower points of 
an unstrained bubble. 
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AIR ENTRAINMENT AND VERTICAL DOWNWARD MOTION OF AERATED FLOWS 



Fig. 1 - EXPERIMENTAL RELATIONSHIP Mg = • d- ) ; STRAIGHT LINES 

AGREEING WITH THE FORMULAE FOR FREE TURBULENT JET; 2 - 
ACCORDING TO KONOVALOV; 3 - ACCORDING TO ABRAMOVITCH. 
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Fig. 2 - STRUCTURE OF A DOWNWARD 
AERATED FLOW IN UNIFORM 
MOTION. 



Fig. 3 - SHAPES OF AIR BUBBLES IN A 
DOWNWARD UNIFORM FLOW. 
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Fig. 4 - RELATIONSHIP BETWEEN RELATIVE SETTLING 

VELOCITY AND VOLUME MIXTURE CONCENTRATION: 

1 AND 2 - MINTS AND SHUBERT EXPERIMENTAL 
CURVES FOR LAMINAR AND TURBULENT FLOWS 
(THE REGION OF TRANSITIONAL CONDITIONS 
SHADED); 3 - LYASHCHENKO EXPERIMENTAL 
CURVE; 4 AND 5 - CURVES FOR LAMINAR AND 
TURBULENT FLOWS ACCORDING TO PALCHEVSKY'S 
FORMULAE; 6 - CURVES AGREEING WITH VELIKANOV'S 
FORMULA; 7 AND 8 - CURVES FOR LAMINAR AND 
TURBULENT FLOWS ACCORDING TO KARAUSHEV'S 
FORMULAE; 9 - CURVE ACCORDING TO ZAGOUSTIN'S 
FORMULA. 
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Fig. 6 


EXPERIMENTAL CURVES: 


W= V’(A) 
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THEORY OF UNDERFLOWS IN STORAGE RESERVOIRS 

❖ 


1. 1. LEVI 

Professor, Dr» Techn, Sciences 
The Leningrad Polytechnical Institute named 
after M, I. Kalinin 


❖ >!< ^ 

SUMMARY 

The paper deals with theoretical problems of underflows origin- 
ating in storage reservoirs as a result of considerable quantities ^ 
of fine alluvium sediment carried by rivers. In actual conditions 
the regime of the underflow is either non-uniform or unsteady, 
but in previous investigations this problem has not received the 
attention it deserves. 

An equation is derived in the paper for the non-uniform underflow 
motion and for plotting curves representing the interface of the 
^ underflow and the water filling the reservoir; the equation is in- 

tegrated for a case of a positive bottom slope and a horizontal 
) cha^nnel, 

u 

The determination of the value of the total coefficient of resis- 
tance Ac for the underflow is of great importance; experimental 
studies that were carried out with this aim in view have allowed 
- the establishment of a regular interdependence of Ac and 

Reynolds number in the range of the latter’s values from Z500 to ^ 
35000, at a considerable turbidity (up to 200 kg. per m?), this fact 
now allowed the obtainment of a more precise value of this co- 
efficient than those so far obtained. 
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Section 2 of the paper is devoted to an investigation of the 
conditions under which an underflow originates* 

On the basis of theoretical considerations, as well as the 
data obtained experimentally a formula is established for 
the velocity gradient at the upper boundary of the underflow 
and equation (22) is derived for determining the minimum of 
turbidity at which the phenomenon becomes unstable. 

Section 3 deals with the problem of unsteady motion of the 
underflow; a system of differential equations is worked out 
for this case and the term is evolved expressing the velocity 
of the wave front under the conditions of the interrupted wave 
which is characteristic of the run-off floods. 


SOMMAIRE 

Dans ce rapport sont examinees quelques questions de la 
theorie des courants de fond (underflows) qui se forment 
dans les reservoirs des riviferes saturees par les plus 
petites particules. 

En r^alite le regime des underflows est ni uniforme ni per- 
manent; n^anmoins dans les recherches pr^cedentes ces 
problfemes n’etaient pas dQment consideres. 

Dans le rapport est d^duite Inequation de mouvement non 
uniforme de l*underflow et sont indiquees les differents formes 
de 1‘interface entre l*underflow et I’eau qui remplit le reservoir; 
I’integration de Inequation est executee pour le cas de la pente 
positive et horizontale. 

Pour calculer le phenom^ne il faut fixer le coefficient total de 
perte de charge dans le cas de ^underflow; les essais entre- 
pris dans ce but ont donne une fonction regulifere du nombre 

de Reynolds Re, la valeur duquel etait dans les limites de 2500 
^ 35000; la concentration de vase etait assez grande (jusqu*^ 200 
Kg/m.^); cela a permis de preciser la valeur du coefficient obtenue 
auparavant. 


Le deuxifeme paragraphe du rapport est consacre aux conditions 
de la formation des underflows. A la base des considerations 
theoriques et des essais,l*auteur est arrive ^ une formule pour 
calculer le gradient des vitesses ^ ^interface de Pinterflow et ^ 
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Inequation pour definir la concentration minimum avec laquelle 
le phenom^ne perd sa stabilite. 

Lie mouvement non permanent de l*underflow est examine dans 
le troisi^me paragraphe; l*auteur a etabli pour ce cas le syst^me 
des equations differentielles et la formule pour determiner la 
vitesse du front d’onde qui se forme dans le cas de I’onde inter- 
rompue qui est caract^ristique pour les grandes crues d^ondee# 


INTRODUCTION 

The experience accumulated in the course of the last twenty years in oper- 
ating storage reservoirs on the rivers possessing a high degree of sedi- 
mentation has led to the discovery of underflows in reservoirs* 

A turbid flow possessing a density considerably greater than that of 
water descends to the lower water layers and moves there in the direction 
of the dam* The peculiarity of this phenomenon is that the underflow is 
capable of transporting the finest particles suspended in water (d < 

10-20 jX) across the total length of the reservoir regardless of how great 
this length is (up to 150 km). 

The described phenomenon, which has been named the density current, 
was made an object of a number of field and laboratory investigations 
conducted in France, USA, Yugoslavia, Chinese People’s Republic 
and USSR^)»^)» ^)i The obtained data were used by Raynaud, Blanchet 
and Villatte in working out an approximate solution to the problem of a 
uniform steady underflow with a homogenous density fluid, rnoving.in tke 
water medium, as the base of the concept of the density currents. Under 
the conditions of turbulence an intermediate layer( y ) is formed between 
this dense fluid and water, within the boundaries of which the velocity 
and density gradients are quite considerable (fig* l). 


H, Duquennois* New methods of sediment control in reservoirs. 
Water Power. May, 1956. 

2 ) 

6tude des courants d’eau boueuse dans les retenues (underflow). 

4 Congr^s des Grands Barrages New-Delhi; 1951, Question n^ 

14. Rapport n^ 48. 

3) I. I, Levi. ’’Laws of a high turbidity flow motion in reservoirs.” 
Reports of the High School. Construction, No. I, 1958. 

4) 

' C. Blanchet et H* Villatte. Courants de density. Etudes exp6ri- 
mentales en canal vitre. Conference technique r^gionale. Tokio, 
1954* Laboratoire Dauphinois, Neyrpic, Grenoble, France, 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


8-C-3 


1 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


Proceeding from these assumptions Blanchet and Villatte produced 
the following equation'^/ 


(!f.- y)iiio - 


ift . /H 
Cl C?(H-h) 


( 1 ) 


In actual conditions the regime of the underflow is unsteady during 
4^5 months; run-offs are of a rain-water variety and are character- 
ized by a few waves of high intensity# The process of sedimentation 
of reservoirs in such or similar hydrological conditions can be con- 
ceived as follows: in the upper parts of the reservoir a cone of allu- 
vium is formed consisting of the largest size sand particles transported 
by the river (mainly of particles d :> 0. 02 - 0# 05 mm. Zone 1) ; dS 
soon as the cone of alluvium reaches sufficient height there appears 
an underflow; the finest sediment (d < 0.01 inm)V particles are 

carried downstream unobstructed^ while intermediate size particles 
are deposited a considerable distance below the cone of alluvium 
(Zone 2) (fig, 2). 


Due to accumulated flood waters, large masses of turbid water with 
high content of the finest particles of silt and clay come together in 
the lower depths of the reservoir; with the end of the rainy season, 
when underflows are formed no longer, these masses of turbid water 
flow down towards the dam, their level becoming horizontal and the 
sediment beginning gradually to settle down. Just this kind of situ- 
ation can be observed in the reservoirs on Lake Mead on the Colorado 
River, in the Kwan-Ting reservoir in China and in a number of res- 
ervoirs in Algeria# 


The presence of a clearly defined change of the shape of the bottom 
slope, in passing from Zone 2 to Zone 3, results in a non-uniform 
current regime coinciding with those periods of the year in which 
the river discharges become sufficiently steady, - a drop appears 
in the horizontal section and a backwater in the section upstream. 

In case of a speedy increase of the discharge and turbidity of the 
stream, the run-off water is propagated in the form of an interrupted 
wave and, in a short period of time, reaches the dam; here the wave 
is reflected, this phenomenon being accompanied by a sharp rise of 
the turbid current at the dam. The narrower the river cross section 




According to the opinion of French scientists the maximum 
diameter of the particles transported by the underflow reaches 
12 - 20 . 
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and the greater the accumulated discharge of the river - the greater 
will be the rise of the turbid current at the dam. 


This general outline of the process of underflow formation in reservoirs 
should be supplemented by an investigation of the conditions of the non- 
uniform and unsteady regime of the density currents, since only such an 
investigation would furnish the data necessary to form a correct idea of 
the current mean velocities for the total length of the reservoir, of the 
speed of the run-off wave front propagation, of the height of the density 
flow rise at the dam at the moment of the wave reflection, and of the 
nature of the accumulation process. Solution to these problems is abso- 
lutely indispensable for the purpose of designing the process of the 
reservoir sedimentation and of determining forces acting upon the struc- 
tures. 

It is at the same time necessary to make a thorough study of the con- 
ditions under which the bottom current is originated; observations prove 
that this phenomenon occurs only when a definite amount of sediment 5fnin 
is present in the river. Sedimentation being a known factor, the dura- 
tion of the period can be determined, in which the underflow will be 
extant, and the finest sediment particles carried down to the tail race. 

The present investigation is devoted to studying the problems indicated. 
The experiments described below were accomplished under the author’s 
guidance in the Laboratory of River Hydraulics of the Leningrad Poli- 
technical Institute by N. P. Kulesh, Master of Tech. Sciences, L.N. 
Alexandrov, Assistant, and A. K. Khapayeva, Senior Engineer, 

1. Non-uniform Motion of the Underflow 


In working out the equation for the non-uniform underflow, we shall 
proceed from the same postulates that Blanchet and Villatte had put at 
the base of their analysis; tangential stresses at the outer boundaries 


r- At tr ^ X 

of the current would be L *= 0 ^ — — while Ao would repre- 

sent the bottom resistance coefficient and A 4 - the upper boundary 
resistance coefficient. 


The slope of the boundary separating the underflow from the rest of 
the fluid mass will be , . 

oh 


1=1,- 


where J - the bottom slope* 


ds 
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Summing up all the acting forces including the force of inertia, and 
referring it to the unit volume (fig. 1, section 1-1 - 2-2) the following 
equation is arrived at: 



where i is the gradient of the reservoir free surface. Member 
is to be excluded from equation (2), For this purpose the theorem of 
pulses is applied to the flow as a whole, including the inverse current 
(section l-A-B-2, Fig, i) and the result referred to the unit volume: 


h A* cJ 

here cL s=s y - y Ll is local velocity. 

irh 

K the inverse current velocity in the upper layers of the reservoir is 
negligible oC differs little from oi; assuming oC - ot. in the first 
approximation and substracting (2) from (3): 


2g ( H - h) 


(4) 


whereupon equation (2) is transformed to appear: 


(5) 


Here 






Substituting iJ by ‘-p^where is discharge per channel unit width 
(rectangular bed being assumed) we have: 2 


ds 


( 6 ) 


y ^ 

Equating the denominator to Zero the critical depth of the bottom flow 
is found: 





(7) 
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Critical Velocity: 






( 8 ) 


Depth o£ reservoir + IoS ; thus, in the given conditions the uni- 

forrh flow regime is impossible* 

dh 


Equating 


to zero the minimum (or maximum) value of function 


h is determined and the following formula obtained: 

Ao+A, ^ 2c>tl< 


(9) 


where h-j - minimum depth, 




H 


H-h, 




To characterize A-c graph 2 (Re) (fig. 3 ) is given from which 

it is clearly seen that the resistance coefficient decreases regularly as 
Reynolds^number grows; in the experiments density was as high as 200 
kg. per m. ; this fact permitted to establish value in a sufficiently wide 
range of Reynolds Number, 

Values of Re > 24000 can be considered as related to the quadratic 
zone. This assumption, however, is subject to additional confirmation. 
Since H ■ f(s), also h© ■ f(s); variation of h© is negligible, however. 

Line h may change its curvature; equating the second derivative 
to zero, the value of.dh^^t the point of inflection D is determined; 


dh 


dl5 


D 


and since 


dH 

d5 

dH 

ds 


->,(h*-h;) 


( 10 ) 


y 0 9 with h ^ h, 


c> 


dh 


ds 


< 0 , 


3) 


It may be inferred from this formula that should a curve of drop be 
formed the curvature of line h is subject to variations. 

Configuration of the interface will vary depending on the correlation of 
h^ and h^, ( and Ic ). 

These configurations for 1^ ^ Ic shown in fig, 4, Lines N-N and 

K“K define three zones of the interface curves: 
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Zone A, where backwater is observed, 

Zone B, where drop is seen, and 

Zone C, where there is backwater. 

In the horizontal section interface curves of only B and C types are 
possible* 

The task is now to integrate the equation for the case of a positive slope 
of the bottom of the reservoir. New variables e h and 0^ r s 

he he 

are introduced. Equation (l) will then appear 


-I, 


where A a 


der 

Ac 

d. lo 

as H ; b 


A 


1 - 


5 


( 11 ) 


"7' 


(b-'nUo 

2cCI,(apl 

a (I + ^ ) 

Ao 


Since H a f(s), quantity A r f >6") ; solving equation (ll) becomes very 
difficult* An approximated solution is, however, possible with A assumed 
being constant within the boundaries of the investigated area of the reservoir* 
It is thus necessary to divide the whole of the upstream water into several 
areas in whose boundaries variations of quantity A are but negligible. 


Separating the variables and integrating equation (7) we have 

dV, 


lo 0' ■»’ K = 4.(A -i) 


T -A 


( 12 ) 


( d\ 

Integral I — — is easily reduced to integral 

well known in the theory of non“uniform motion, by substituting 




/ 


_ i 

'7’- A ^ 


cp(z) 


( 13 ) 
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We obtain the final solution in the following form; 


T 5a ~ S< 
he 




A-1 




( 14 ) 


In the process of the silting of the reservoir the river bed in front of 
the dam spillway openings is levelled and becomes horizontal; the case 
of the bottom flow moving over the horizontal bed is, therefore, of a 
particular interest from the practical viewpoint. The equation of motion 
will in the case take the following form: 




In the studied problem quantity H can be regarded as being constant 
since the river bottom is horizontal while the gradientj. of the free 
surface is insignificant. 


Separating the variables from the following equation and integrating 
equation (l5) we have: 





(■^^-1) {a-ij 
b - 





where 






b 


In the first approximation the method of summation can be applied in 
designing 

f sub-integral function 


Fig. 5 illustrates an experimentally plotted curve of a drop that was 
observed in one of the experiments conducted in the laboratory of the 
Leningrad Polytechnical Institute; there is, in the same graph, another 
curve plotted by designing with 

A<,= 0.005 and A, - 0.006. 
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2, Factors contributing to the underflow formation 

A steady bottom current can be formed only if certain definite conditions 
are extant; these prerequisite conditions are: depth h must be consider- 
ably less than H; since quantity h depends upon the flow discharge and 
its turbidity» it can be asserted that the turbidity should exceed a certain 
minimum quantity if underflow is to be formed. 

Should depth h of the flow be comparatively small, velocity u of the in- 
verse current will be low as compared to velocity of the density flow, 
while the intermediate zone will become negligible (fig. 1). With de- 
creased turbidity and increased h the inverse current velocity will grow, 
gradually approaching if ; simultaneously, the intermediate zone thick- 
ness will increase. Finally, the inverse current velocity gradients 
will approximate each other and the turbulence will spread to the total 
depth of the reservoir; at this moment the underflow will become un- 
steady, On the basis of this condition we shall proceed now to deter- 
mine the minimum density at which the underflow is formed. 

The experiments conducted by Blanchet and Villatte, on one hand, and 
by Kulesh on the other prove that the velocity curve in the intermediate 
zone is approximating the straight line; this is in full accord with the 
deductions of the turbulent stream theory, so that it can be assumed 
by analogy that 


= Const = 

dy y 


(19) 


since the intermediate layer thickness y_, provided » u, is equal to 

H-h . du 
T? 


H-h 


; here the velocity at the upper 


3 I I I I 

boundary of the underflow (at height h). 

The distribution of velocities in the underflow is ruled by the Logarith- 
mic Lawl); the variation of velocity in the central body of the cross- 
section is therefore insignificant so that tin can be as sumed to b e pro- 
portional to the mean velocity y , which is equal to « / Y * lo 

du V A( 


On equal grounds quantity 


tional to \/r’ and inverse^ proportional to VF 


can be considered as directly propor* 


1) 

Idem, 


p. 3, N. 4 
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In order to verify these assumptions special experiments were carried 
out at the Laboratory of River Hydraulics of the Leningrad Polytech- 
nical Institute. The results obtained are shown in fig, 6; two points 
related to the experiments of Blanchet and Villatte and borrowed from 
their paper have been plotted against this graph. 

On the basis of these data the rated correlation for 
represented as 


dh 


can be 


dh 


2,5 




( 20 ) 


We shall now introduce (20) into (l9) and express 


in the term of O' : 


H-h 


CO , where C K. 



1 . 


2.5 {it - 0.03)^[T 


reducing the expression we have 


h 

H-h 



( 21 ) 


Using fornaula (9) to determine h we write down the final equation to 

find min • 




0,8 C f 




loJT 


mih 




A( 


= 0,03 


( 22 ) 


r 


min 


where f - 



In this equation q, H, 
is a function of q, H, 
for two values of 


L 


)rTun 


and Ac are given quantities; hence 
and A.C • Values of jfmin and Smin 
: 1, O and 0, 1 are given in Table No, 
quantity *^0’* is taken as being equal to unity, which constitute a calcu- 
lation reserve* 


1 ^" 


Ac 


1 : 
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TABLE No. 1 


f 

lo/Ac 

B 1 

tt 

0. 1 


Tmln. 

^min.^s/m.^ 

t min- 

q 

min. 

0.01 

0.0011 

1. 85 

0.001 

1.7 

0.02 

0.0014 

2.4 

0.0013 

2.2 

0.03 

0.0017 

2. 9 

0,0017 

2.9 

4. 04 

0.0022 

3.7 

0.0025 

4.25 

0.05 

0.0029 

5 

0.0038 

6. 5 

0.06 

0.0039 

6.6 

0.0056 

9.5 


We are thus convinced that with an increase of i, e. with an increase 
of discharge ^ or with a decrease of K, turbidity necessary to maintain 
the steady motion of the underflow must grow; the effect of lo/Ac 
is slight. 

To corroborate these theoretical assumptions corresponding experi- 
ments were conducted, of which the first series was for the purpose 
of qualitative verification of the phenomenon and provided an evidence 
confirming the variation of the regime of motion in the spatial conditions 
on a plot of 10 x 13 m, (S, J, Goriunov, Master of Techn, Science, The 

A. 11- Union Scientific Research Institute of Hydrotechnics, named after 

B. E, Vedeneev), turbidity being 2-4 kg. per Tn? and _f being equal to 
0.02 - 0,04. The experiments of the second series (N. P. Kulesh, 
Master of Techn. Sciences), were to establish the greatest depth of the 
underflow, at which it becomes unsteady; value of "^^luctuated within 
0.35 - 0.45, turbidity being about 2,5-5 kg. per m?; obtained by calcu- 
lation, are shown in Table No. 2, 


TABLE No. 2 


la. 

Ac 

f 

^exp 

q ^ calc 

"^calc c 

^exp. 

1 

0.04-0.045 

2.5-3 

3.7-5 1.5-1. 7 

2.5 

0,065-0.08 

4-5 

6.5-8. 5 1.625-1.7 
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From the data shown it can be seen that with Cal the calculation fur- 
nishes a value of Smin which is in excess by 50-70%, this constituting 
a calculation reserve; it would naturally be impossible to expect the 
calculated and experimental values of Sniin to tally compl'etely in the 
case under investigation. 


3, Unsteady Motion of the Underflow 

The system of equations for unsteady motion of fluid as applied to the 
conditions of the underflow in reservoirs should be corrected in view 
of the fact that the density of a flow is a quantity variable in time. 


With variable ^ introduced in the calculation our equations will 
take the following form; 





g ds 






(23) 


Characteristic equation a f(t, s), which determines density varia- 
tion in time and over the current field, should be included in the system; 
as a rule, the density in the upper boundary cross- section is predeter- 
mined. ^ 


Going over to the new independent variables h andjt (as is usual with the 
analysis of an unsteady regime) and transforming the equations we . 
arrive at the formula for the velocity of the wave fron propagation : 

hence. 


c 


ds 

dt 



(24) 


1 ) 


S. Angelin and K, Flikstad determine quantity ’’C” by an analogical 
formula, but without the first, summand, 

S, Angelin and K, Flikstad, An investigation of intake arrangements 
for cooling water supply stratified Seawater, Alokarleby, Sweden 
VII Congress lAHR (Lisbon) C-13. 
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The second equation will take the following formj 

( !r, - y)!. gh ^ U5) 


(equation (25) is given for a rectangular channel). 

In case of turbidity being constant the equation is reduced to complete 
differentials and is solved in finite differences in conjunction with (24.) : 


dM- ±d(Ji,h)\P^+ Ac dt 

where jTd O' and h are variables* 


( 26 ) 


As can be seen from formula (24) the velocity of propagation of the wave 
front depends, to a considerable measure, on the turbidity of the flow. 
Together with the increased flood discharge and turbidity the velocity 
of the flood wave propagation will also be increased. 


Since the turbidity of the flow increases several times during rain falls 
(5-10 times on the Huang Ho river for instance) while at the same time 
the depth of the underflow substantially increases, the flood wave, gently 
sloping in the beginning, will rapidly increase its slope and will break. 
This entitles us to regard the rain flood wave as an interrupted wave. 

The above was confirmed by calculations carried out for a large storage 
reservoir (150 km. in length), on the Huang Ho River} these calculations 
proved that the flood wave is transformed but little and reaches the dam 
inrio less than 12 hours, the duration of the flood being 5-8 days. 

In analysing the problem of the underflow motion in the form of an inter- 
rupted wave (fig, 7) we shall apply the Law of Quantities of Motion to 
fluid sectionI~X -I-I ; it can be assumed then that during the period 
of time^lt the section will shift to position AA - BB, 


As a result of this shifting the increment to the fluid mass will be equal 
to AA' - CC (shaded): 

where ■ 0'4 and hj a initial velocity and depth of fluid respectively, 
iTjl and h2 ® velocity of wave front propagation, 

C -IT, s a, qi - -CTj hj , q^ - O' 2h2» 
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The velocity acquired by this mass • C hence the incre- 

ment o£ the section quantity of motion during At will take the form ofj 


A li AU = ^Q(h2'K)(c-t)-2)At 


Substituting ^ and b2“ hj for C-U 2 we have: 


C - Uo = 


2 -j— . 

^2 


KT2 “ •!/ J a a 


h2 - hj 


Hence, the increment of the quantity of motion will appear as 

' “ 1 ) h, a2 At 

gE^ I 


The impulse of the acting forces, with the bottom of the reservoir being 
h orizontal, will be: 

V 2 2 

^7 - hi 


P < At 


- ( 'f, ■* ^ ) 


.t 


Equating it to the increment of the quantity of motion we arrive at the 
final relationship: 




2 • 2 . 


h,-h 


(27) 


and hence 


= Ut ^ \/gh,lf (*+f F 


(28) 


Witl^^g being small, as compared with h^ ( Z ^0»4), qu^intity 


2hi2 


can be ignored. 


Formula (28) differs from that of the interrupted wave velocity in open 
channels by factor only. 

Formula (28) has been derived with the assumption that the thickness 
of the intermediate layer betv/een the underflow and inverse current of 
clear water is relatively small. In cases when this assumption is im- 
possible the calculated formula for the interrupted wave velocity will 
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take the following form: 



in deriving formula (29) the distribution of velocities and density in the 
interrupted layer is taken as variable in accordance with the Linear 
Law* 

In the event of the fluid discharge accumulation, the flood wave, reach- 
ing the dam, will be reflected from the construction, the height of the 
return wave will depend on the quantity of the accumulated discharge* 

In order to determine the height of the reflected wave we shall use the 
method of approximation offered by M, D. Chertoussov; the velocity of 
the front wave propagation will be determined in this case from formu- 
lae (28) or ( 29 ), the minus sign to be placed before the term which is 
under the radical sign (photograph 8)* The accumulated discharge A Q 
being known, we can determine the height of the wave in front of the 
construction at the initial moment of its formation from formula 

aQ.=»Q-Q. 2=' + “^2 j Bz (30) 

where B is the channel width in front of ihe construction at the height 
and B]^^ — width across the bottom. Velocity t / 2 is dependent upon the 
conditions of the flow approach to the construction; since the flood is 
propagated in the form of an interrupted wave, we have to calculate 
2 from formula (28); this calculation should be performed before 
the reflected wave is designed. It is obvious that this case will give 
the maximum possible value for Z, 

It should be borne in mind that the channel in front of the construction 
is usually considerably narrower than the channel of the upper pool; 
the reflected wave, therefore, in its subsequent movement, widens the 
plan and its height decreases. The filling of the bottom layers of the 
reservoir is rapid due to repeated reflections; the final level of the 
underflow depends on the total volume of the accumulated discharge 
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W and the form of the reservoir between the initial and final level ^ of 
the underflow which may be considered as being horizontal. The level 
of the free surface of the reservoir, which can also be regarded as being 
horizontal, will rise at the same time. 

In order to compute and Ay the balance equation is to be used in 

its integrated form, worked out separately for filling the bottom and 
surface layers of the reservoir. For the bottom layers (turbid flow) 
we have: 



where « bottom area at level Z, variable with 

increase of Z 

Q •• inflow discharge 

Q s discharge to the tail race 


Similarly, for the surface layers: 

T.adt-j Sidy 

ic 

where S s area of the reservoir water surface at level y • 

Since Q is considerably greater than ^ 2 * level rise of the free 
surface is much smaller than the rise of the bottom .current level* With 
the channel in front of the construction sharply narrowing, the initial risf 
of the bottom current level may prove to be the greatest and will not be 
exceeded in the process of discharge accumulation. In view of this fact 
the problem of the reflected wave calculation at the moment of its incep” 
tion is of great practical interest. To ascertain the validity of the assum- 
ptions the Laboratory of River Hydraulics of the Leningrad Polytechnical 
Institute undertook a number of experimental studies in the course of which 
the relationship between the accumulated ( Q) and spilt (Qz) discharges was 
modified and so was the turbidity of the flow. Distribution of turbidity by 
depth was non-uniform; due to this the turbidity in the reflected wave was 
much lower than it was in the direct, and its velocity was low (fig, 9); by 


Initial level is determined by calculation, in accordance with the 
formulae of non-uniform motion at the discharge proceeding the 
flood. 
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this fact high values of 


observed in the experiments, can be 


explained. With uniform distribution of turbidity by depth, values of 

z 

are found to be within the range of 0* 5 “ 1#5> depending on quantity 
•Ua 


value of z 


The results of the comparison of the calculated 
with the results obtained experimentally are given in (fig» 

10) where calculated and experimental values of ^ are shown on the 

abscissa and ordinate axes respectively; the velocity of the flow at the 
approach was near the critical velocity* Thus,, the data obtained by 
calculations based on formula (28) are satisfactory enough, accuracy of 
the result ebing within 10%» In most cases the values obtained by com- 
putation are exaggerated, the excess constituting a reserve of accuracy* 


CONCLUSIONS 

The obtained calculated formulas supply means for determining the main 
hydraulic characteristics of the underflow originating in reservoirs, 
provided there is a certain definite amount of sediment present in the 
river. 

These characteristics are: mean values of velocity (regime of the flow 
being steady), velocity of the flood wave front propagation, height of the 
rise wave in front of the dam with the flood waters partial accumulation, 
and the general rise of the underflow levels in the process of the flood 
waters accumulation* 

These characteristics being known one is able to design the process of 
the reservoir sedimentation proceeding from the scheme described 
earlier and to determine the quantity of fine suspended sediment that 
can be discharged into the tail race of the hydro-electric statioii* 

For this purpose it is first of all necessary to establish the minimum 
turbidity at which formation of underflow is possible; the duration of the 
period will thereby be determined, in which the fine suspended sediment . 
may be discharged into the tail race# In the remaining part of the year 
all the sediment will settle in the reservoir proper since velocities of 
the current will be very low. 

It is further necessary to analyse the problem of the volume of the 
accumulating flood waters, which is of importance in controlling floods 
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on one hand and in the matter of a more complete utilization of the water 
discharge at the hydro-electric stations on the other , and for the purposes 
of irrigation* 

The underflow velocities in flood periods are perfectly sufficient for silt 
and clay particles to be easily transported by the current; in the time 
intervals between the floods and at the time following their passage the 
current velocities will be considerably lower and especially so at the 
time when deposits of the finest sediment particles will be settling in front 
of the dam (Zone 3)* The deposits will cause a backwater in the upper 
areas, where current velocities may become so low that even the finest 
sediment particles will begin settling; it has been observed that settling . 
begins at velocities lower than 10-15 cm per sec# 

Making use of the equation for the non-uniform motion of the underflow 
in the horizontal area of the reservoir the zones and volumes can be 
designed of probable deposition of silt and clay particles, forming in 
the result of the backwater mentioned* 

Owing to the above calculation methods it becomes possible to draw up 
a general picture of the underflow formation, to establish the degree of 
intensity of the wave phenomena and the volume of sedimentation that will 
be carried by these currents into the river tail race# 

Designers will thus be enabled to determine, with sufficient degree of 
certainty, the date of the reservoir sedimentation and pressures borne 
by hydrotechnical constructions when waves of the rising turbidity flow 
appear at the site of the hydro-electric station# 


SYMBOLS 


h 

H 




Cq and Cj 
Aq and 


- mean velocity of underflow 

- depth of underflow 

- depth of reservoir 

- specific gravity of underflow 

- specific gravity of water 

- Chezy Coefficients for lower and upper boundaries 
of underflow 

- coefficients of resistance for the same boundaries 

t 

- tangential stresses 
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N 


CL and oL 
hp and "CTj 

" 1 = * 


Si 

V I and 1/ 2 


C 

Z 


- correctives of quantity of motion 

- critical depth and velocity 

- (a = ^ , q, specific discharge 

"c 

- channel width 
reservoir bottom area 

- reservoir water face area 

- velocities of current in the range 
wave front 

- velocity of interrupted wave front 

- depth after wave passage 

- height of wave: Z =h2 - hj 



of interrupted 
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Fig. 7 

Computed table for interrupted wave of filling. 








Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 



Fig. 9 

Dependence of the wave rising on discha^e 
accumulation (experiments). — = ) 



Fig. 10 

Comparison of the empiric data with 
designed wave of rising. 
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SUMMARY 


The majority of navigation locks built during the last years in the 
Volga River basin were preliminarily tested in hydraulic laborat- 
ories, as many problems of their hydraulics cannot be solved only 
by calculations and on the basis of theoretical analysis. 

Laboratory tests included mainly the following problems; 

; aj Determination of hydraulic characteristics of the lock culvert 

• ' system; 

\ b) Determination of hydraulic and operation conditions of vessel 

lying at the time of the locking and in approach canals; 

^ c) Improvement of the initial design supposition and lock struct- 

ures on the basis of laboratory testing and checking. 

' Field observations have been performed on several constructed and 

working locks, having as their aim the following; 

a) Determination of actual hydraulic and operation characteristics 
of built structures for their further use in future projects; 

b) Derivation of data for comparison of results of laboratory 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


9-B-l 



Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


tests and theoretical calculations with actual values for im- 
proving the procedure of laboratory tests and theoretical 
calculations. 

All laboratory and field observations carried out on locks were per- 
formed according to specially worked out procedure with the applic- 
ation of unique measuring apparatus. 

The paper contains the main results of tests of upper head and dis- 
tribution systems for filling and emptying lock chambers, as well 
as their comparison with the .results of field observations. Besides, 
the paper describes the methods of laboratory tests and field obser- 
vations, In particular it is shown that during dynamometric measure 
ment in mooring tackle, variation in mooring conditions introduces 
many uncertainties , while the application of the longitudinal trim of 
the vessel, as one of the principal factors of lying conditions gives 
reliable and certain data. 

The paper shows that as a result of research work certain relations 
have been determined between the conditions of vessel lying and the 
nature of wave phenomena, especially in approach canals, as well 
as several general conclusions on operation and hydraulic charact- 
eristics for various lock culvert systems. 


SOMMAIRE 

La majority des ecluses de navigation r^cemment construites dans 
le bass in de la Volga furent l*objet de recherches pr6alables dans 
des laboratoires hydrotechniques, certains aspects hydrauliques 
des problbmes les concernant n’ayant pu gtre r6solus uniquement 
par des calculs. 

Le programme des recherches en laboratoire portait essentielle- 
ment sur les questions suivantes: 

a) Determination des caract4ristiques hydrauliques des systbmes 
d*alimentation de l’6cluse; 

b) Determination des conditions hydrauliques et d*exploitation 
des bateaux en attente durant Declusage et dans les larges 
aval ou amont; 

c) Amelioration, ^ la base de la documentation obtenue en labor- 
atoire, des imperatifs d*etude et des systbmes constructifs 
des eduses. 

Des recherches sur place furent entreprises sur certaines eduses 
realisees et exploitees, ayant pour but de: 

a) Fixer des caracteristiques hydrauliques et d* exploitation 
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reelles des ouvrages dej^ r6alis6s afin de les utiliser dans 
les projets ulterieurs; 

b) Acquerir des donnees afin de conf renter les resultats des 

recherches obtenues en laboratoire et des calculs theoriques 
avec les resultats reels pour ameliorer les methodes des 
recherches en laboratoire et des calculs theoriques. 

Toutes les recherches sur les 6cluses en laboratoire et sur place 
furent effectuees d’aprbs des methodes specialement etudi6es avec 
emploi d*appareils de mesure originaux. 

Le rapport expose les principaux resultats des recherches sur les 
systbmes de tgte et de distribution, de remplissage et de vidange 
des sas et les confronte avec les resultats des recherches effectu4es 
sur place. Le rapport expose en outre les methodes des recherches 
faites en laboratoire et sur place. II est not6 en particulier que 
lors des relev^s dynamometriques dans les amarres, la diversity 
des conditions d’amarrage introduit de nombreux facteurs d* indeter- 
mination et la difference du bateau, en tant qu*un des indices 
fondamentaux de son comportement , permet d*obtenir des donnees 
plus concretes auxquelles on peut se fier davantage. 

Le rapport note que les recherches ont permis d*etablir des relations 
determin6es entre le comportement du bateau et le caractfere des 
ondes , surtout dans les larges. Ces recherches ont permis de 
tirer certaines conclusions generales sur les caract^ristiques de 
service et les caracteristiques hydrauliques de differents systfemes 
d’alimentation des ecluses. 


INTRODUCTION 

During the last 30 years in the basin of the Volga River and in artificial 
water-ways connecting it with the White, Baltic and Black Seas, over 
50 navigation locks have been built, including eleven locks on the Volga 
River proper. Ten more locks have already been designed and are now 
under construction. 

Several hydraulic problems which have arisen during the design of these 
locks cannot be solved only by theoretical calculation methods. There- 
fore, all culvert systems of these locks were tested (to different extent) 
on models in hydraulic laboratories. 

To find the degree of conformity of the model testing results with actual 
conditions, and also to obtain data on the actual hydraulic conditions of 
lock operation and conditions of vessel lying in chambers and in approach 
canals, field observations were performed on several locks with specially 
designed and built apparatus. 
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Because of the wide amount of investigation material this paper is limited 
by the description of procedure and principal results of hydraulic tests 
carried out for the typical feeding systems of locks, not touching upon 
the investigation of operation of hydro -mechanical equipment (lock gates, 
culvert gates, guard gates, etc.) 

All locks working with water heads in the chamber up to 13 m were de- 
signed with head feeding systems. 

In case of small variations of water level in the upstream of the ment- 
ioned locks, it was typical to fill the chambers through lifting lock gates 
located on the drop wall. The energy of the stream is dispersed in the 
chamber formed by shaped baffles (Fig. 1). 

For large variations in upstream water level and low drop walls the 
chambers are filled in most cases through short by-pass culverts 
located in the abutments leading into the baffle chambers and being 
closed by flat or cylindrical gates. 

At lock heads over 9-10 m, either multi- speed or variable acceler- 
ated opening of culvert outlets is widely used. 

The chambers are emptied, as a rule, through short by-pass culverts 
closed by flat gates and water discharges through outlet culvert sect- 
ions, partially or completely embedded into the lock bottom (see Fig. 1). 

At high water heads reaching 18 to 21 m in single -chamber locks and 
up to 27 m in the middle head of the double -chamber locks the large 
two -line Volga locks have distribution culvert feeding systems. 

Most of these locks are filled through longitudinal bottom culverts 
(Fig. 2). When these locks are located on derivation canals the water 
is also discharged downstream through a system of distributing culverts. 
When the locks are placed right next to the river bed there are two cases 
of water intake and discharge through the side structure from the chamber: 
in one case separately, i. e. the water intake is separated from water 
discharge; in the other, combined in a separate structure adjoining 
the spillway dam. 

The general problem of hydraulic laboratory tests of each lock is the 
improvement of the designed feeding system and the determination of 
such dimensions of separate parts of this system, which ensure filling 
and emptying of the chambers during a certain time conforming with 
the required capacity of lock. In this case satisfactory operation con- 
ditions of vessel lying during locking and satisfactory operation of 
individual parts of the lock should be fulfilled. 
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The conditions for vessel lying in the chambers during locking and at 
the approaches are estimated by the hydrodynamical forces acting on 
tlie vessel and ^lctually taken by steel ropes, used for mooring the 
vessel in the lock. 

The allowable values of the longitudinal and transverse components of 
hydrodynamical forces are standardized in the USSR in relation to the 
displacement of the vessels passing through the locks 

/S _ 1 to 1 \ 

Iw 1.000 2,500/' 

The maximum (average by section) water velocities in approach canals 
to the locks are limited by a value of 0,7 to 0. 8 m per sec. for ensur- 
ing safe movement and manoeuvrability of vessels during filling and 
emptying of the chambers. 

On the basis of these requirements for the lock feeding systems the 
following principal problems arise, which may be solved by hydraulic 
tests on lock models: 

1. Conditions of lock filling (direction and velocity of flow, shape of 
guide structures, limits of reinforcing and structure of the water 
intake itself). 

2. Duration of filling and emptying of chambers for the one or the 
other water distribution system in the chamber for various 
velocities and graphs of gate opening. 

3. Hydrodynamical forces acting on the vessel during filling and 
emptying of the chambers and their dependence on the speed and 
graphs of gate opening, discharge distribution system by length 
and width of the chamber, location and displacement of the vessel, 

4. Conditions for discharge of water downstream (distribution 

and value of flow velocity, shape of bank joints reinforcing limits, 
water discharge structure). 

5. Conditions of vessels lying at downstream approach of lock. 

6. Wave phenomena in approach canals and between lock canals. 

2. LABORATORY TESTING METHODS AND MEASURING APPARATUS 

A. Testing methods. Laboratory tests of navigation locks and approach 
canals were performed in the hydrotechnical laboratory of GIDROPROJECT 
and in the laboratories of other scientific - research institutes on models 
of various scale. 

For solving individual problems by testing locks, for example, determin- 
ation of the outline of the culvert and baffle chambers, determination of 
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the flow discharge coefficients of the system, etc, , tests were performed 
on fragmentary models reproducing usually to the scale of 1:40 the re- 
quired parts of the structure. This allowed a timely introduction of the 
necessary correctives into the design of the lock, using them to build 
the complete lock model. 

The scale of complete models for testing locks were taken in the range 
of 1:20 to 1:40. Models built according to these scales were employed 
for finding the principal operation characteristics: duration of water 
levelling, conditions of vessel lying in the chamber and at the approaches, 
operation of lock gates and culvert gates. Smaller models built to the 
scale of 1:100 were used for preliminary testing of the operation of dis- 
charge outlets downstream and water intakes upstream. 

Together with the development of laboratory research work, testing 
methods and measuring apparatus were also constantly improved. 

At present the hydraulic laboratory of the Scientific - Research Sector 
of GIDROPROJECT employs the following methods for testing locks: 

1. Determination of the hydraulic characteristics of the designed feed- 
ing system. Thus, the graph for water level variations is plotted, 
as well as the graph for gate operation and flow discharge co- 
efficients changing with time. 

2. Evaluation of vessel lying conditions in the chamber during locking, 
being performed by two mutually controlled methods: determination 
of forces by the dynamometric device of Prof. V.M. Mackaveyev 

and determination of vessel trim by a trimograph. The latter method 
proposed by engineer A.D. Khalturin has been found very successful 
both for laboratory and especially for field observations. 

3. Determination of wave phenomena in approach canals during locking. 

4. Introduction of required changes in the lock feeding system, determ- 
ination of the best conditions for gate operation from the point of 
view of favourable conditions for lying of vessels in the chamber and 
at the moorage on the approach sections. 

Laboratories are equipped with the installations usually employed for 
such tests and devices containing a number of weir chutes with large 
spillway front. Upstream and downstream are constantly fed for pre- 
venting variations in water level during intake and discharge of flow in 
the process of locking. 

In the process of model tests it was found that the water outlet conditions 
reproduced on models with different length of modeling sections for down- 
stream approach canals show different effect on conditions of vessel 
lying not only in the canal but also in the chamber when it is being emptied, 
due to the fact that the wave phenomena in the canal depends on the length 
of the latter. 
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In view of this a series of purely methodical tests were performed and 
measures were chosen allowing to get closer agreement of data and 
to create in the approach canals shortened on the model the conditions 
of wave phenomena approaching the actual conditions. These tests 
allowed by means of inclination of the chute water level regulator to 
obtain on the model any coefficient of wave reflection, both positive 
and negative. With the aid of these chutes it was found possible to 
disperse the discharge wave, having obtained on the model conditions 
corresponding to a canal of infinite (great) length, i.e. the case when 
this wave reflected from the canal outlet into the river bed does not 
have enough time to return to the lock during locking. 

The results of tests are preliminarily plotted as curves for water level 
time as well as flow discharge -time relations, hydrodynamical forces 
and speed of gate lifting, vessel position in a chamber, its displace- 
ment, initial depth in a chamber. This procedure for performing tests 
excludes possibility of accidental mistakes, as individual faulty tests 
are immediately discovered by the deviation of test points from their 
general disposal order. 

The comparison of records received on the trimograph and on the 
Mackaveyev device have shown that the vessel trim value, while the 
chamber is being filled changes similar to the readings on the Mackaveyev 
device, but with the advantage that the trim does not depend on the moor- 
ing system; this permits to use the trimograph readings for determ- 
ination of the longitudinal forces for any system of vessel mooring. 

When choosing the water discharge method into approach canals tests 
were performed on models of two scales; preliminary tests on models 
made to a scale 1;100 and then on the main model of the lock to a scale 
of 1:40. Velocity in the approach canals was measured both for steady 
hydraulic conditions conforming with the maximum flow discharge and 
for unsteady conditions during emptying of the chamber. On the basis 
of these measurements the water discharge system was chosen and the 
required reinforcing limits of the canal bottom were determined. 

The comparison of test results obtained on models of various scales 
allowed to find the acceptable scale of the model for solving one or the 
other problem. Thus, it has been found that the solution of problems 
connected with wave phenomena, distribution of water velocity in canals 
is possible on models of a scale 1:100 and even 1:150, while the determ- 
ination of operation characteristics of the lock (gate operating conditions, 
filling duration, conditions for vessel lying) are best of all found on models 
of a scale 1:40. 

B. Measuring apparatus. For studying the above mentioned phenomena 
the GIDROPROJECT laboratory has new unique apparatus allowing to 
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carry out automatic and high accuracy recording of rapid water level 
variations, changes in forces on vessels; values of instantaneous and 
average velocities, hydrodynamic pressures, etc. 

Most of the apparatus newly applied in the hydraulic laboratory have 
been worked out and made in the design department and shops of the 
Scientific - Research Sector of GIDRQPROJECT under the guidance 
of V. P. Bombtchinsky. 

For recording the water level variations in the chambers float lymni- 
graphs have been used, selfwriting on paper with the aid of a constantly 
beating high-frequency striker. Similar apparatus register the lifting 
of gates. Small fluctuations of water level are registered by special 
optical or pneumatic wave -meters. 

The forces acting on vessels while lying were previously measured by 
the Mackaveyev devices with spring-type dynamometers. At present 
a new device is being employed for evaluation of the conditions of vessel 
lying with dynamometers having induction pick-up units, their readings 
being recorded on an oscillograph (Fig. 3). The new device has stiff 
springs of low strain, this almost entirely excluding longitudinal move- 
ment of the vessel during operation of the dynamometers, which could 
reach a significant value (10 - 20 mm) when measuring by the Mackaveyev 
device. 

The velocity of steady water movement is registered by water -vanes 
and the usual plan photographing of illuminated floats. For unsteady 
water movement the surface velocity is fixed by a moving picture camera, 
or rapid photography with multi-frame film cameras, while the velocity 
in the water is fixed by a special ball-type probe with optical regist- 
ration of the two plane velocity components on the moving picture camera 
film. 

For recording the trim of the vessel the laboratory employs a mechan- 
ical trimograph, which registers on smoked paper in an increased scale 
the half difference of the altitude position of the vessel fore and aft, 
being transmitted to the instrument by thin metal threads. 

Hydraulic hoists are employed for lifting the gates, allowing to change 
the lifting speed in a wide range. 

All apparatus are controlled automatically from a special control board. 

This control board has two attachments with stop watches and a time- 
marker, allowing to synchronize registration on all apparatus, obtain 
special signals at the important moments of the test and to register 
automatically the time of gate lifting and equalization of water levels. 
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3. PRINCIPAL RESULTS OF LABORATORY TESTING 

A, Head water feeding systems. Initially, by testing the first locks 
on general models these systems were chosen both according to the pre- 
determined conditions of lying of vessels during locking as the main 
dimensions of the designed water feeding system, duration and speed 
of opening the culvert gates, and also according to the type and dimensions 
of dissipating devices, which determine as well as the hydraulic charact- 
eristics the time of filling -emptying of the chambers. This required, 
with a large amount of variables, many tests, which reached over a 
thousand in case of multispeed opening of the gates. 

Later, on the basis of analysis of extensive test data and obtained re- 
lations between the maximum hydrodynamic force on vessels in locks 
and at approaches to them and the hydraulic values characterizing filling 
and emptying of chambers, A. V. Mikhailov had worked out methods of 
hydraulic calculations of lock head feeding systems according to the 
allowable conditions of vessel lying during locking,^^ These methods of 
calculation, confirmed by a series of methodical tests, allow theoretically 
to determine the principal dimensions and maximum speed of opening 
culvert gates for the designed water feeding system, to plot their hydraulic 
characteristics and, finally, to find the minimum time required for filling 
and emptying chambers depending on operation conditions. 

The formulas derived by V.M, Mackaveyev and I.M, Konovalov and 
confirmed during tests also allow to determine the pressure drop behind 
water-tight culvert gates and those located in open shafts. 

Due to this, hydraulic testing of the Volga-Don locks performed by the 
laboratory of GIDROPROJECT from 1949 to 1951 corroborated the prin- 
cipal dimensions of the water feeding system and speed of culvert gate 
opening, which were determined by theoretical calculations on the basis 
of the predetermined time of filling -emptying the chambers and allow- 
able longitudinal hydrodynamic forces on the vessel in the chambers and 
at approaches. These overall dimensions were used for designing the 
structures proper and main equipment and mechanisms for. lifting and 
lowering upper lock gates and for culvert gates of lower lock head were 
ordered on the basis of the designed speed for opening the water outlets. 

In the hydraulic laboratory fragmentary models reproducing parts of 
locks and general models were also used for finding only the shape and 


2) A. V. Mikhailov, Head Feeding Systems for Navigation Locks and 
their Calculation, 1951, 

3) V.M. Mackaveyev and I. M. Konovalov, Hydraulics, 1948. 
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size of baffling and stilling structures when filling the chambers, for 
which the hydrodynamical forces on vessels in the chamber during water 
discharges increasing with time do not exceed the designed values of 
the first peak, caused at the beginning of filling mainly by wave phenomena 
in the chamber. 

The problem was studied of replacing the curvilinear and multi-speed 
opening of gates by the installation of slot sills before the gates lifted 
at uniform speed (proposed by the personnel of the hydraulic laboratory 
A.D. Khalturin, V.S. Mikhailovsky and N. V. Khalturina). The piers 
of the slot sill during the first period of chamber filling obstruct the 
opening (disclosed under the gate), thus allowing to receive the effect 
of the multi- speed curve for gate lifting by a simple single -speed mechanism 
(see Fig. 4), Making attempts to lower the time for filling the chambers 
in comparison with the designed height of the slot sill piers (0. 75 m), the 
investigators began to lower the pier height on the general model of the 
lock without changing the speed of gate lifting. It turned out that lowering 
of the pier height to 0.25 m did not increase the maximum longitudinal 
hydrodynamical forces on the designed vessels during locking, while 
the time for filling the chamber decreased significantly (up to 17%). 

Further lowering of the pier height lead to rapid increase of forces 
(by 1. 7 times) at the same gate lifting speed without the slot sill. 

The principal reason causing the rise of direct longitudinal forces act- 
ing on the vessel during lock filling is the wave phenomena in the chamber. 
The analysis of laboratory test data has shown that the smallest required 
time for filling the chambers may be obtained when the gate reaches the 
top of the piers when the wave front, being formed at the beginning of 
water feeding into the chamber, moving along the chamber and vessel 
inside of it is reflected from the lower lock gate. 

The time of wave run through the chamber with the vessel in it determined 
theoretically agrees well with test data. This has allowed in the future 
to find by calculation the pier height conforming to the given gate lifting 
speed. 

Slot sills are made on all standard locks on the Lenin Volga-Don Canal. 

The design of the baffle chamber at the lock head, consisting of a screen 
and beam -type baffle wall, was chosen mainly according to the conditions 
of lowering the maximum of reverse forces to a value not exceeding the 
first maximum of direct forces. 

When emptying the chamber through designed shape and size by-pass 
culverts and designed speed of gate lifting, in the laboratory were chosen 
the grate obstructions increasing towards the lock center line above the 
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discharge culverts, and also the shape of the mooring structures at the 
approaches providing the best conditions of lying and manoeuvring of 
vessels at the approach to the lock. 

A solid wall is built across the water discharge culverts along the lock 
center line. 

This measure allows to obtain satisfactory water spreading in the canal 
even when emptying the chamber through one by-pass culvert. 

It should be noted that the performance of model investigations of the 
conditions of vessel lying in the chamber and downstream approaches 
when emptying the chamber has allowed to come to the conclusion that 
opening of the culvert gates should be always designed with uniform 
speed for the head feeding systems, being opened duringi^^ =(0. 5-0. 7)T - 
duration of chamber emptying. 

B. Distribution feeding systems. The design of distribution feeding 
systems for two-line large -size Volga locks with heads in the chamber up 
to 18 m, due to the presence of large and rapidly growing load capacity 
should conform with the following important operation requirements: 
a) filling and emptying of chambers should not exceed 8. minutes and b) 
the maximum longitudinal hydrodynamical forces on the largest vessels 
should not exceed 1/2,500 of their displacement and for medium-size 
vessels correspondingly 1/1,500 of their displacement; the highest trans- 
verse forces should not exceed 50% of the longitudinal forces. 

Conformance with these requirements in the conditions when the maximum 
discharges of water entering the chamber during filling and discharged 
during its emptying exceeded 500 cu.m. per sec. with gradual and slow 
opening of culvert gates required long duration and extensive laboratory 
tests. These tests were performed initially in relation to the conditions 
of vessel lying during filling and emptying of chambers in the hydraulic 
laboratory of the Leningrad Polytechnic Institute (under the guidance of 
M. D, Chertousov and B. D. Khachanovsky) , and then on the general models 
of locks and approach canals in the hydraulic laboratory of GIDROPROJECT 
under the guidance of A. D. Khalturin. 

The performed tests of several different schemes of distribution feeding 
systems have shown that at great length and width of the lock chambers; 

a) the simple types of these schemes with short outlets from the long- 
itudinal culverts distributed along their entire length may comparatively 
slightly improve the conditions of lying of vessels during locking and 
slightly lower the time required for filling the chambers in comparison 
with the best modern schemes of head water feeding: 

b) for filling the chambers of the Volga locks with the volume of water 
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in the chambers up to 170,000 cu.m. for a short period of time at rigid 
conditions of vessel lying in the locks each of the longitudinal culverts 
should feed water only for a certain part of the length of the chamber 
through the discharge openings of the short transverse culverts. 

As a result of further laboratory tests the system for flow discharge 
distributing along the length and width of the chamber ensuring the 
given conditions was chosen (see Fig. 2). 

Regulation of water inflow along the length of the chamber is made in 
these locks by the corresponding location of inlet and outlet openings 
along the length of each culvert and time shifting of the beginning of 
gate lifting in various culverts (gate lifting advance). 

The distribution of flow across the chamber is performed by a system 
of transverse outlets formed between each pair of concrete distribution 
I-. beams (see Fig. 2). 

Water discharge into the downstream approach canals was also designed 
of the distribution type. 

The measurement of velocities in the approach canals was performed 
during laboratory tests both for steady-hydraulic conditions conforming 
with maximum discharge and unsteady conditions when emptying the 
chambers. 

On the basis of these tests a scheme of discharge culverts was chosen, 
as well as the required limits of reinforcement of the canal bed determined, 

4. FIELD OBSERVATION METHODS AND MEASURING APPARATUS 

For obtaining data on the actual hydraulic conditions of lock operation 
and vessel lying in chambers and approach canals, and also for check- 
ing the degree of conformance of model test results with actual data, 
the GIDROPROJECT laboratory had carried out field observations of 
some of those navigation locks which were designed on the basis of lab- 
oratory tests. 

For these field observations a method was worked out similar to that 
of laboratory tests; special apparatus have been designed and made. 

During field observations, work was performed on the following prin- 
cipal problems: 

1. Determination of hydraulic characteristics of lock operation, by 
recording of water level change in the chamber with time and cal- 
culation of flow discharges according to the increment of water 
level for each short period of time; 
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2. Investigation of conditions of vessel lying in the chamber during 
filling and emptying and also in the approach canal by the determ- 
ination of forces on the moorings, trim and longitudinal shifting 
of the vessel; 

3. Study of the mutual influence of two -line lock chamber operation 
with a common emptying system; 

4. Study of wave phenomena in the canals between the locks and at 
the downstream approach canals of the locks; 

5. Measurement of flow velocity in the canal at the end reinforced 
section. 

For measuring of water levels in chambers and canals the following 
apparatus were used: 

a) Lymnigraphs with four scales for measurement of water level 
and with four time scales; 

b) Suspended 15 -meter metal rods for readings of the position of 
water level in the chamber; the time of the reading was registered 
by radio signals sent every 15 seconds; 

c) Portable 1. 5 m fixed rods installed on the canal slopes for study- 
ing the wave phenomena in the canal; the time of the reading was re- 
gistered by synchronous stop watches, also every 15 seconds. 

The use of rods was very reliable and sufficiently precise; they allowed 
to organize observations in the most variable conditions, while their 
results completely coincided with the control records of the lymni- 
graphs. 

The evaluation of the conditions of vessel lying in locks was carried 
out by two methods: by determination of the vessel trim, and by 
direct measurement of the forces on the mooring ropes. 

The trim was registered by means of two systems of measuring units: 

a) Trimograph with heavy-weight pendulum and optical mirror system 
of registration of deviations in relation to its base plate; 

b) Trimograph with precise level guage and continuous photography 
of its position on a moving film. 

The comparison of records on these apparatus is given in Fig. 5. Both 
apparatus are placed on the deck of the vessel. 

The greatest difficulties were met with when measuring the forces on 
the vessels and transmitted to the mooring ropes. The actual system 
of vessel mooring cannot be performed in the same way as in the lab- 
oratory and the conditions of measuring these forces must be adapted 
to the mooring conditions usually employed during operation (mooring 
by two ropes at one board). The forces were measured by spring-type 
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dynamometers included in the mooring tackle with remote registration 
of their readings on an oscillograph; the trim of the vessel was re- 
gistered at the same time as the forces measured by the dynamometers. 

The trim, not depending on the vessel mooring system, may be direct- 
ly compared with the trim recorded on the model. Simultaneously, the 
trim in combination with the data on longitudinal vessel shifting allows 
to obtain material for control of the dynamometer readings in the moor- 
ing tackle, thus permitting to analyze the causes of separate jerks marked 
by the dynamometers. This comparison allows to judge whether the 
jerks were caused by filling or emptying of the chamber or whether it 
was caused by insufficiently tight mooring and sagging in the ropes or 
jamming of the mooring rings. In the first case maximum trim coin- 
cides with maximum force; in the second case the forces increase in 
the ropes only at small changes of the trim and significant shifting of 
the vessel. 

The device for recording the longitudinal shifting of the vessel during 
locking was either placed on the parapet of the chamber or on the board 
of the vessel. The device has a movable sighting carriage and rod with 
a centimeter scale. 

Common water wheels were used for registering the velocity in the 
canal and lock culverts. The direction of flow is registered by indic- 
ators set on the water wheel rods designed in the same way as wind 
vanes with electrical contacts. The lighting of lamp on the control 
panel allows to judge the presence of direct or return velocity. 

5. RESULTS OF FIELD OBSERVATIONS 

These observations have allowed to find the actual values of the flow 
discharge coefficient for various systems of filling and emptying of 
chambers and to compare them with the data of laboratory tests; 
besides, it enables to compare the data on hydrodynamic forces re- 
gistered by the dynamometers in the mooring ropes with the forces cal- 
culated according to the trim, and directly compare the trim records 
obtained in the field and in the laboratory. An example of these compar- 
isons is given in Fig. 6. Though the last curve of vessel displacement 
on the model and in the field somewhat differs, the general conformance 
of the trim changing can be well seen. When filling Volga-Don locks 
there is some difference probably connected with the different type of 
seal work when the gate passed through the slot sill piers on the model 
and in the field and with the presence of slight distortion of the gate in 
the field. 

During the testing of one of the locks, due to the presence of a perm- 
anent test vessel with the principal apparatus located on it, the con- 


9-B-14 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 



Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


ditions of operation of mooring tackle were checked when mooring the 
vessel by two ropes at one board and four ropes at both boards. In the 
latter case the lengths of all ropes were made exactly equal and two 
ropes were placed on winches for adjusting the rope tension. 

These tests have shown that for normal operating conditions two-sided 
mooring by four ropes does not achieve any advantage, first of all, 
because it is impossible to place the center line of the vessel parallel 
with the center line of the chamber, the forces are taken by any one 
fore or aft rope, but the friction of the vessel board against the chamber 
wall, which improves the lying conditions, is not used. For methodical 
mooring tests on four ropes there are some advantages, as this avoids 
the friction of the vessel against the side of the chamber. 

Besides the field observations on filling of locks and vessel lying con- 
ditions in the chambers and approach canals, the wave phenomena 
were studied in canals adjoining the locks. These observations were 
performed for canals of various length and different shapes in plane 
with one end at the river bed or reservoir or for those having locks 
at both ends. 

The obtained field data on waves phenomena in locks of canals of diff- 
erent length and shape allowed critically to evaluate the conditions of 
laboratory tests of locks, where due to limited space the approach canals 
were not fully reproduced on models. 

These field observations have shown that during laboratory testing it 
is necessary to consider the nature of wave phenomena which actually 
arise in the approach canals of the model, as these phenomena may in 
some cases be superimposed and alter the conditions of vessel lying 
not only in the approach canals but also in the lock chambers. 

When performing field observations of wave phenomena in navigation 
canals the coefficient of wave deflections was determined at the entrance 
of the navigation canal into the wide river bed at an angle of 30® to 40®, 
The obtained deflection coefficient was K=0. 35-0.37. The actual vel- 
ocities of the waves were obtained for canals having several small shallow 
one-sided widenings by leno'th. The speed of wave propagation in canals 
conforms with the Lagrange formula. 

The most suitable conditions for lock operation were determined in the 
field with the locks being located at the end of short canals between locks, 
eliminating the possibility of swinging them, measures being advised 
allowing during lock operation to take these swingings into account. 

On structures observations took place with oscillations varying in the 
range of 0,7 to 0,5 m. It should be noted that the evaluation of the con- 
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ditions of vessel lying in the field is more complex than in the laboratory. 
The dynamometric measurement of forces on mooring tackle takes into 
account phenomena which do not take place in the laboratory tests 
(speeding-up of the vessel, friction against the chamber walls, jamming 
of floating mooring rings insufficient rigidity of mooring rope fasten- 
ings on the vessel cleats, etc.). 

Therefore, field dynamometric measurements should be accomplished 
simultaneously by registration of the trim of the vessel and combined 
consideration of the obtained results, each of them correcting the other. 

The longitudinal incline of the vessel (the trim) may be recommended 
as the main qualitative index for field study of the conditions of vessel 
lying in lock chambers and approach canals. 

Of special importance for the trim is its independence of the conditions 
of vessel mooring and the possibility of direct comparison of model and 
field data. 

The field observations showed some of the complications in lock oper- 
ation obtained due to the presence of an emptying system common for 
both chambers. The combined system of water discharge was accepted 
owing to the considerations of obtaining the most uniform distribution 
of velocities in the canal when emptying any of the chambers. 

However, experience has shown that the combined water discharge 
system for lock lines introduces more complications into lock operation 
than unequal distribution of velocity in plane and it is not reasonable to 
achieve uniform distribution of velocity in this manner. 
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SYMBOLS 

H - Lock head 

ho - Initial water depth in chamber 

h - Head before upper lock gate 

S - Total longitudinal forces acting on vessel, 
tons metric 

W - Vessel displacement, tons metric 

tj - Duration of gate lifting 

T - Duration of filling or emptying 

I - Longitudinal incline of vessel 

t - Variable time during process of filling or 
emptying chamber 

K - Coefficient of wave reflection when reflected from 
wide river bed 

u - Gate lifting speed, m.per min. 
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SUMMARY 


We present here an approximated theoretical investigation of the 
fluctuations of water and of the conditions of ship stay in the chambers 
of locks and inclined ship elevators. Both these problems have rnuch 
in common and are solved here by the same method. 

The problem of fluctuations in lock chambers and in approach channels 
is treated in the first part of the report. Differential equations , 
with some restrictions, have been derived for water surface fluct- 
uations in a linearised form, with an approximate account for the 
movement of the ship itself. 'An example is worked out for the initial 
stage of the filling up of a head-fed lock chamber. As a result, the 
wave force acting upon the ship is determined. 

The problem of the inclined ship elevator is considered in the second 
part of the report. Approximated differential equations have been 
derived for the water surface fluctuations in the boat- carrying tank. 
The equations are integrated and the forces acting upon the ship are 
analysed. 


SQMMAIRE 

Cette communication est consacree a I'etude theorique des flottements 
de I'eau et aux conditions de stationnement des bateaux dans un sas de 
I'ecluse ou d'un ascenseur en pente. Ces deux problbmes 6troitement 
lies sont r6solus par I'auteur k I'aide d'une mSme methode. 
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Le probllbme des flottements dans un sas (chambre) et dans les 
canaux de l*ecluse est expose dans le premier chapTtre. On donne 
ici des equations differentielles linearis6es de flottements de sur- 
face de 1’e‘au en tenant compte de la presence du bateau. A la fin 
du chapiti’e est donne un exemple de la solution de ces equations 
pour le premier moment de remplissage du sas de l*ecluse (par 
amont) ^ Comme resultat de calculs on obtient la force qui agite 
le bateau. 

Le probleme hydraulique de Lascenseur de bateaux est etudie dans ^ 

le deuxifeme chapitre. En se servant de la m^me methode d*analyse <j 

on obtient des equations differentielles approximative s des fluctu- 
ations de niveau dans la chambre de I’ascenseur, On effectue 
^integration de ces equations et on analyse la valeur des forges 
agissant sur le bateau stationnant dans la chambre de Lascenseur. 


Many works have been dedicated throughout the world to the hydraulic 
theory of locks; however, up to the present time it is still insufficiently 
elaborated. The theory of inclined ship elevators was hardly even touched 
upon. At the present time the latter problem is of considerable interest 
to Soviet hydraulic engineers, in connection with the design of large ship 
elevators for the Bratsk and Krasnoyarsk hydroelectric power projects* 


I* THE PROBLEM OF OSCILLATIONS IN LOCK CHAMBERS 

1, Approximate Differential Equations for Water Surface Oscillations in 
Chambers and in Approach Channels 

It is known that the tension of a ship*s mooring ropes is, in the main, 
determined by the wave action of a fluid upon the ship, due to the oscilla- 
tions of the surface of the fluid when the chamber is either filled or 
emptied. 

When the ship is moored in a lock-chamber, the maximal tension of a 
mooring rope is registered usually in the initial period of chamber 
filling. 1/ To simplify calculations, it may be assumed that for the 
initial period of chamber filling (or emptying) , which is of the greatest 
practical interest, the depth in the chamber varies but little from the 
initial depth (h) . 

In analysing the unsteady movement of water and the conditions of ship stay 
in approach channels the surface fluctuations may also be considered to be 
small as compared with the initial depth, in this case, during the entire 
locking period. 


1/ This has been shown in the works of A, V, Mikhailov. 


8-B-2 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 



Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


These circumstances permit us to simplify substantially the in- 
vestigation of the most important phases of lock operation. 

For the calculation of the unsteady movementof water both in locks and 
in approach channels, with ships in them, A. V, Mikhailov applied the 
theory of shallow - water waves in its general form, i. e. , using non- 
linear equations, which were integrated by the method of characteristics 
We think that, in these most important cases of hydraulic calculation of 
locks, the shallow - water wave theory may be applied in its simplified 
form, by the use of linearised equations, which greatly simplify calcu- 
lations. That is all the more so as the non-linear theory is only an 
approximation, since the calculations of the ship vibrations oscillations 
are approximate. 

Let us consider the lock-chamber to be a prism with a horizontal floor. 
Let us also consider the approach channels of prismatic sections with 
horizontal floors. 


Let us so arrange fhe xyz system of coordinates that the xy plane should 
coincide with the fluid surface in a state of rest; let us place the x-axis 
parallel to the lock axis (or to the approach channel axis) and let us 
arrange the origin of coordinates •*0** at that lock head, from which the 
lock (Fig, l) is being filled (or emptied). 


Since in both the lock and in the approach lock channel water surface 
fluctuations may be only small, the solution may be based on a system 
of assumptions which are usually made in the theory of long wave of 
small amplitude. For this reason we will consider that ship oscilla- 
tions are also small. Considering only longitudinal waves, we may re- 
gard the component of the fluid particle velocity u^ to be considerably 
larger than the u^ and u^ components. 


Ignoring water viscosity and the lesser members in the equations of 
fluid dynamics we obtain: 


f dt ^ /i.i/ 




0 

a)z 


/l.z/ 


(quadratic velocity members are cast off, since the undulations are con 
sidered to be small). 
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From equation /1.2/ we obtain: 


'= Po + -z) po+ U) -2^)* 


Substituting this expression into the equation /l.l/ we obtain: 

^ _ Cl ^ 

?i t J 'd X 


/I. 3/ 


/I. 4/ 


In addition to the dynamic equation obtained let us derive an equation of 
continuity. 

Let us separate out a section by transverse cross-sections I and II 
(Fig. 2), with a distance dx between them. Let the free surface curve 
occupy the position a-a in the time instant t and the position a* -a* in the 
time instant t -h dt; at the same time the ship bottom gets displaced 
from position b-b into position b*-b*. 


During the time dt into the section under consideration through cross- 
section I there is an inflow with a volume u Cx)dt, At the same time 
through cross-section II there is an outflow with a volume 

dt. 


Due to the inflow of water through the longitudinal culvert, there will 
be delivered a volume qdxdt, where q " q (x, t) is the rate of flow enter-’ 
ing the chamber (or channel) per unit of time and per unit of length, 1/ 

Thus, the increment of the volume of water in the section during the time 
interval dt will be: 


IpiM. clxdt 
ax 


+ ^1 


clxdt. 


/1.5/ 


This volume increment in the section of the lock free from the hull is 
due to the surface undulation while in the section of the lock occupied by 
the boat’s hull also by the change of its position. It is equal to: 


dt dx 


71 . 6 / 


Equating expressions /l, 5/ and /l, 6/ we get the continuity equation: 

2)u} ^(Uxb3), ^ , 




/I. 7/ 


The longitudinal gallery outlets are usually arranged sufficiently 
frequently, and therefore the rate of flow entering along the length may 
be presented as a continuous function q of x. When water flows out 
through the gallery q is negative. 
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Let us note that in the free sections of the prismatic chamber 
while the sections occupied by ships OO = CO (x, , '^ ) • where 

is the ordinate of the ship's waterline equal to the rise of the waterline 
with relation to its initial position. Therefore, for free sections 


^co _ Oco _ L 

0 t 0^ 0 1 ^ 0^ 

and for sections occupied by the ship's hull 


/l. 8/ 


3 cd dtd 
0T ~ 


c)? . 0C0 _U 


at 

yt 


+ be 


0t 


/I. 9/ 


where b = b, ~ b : the width of the free surface of the fluid on a 
Ic c 

section occupied by the hull. 


It is easy to notice that the expression /l. 9/ is more general than /l. 8/ 
and contains it as a special case. Therefore, let us perform the 
further transformations with the aid of the expression /l. 9/» i. e. , for a 
section of a chamber occupied by a ship. 


The continuity equation will be represented thus: 



be 




dx 




/i.io/ 


In analysing the relative motion of a ship moored in a lock on ropes, the 
longitudinal shifting of the boat may be disregarded, since they are of no 
great importance but the mooring ropes have practically no effect either 
on the longitudinal oscillations of the boat about its transverse axis (the 
pitching motions) or the vertical oscillations. Therefore, the vertical 
and longitudinal oscillations of the boat have to be taken into account. 


In calculating vertical and longitudinal ship oscillations, which, as afore- 
said, are considered to be small we shall regard each of these oscillations 
as a sum of two kinds of oscillations: 


1) forced - with a period equal to 
the period of oscillation of 
fluid in the lock; 

2) natural - with a period equal 
to the period of the free 
oscillations of the ship. 

The period of the natural oscillations of boats may be evaluated by 
formulae known from the general theory of ship design (see, for inst. , 
/3/, pages 13-15 ) 
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Both theoretical and experimental analyses have shown that for 
ordinary practical conditions the oscillation period of a fluid in a chamber 
exceeds many times the period of the natural oscillation (both vertical 
and longitudinal) of the ship itself. For such a correlation of periods of 
oscillation of the fluid surface and of the ship natural oscillation, as is 
regarded in the theory of ship design (see, for inst. , /3/, pages 13-15 ), 
the amplitude of ship oscillations is close to the amplitude of water 
surface oscillations, i. e. , the ship follows after the movement of the 
water surface. 


That circumstance logically leads to an assumption, which greatly 
facilitates the solution of the problem: LET US ASSUME THAT THE 
ORDINATE OF THE WATER LINE AT ANY INSTANT OF TIME IS 
ON THE AVERAGE EQUAL TO THE ORDINATE OF THE FREE 
SURFACE OF THE FLUID : U 




/ 1 . 11 / 


By bringing /l. 11/ into the continuity equation 
approximate expression: 


W 




c)x 


/l. 10/ we will receive its 


/ 1 . 12 /. 


Let us note that, for free sections of the chamber, this equation, not being 
tied up with the assertion /l.ll/, may be regarded to be precise. 


Let us exclude u^ from the equations /1.4/ and /l.lZ/, 
we will represent the first equation in the form 



^ (Ux(^) 
C)t 




for which purpose 


/1. 13/ 


differentiate it along x and add to the second, 
t. As a result we will obtain: 


bk 







which is differentiated along 




/1.14/ 


If the area of the cross-section of the fluid in the chamber, which is in a 
state of rest, is a constant quantity ( CO = constant), then the equation 
/1.14/ gets transformed into an ordinary wave equation, 

I < 'S'V 

bn ’ n.is/ 


where the celerity of wave distribution along the chamber length 



1/ A more rigid analysis of this assumption shows that it practically is 
equivalent to the omission of the inertia of the ship. A more precise 
theory of fluid and ship oscillations without any assumptions /l. 11/ , is 
now being developed by the author. In this case the problem is resolved 
to integral-differential equations. 
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The equation /1. 15/ may be applied directly to free lock sections; 
here cO = and the celerity of waves 



If, in order to simplify the solution of the problem, we consider that 
the transverse cross-section of the submerged portion of the ship does 
not change along its length and that COo = const, and b^, “ const. , then 
the equation /1. 15/ may also be applied to lock sections occupied by 
the hull; here CO = CO|^ — CJc the celerity of waves 

The relation of the kind of /l. 18/, V. M, Makkaveev and A. V. Mikhailov 
obtained by another method. In the present work this relation, in the 
author's opinion, was fot;nded with greater precision. 

When the differential equations obtained here are applied to lock chambers 
(or channels) with concentrated systems of supply, these equations get 
simplified, since then q = 0. In calculating oscillations in locks with 
distributional sluicing systems the function q = q(x, t) must be first 
established, based on the calculation of the insteady flow in the distribu- 
tion system itself (in the longitudinal culvert). 

Let US formulate the boundary conditions for the integration of the equation 
/l, 14/ in analysing oscillations in lock chambers, whereby, for simplicity, 
we shall limit ourselves to the case when there is only one boat in the 
chamber (Fig, l). 


1 / 


In deriving equations /l, 14/ and/ 1,15/ and transforming /1.13/, we 

considered ,, ^ ri. 

Ux U) 

based on the accepted assumption that the oscillations are small. Let 
us note that the range of application of these equations in calculating 
oscillations in lock chambers, due to the smallness of the quantity 
(because of the slowness of the process of filling the chamber), may 
be somewhat extended, if we consider the celerity of the distribution 
of waves C to be a function slowly varying with time, for which COy<; 
should be considered to be a time variable: =f/h(t) /• 


Nevertheless, these equations will only be true till 
^Ux ^ (j^ ^00 „ ^ ^0- w 1 1 . ^ 


Cx) 


or 


■» Ux 


3^ ’ ■■ at "3 t ' 

where Q is the discharge of flow in the transverse cross-section 
of the chamber. 
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With relation to the distributional sluicing systems we may write down 
for the extreme cross-sections: 


Ux (0,t)= Ux (Ik, t) - 0 


5 


whence, with the aid of /l, 14 /, may be obtained 

(0,t) (iK.t) 


■^x 




= 0. 


/I. 19 / 


In the head supply system for the extreme cross-sections we may write 
down: 

a)K Ux (0,t) = a(t), Ux (1 kiU = 0, 

where Q (t) is the discharge of water delivered into the chamber. 


Hence, it follows: 

f \ { O't) _ I o' (f ) ^ Ik Q 


9 


/I. 20/ 


Let us denote the values ^ and on the section ( Ij) through 
*^1 and cOi , on the section (ij^, 1^) through and 0J2 * on the 
section (I2* Ij^) through ^ ^ ^ 3* Then the conditions, which must 

be obseryed in the intermediate sections x = Ij and x = I2 may be written 
down in the following way: 

(i==1o2.) /1. 21/ 

The first of these conditions follows from the Bernoulli equation, if it 
be written for cross-sections arranged at a small distance either to the 
left or to the right from the separating x = 1^^ and if velocity heads are 
neglected because of their smallness. The second condition may be 
obtained by the aid of /I. 4 ./ from the continuity equation written down 
for the same cross-sections. 

If a concentrated supply is realized in any one of the intermediate cross- 
sections of the chamber x = Ij, where the discharge of water Q.(t) is 
delivered, then this cross-section is subject to the following conditions: 




2 C)X 




Q.i 


/ 1 . 22 / 
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where'^ j and (Dj belong to the section to the left of the cross-section 
under consideration, while"^' ^^^and U)" ^belong to the section to the right. 


Since at the initial instant of time the fluid is at rest, the initial conditions 
are: 




n 


/I. 23/ 


When problems are solved for approach channels, boundary conditions 
are set up similarly. Let us note, however, that for the open end of 
the channel with a length 1, opening into a large, water basin, the 
boundary condition is / 2 24/ 

When local changes of the cross- section of the channel are confronted with, 
in particular, directly in front of the lock, then / 1,21/ must also be applied. 

Thus, the calculation of unsteady fluid movements _in locks and their approach 
channels in the indicated characteristical cases is brought down to linear 
wave problems of the simplest type, the mathematical apparatus for the 
solution of which has been well developed and is classical in essence 
(problems on the transverse vibrations of a string, longitudinal vibrations 
of a bar etc, ), 

In analysing fluid oscillations in locks it often becomes necessary to 
determine the distribution of longitudinal velocities along the length of 
the lock (for inst, , for the determination of the forces acting upon the 
ship). To that end, after the expression ^ (x, t) is established, equations 
/I, 4/, or /I, 12/ may be intergrated. 

Let us set up an expression for the wave force acting on the ship. In 
practical lock calculation, the wave force is known as the longitudinal 
horizontal component of forces of fluid pressure acting upon the ship: 


Pb- 


pcos (n,x) ds, 


Sc 


/L25/ 

where is the submerged surface of the hull, n is the direction of a 
line interior normal to the surface , The analysis shows that as the 
oscillations, both of the fluid surface and of the ship, are small the wave 
force may be determined integrating the formula / 1, 25/ for the 
submerged hull surface in a state of rest (l) (the error is of the 

second order of smallness). Considering, in addition, here and further 
on the submerged portion of the hull to be close to a rectangular 

(l) In the case of calculation of locking in a chamber, it is appropriate 
to consider the boat's water-line as coinciding with the mean position of 
the level in the chamber. 
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parallelepiped in shape, we will obtain: 





thence, utilising /1,3/, we get: 

Pb = li (l 2 ot]d^ = 


/I. 26 / 


Let us note that in the expressions presented, a reduced length of the 
ship is being considered: 


p = ^ 

00c 


•) 


where OJ is the area of the middle cross-section of the ship. 


/I. 27/ 


We will demonstrate the application of the expressions obtained here by 
a practical solution in the case of a head-^qn type of water suj^ly^ 


2, Calculation for the Initial Stage of Lock Chamber Filling with an 
End^on System of Supply, 

Let there be one boat in the chamber. In this case the level fluctuations 
are determined by integrating / 1. 14/ with q =. O, boundary conditions 
/1. 20/, /l,2l/ and initial conditions /1.23/, 

Practically it may be considered / l/, that at the initial stage of chamber 
filling 

Q*(t) ^ Q* (O) — constant, /2, 1/ 


Integrating by the Fourier method and presenting 
aid of dimensionless quantities, we get 


V a' CO) I 

V Cc > 


— ( 
Z 


the solution with the 

-if-l- X (x) - 


-2 



i_C05-^ Vn IXn m \ ^ 


72 . 2 / 
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where +A\n f , 


COS 


ai Q^x ^ 


X^W=J Ag^cos V^(x-{J+B2nSin a\ U^x^Iq 


■ cosS^n (x “?£) ^ BgpSin6^ (x-Ep) ^2 ^K» 


V 


For these same three sections (i=: 1,2, 3) : — 

Si (X) = S (Ti-,)- (x-U.,1 (1-i^ - 

l-K 


2ik r _ 

where '^(0)=2 {k 


oo 


n=i ’^n '^r> 


Coefficients 
A2n = C0S(P 


^30 = 


Szh — sin(D '))p 






sinV„ OOSVn 


R = — <3' 
•> '^30 


A2h 

COsVp Sin ))rt 


V 




The spectrum of values is determined by the transcendental 
frequence equation: 

<0^+ s' cotv (tarKoIiV + tan O' si) - 


- tanerl,')' tanSlaV 
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where values are connected with natural frequencies ^ n and 
corresponding periods of water oscillation in the chamber by 
the following relation: 



/Z. If 


Substituting from formula /2,2/ ^ into the expression / 1,26/ and 
utilizing / 1. 27/, we will obtain the following law for the action of 
the wave force in the initial locking period: 


p - g’ (0) 


yw 



i + 2 . 1 , 




+ 


OO 



in= 1 


A3h-'A2 


^ COS 





/ 2 . 8 / 


11, THE PROBLEM OF OSCILLATIONS IN THE BOAT - CARRYING 
TANK OF AN INCLINED SHIP ELEVATOR. 

Oscillations of the water surface and of the boat in the tank of ari 
inclined ship elevator may be investigated on the basis of two 
assumptions: 

1, the acceleration or stoppage of the tank occurs during 
some finite time interval; 

2, the moving tank stops in such a brief interval of time, 
which may be considered as instantaneous (sudden stoppage). 

The solution of the problem of an instantaneous chamber stoppage is of 
interest for the clarification of maximally possible undulations of the 
surface in emergency cases. It is treated in article / 5/ and therefore 
is not presented here. We wish to note, however, that for this case, 
when the assumed operating speeds of the tank are of the order of 
1 m,per sec, the undulations are quite admissible (this conclusion 
has later been confirmed experimentally). Now let us turn to the 
more general first problem. 
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3, Approximate Differential Equations of Water Surface 
Oscillations in the Chamber . 

In this case differential equations are derived in a manner similar to 
the method given above for lockst The complication consists in the 
following: in utilizing dynamic equations of fluid motion we have to 
regard as mass forces not only gravity, but also the inertia due to 
the transition of the system since it is expedient to examine this 
problem in a moving system of coordinates rigidly connected with 
the tank. 

With relation to the shape of the tank we will make the same assumptions 
as above in the case of the lock“chamber. We shall assume that the 
tank is moving along a straight rail“line, inclined at an angle ®. 

The motion of the fluid and of the boat in the chamber we shall refer 
to a moving system of coordinate xyz (Fig. 3). 

Since in the ship elevator tank may be permitted undulations only 
of a small height, the theory of long waves of small amplitude may 
be used as a basis for the solution of the problem in this case also. 

Now let us write down the dynamic equations for the fluid in motion; 




a X 

/3. 1/ 

1 

>< 

1 

c>X ~ 



( ~6p 

? ?)Z 

0 

/3.2/ 


where j^= j cos ©, _ j sin © , a j (t) - acceleration of 

the tank. The component of the acceleration is usually practically 
insignificant as compared with the acceleration of gravity g; therefore, 
we will disregard it in the equation /3, 2/, Then for a pressure p we 
will again obtain the formula / 1. 3/, 

Instead of / 1,4/ we will now have 



/3.3/ 


The continuity equation, derived with the same assumptions as for 
the lock, will result in 


, (Ux^) ^ 

"d t 


/3,4/ 
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Excluding from 1 3.31 and from /3,4/ we get (l) 



while at 


cO = constant 


73.5/ 


is 2 ^ 


/3.6/ 


where C, the celerity of wave distribution, is determined by the 
former expression / 1. 16/, / 1. 17/ and/ 1.18/. 


In the extreme cross-sections of the tank the bounchry conditions are 

~ 0 ) 

due- to which 




ri 


/3.7/ 


For intermediate cross-sections x«lj and x=lj by analogy with 
/ 1.21/ we now will obtain: 










/3.8/ 


1 £ the fluid in the chamber is in a state of relative rest in the initial 
instant of time t=0, then the initial conditions take on the shape 
/I. 23/. 


4. Forces Acting Upon the Boa t 

In regard to the boat placed in a moving boat-carrying tank, amongst 
/ 1 \ . 

Here g again an assumption is made that Ux is small. 
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the forces which must be taken into account in determining the 
tension in the mooring ropes, by which the ship is fastened to 
the tank are: 

1) the force of inertia (the transition inertia); 

2) the wave force (the wave pressure); 

3) the force of friction between the fluid and the ship 
(the surface resistance). As the analysis indicates, 
it is negligible. 

If we disregard as negligible the inertia of the fluid mass perturbed 
by the ship (the so-called virtual mass), then the inertial force may 
be considered equal to (l) 


= — — ^)rw. /4. 1/ 

The wave force is determined by formula /1.25/ and / 1 . Z 6/. The 
total longitudinal force acting upon the ship will be equal to 


P-Pu -t- Pb 

5, A Solution for the Simplest Cases . 

Let us present here the results obtained for two simplest cases: 

1) in the tank there either is no boat at all, or the boat 
is of such small dimensions that it does not influence 
to any considerable degree water fluctuations (cscj^)/ 

2) the length of the boat is close to the tank length 

, ^k)» and the cross-section of its submerged 

portion does not change throughout (C = C(;), 

In these cases the equation /3,6/ is applicable. Let us accept 
the initial conditions in the form of / 1. 23/, the boundary 
conditions - in the form of / 3. 7/, 

(1) Tests conducted by us on a ship elevator model have shown 
that in calculating the virtual mass may be omitted. 
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We shair considei' the motion to be uniformly accelerated ; j = constant. 
Integrating by the method of Fourier, we will obtain the following 
solution: n— ■( 




X - 


,Z n 


n- U 


. sin JUi — x' 


/5. 1/ 


where x* = X — — 

The second solution may be found by the method of D’Alembert: 


\=(x',t)-= 


tf' 


where 


hf(x'-ct)+f(x' + ct) 


/5.2/ 
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and so forth, with either the increase or the decrease of . 
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It follows from formula /5. 1/ and /5.2/ that the period of oscillations 
and the duration of the oscillation phase T" are determined by 
the relation ° 


'y' To ik . / n / 

Lo — C — ^ 

however, when there is a boat in the tank then evidently the following 
relation will be more precise: 


To=-|^ = -^+-^r-^ /5.4/ 

^ Ck 

For indivj.dual time instants the longitudinal profile of the free surface 
is presented on Fig, 4: 

0- 0 - at instants t=0 and T , 

1- 1 - at instants^ To T^, 

2 - 2 - at instants -T— T and - 2 _t 

4 o 4 . o, 

3- 3 - at instants -^T and ■ T 

8 o do, 

4- 4 - at the instant — ^ — T . 

2 o 

Thus, in the process of uniformly accelerated motion of the chamber 
the fluid surface oscillates about the position (profile 2-2 on Fig. 4), 
which it occupies in the case of relative rest. In that position the 
surface is a plane, inclined to the horizon under an angle ip , 
where 


tan 'P = — 




/5.5/ 


In the extreme positions the free surface also has the form of planes - 
a horizontal and an inclined, under an angle tp to the tangent 


tan 




max 


3 


/5.6/ 
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The extreme position of the surface 4-4 (Fig. 4) will be reached, and 
in that connection the formula (5, 6) will also have a meaning only in 
the case when the time of acceleration 


T> ^ 


15.71 


Due to the smallness of accelerations, which may be conveyed to 
the boat carrying tank under normal operating conditions, because 
of the need of restricting the tension of the mooring ropes, the 
value of the acceleration period T will be quite considerable; there- 
fore, T will meet the condition 1 5.7 1 . 


Correspondingly the greatest oscillations of the surface 
front and rear tank gates will be equal to 


ai Ducn Liie 




max 




JdL 


tan4^^ax=+-^lK 


/5. 8/ 


By expressing j through the time T, in which the acceleration or 
stoppage of the tank takes place, and through the speed of uniform motion 
V, we will get 

- vl cos e 

+ /5.9/ 




max 


The maximum value of the wave force will be reached at such 
instants of time, which correspond to the odd number of phases 
of oscillation T , when the free surface attains the maximum incline, 
determined by the formula /5,6/, Here it follows from / 1,26/ and 
/ 5. 8/ that 


b max 




max 


t^c=-2-4^^w. 

a 


/5. 10/ 


At the beginning of the period of acceleration (or instant corresponding 
to the even number of phases ^q) the longitudinal force acting on the 
ship is reduced to the single inertial component 

P=rP =- ^ r W. /5. 11/ 

6 

In the process of acceleration at the instant t = (or the corresponding 

odd number of phases HC 


P = P + p, 

u ^ -^b max 
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But if the accelerated motion of the chamber ceases at such an instant, 
then the action of the wave force is displayed fully: 




max 


2 


i^irw 


/5. 13/ 


Of all the possible magnitudes of the longitudinal force this one is the 
greatest. 

The allowable longitudinal tension of mooring ropes, by which the boat 
is fastened to the tank, may be determined for normal operating 
conditions with the aid of A. V. Mikhaylov’s formula, which he received 
on the basis of an analysis of the conditions of ship stay in lock chambers: 

p = y , 

allow. Hallow. 

for which Mikhaylov gives both a formula and a diagram for the yalue 
^allow. • the value of allowable tensioning 

of ropes was lessened by introducing the coefficient ^ 2 ^ 0*9* which 
takes into account the increase of the capacity of the lock-filling system 
as compared with the data of their model tests (reaching in some cases 
10 %). 

Since this has no relation to ship elevator operation the coefficient K^, 
must be excluded from the results. 


The allowable (in calculations and model investigations) longitudinal 
tension of mooring ropes in application to ship-elevator conditions 
may be presented in the form 


P 


allow. 


y w 

rt * 

' allow. 


where 


/5. 14/ 


n' 

allow. 


^ allow. — 0.9 n 

allow. 


If we disregard the damping of oscillations and assume that the tank 
acceleration ends at the instant, which corresponds to the odd number 
of phases , then equating expressions /5. 13/ and /5. 14/, we 

will receive tlie maximum allowable normal working acceleration 


) 


allow. 


1-8 


n 


9 

COS0 

allow. 


/5. 15/ 
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6 . A More General Solution 


K the assumption 1 1 is not introduced in the presence of a boat 

c k 

in the tank, then the equation /3.5/ should be integrated for boundary- 
conditions /3,7/ and 13.8/ and initial conditions /1.23/. Considering 
that the acceleration is constant (j rr constant) and integrating by the 
G. A. Grinberg method, we will get for t ^ T: 




3 s- 


to 


71=1 


-dll^(t-C054sV„t)X„(xW6.1/ 

V Nn U 


or 


at t ^ T: 


cjx 9 

% ^ 




OO 






4<ix 

S (o' 


cP 


n = 1 


Here 


'iP- Sin 


T T7 - t) (^ )• 

/b.i/ 




while N,^ , X (x) and A. > B. are determined by the previous 
n in in 

expressions. Natural frequencies are also subjected to the previous 
expressions /2. 6 / and/2,7/ . 


Applying / 1.26/ and / 1,27/ , we will get the following law of action of 
the wave and total forces. For t < T the wave force: 


P =: ^ ^ 

b Q (32 


Kn (A2n-A3n l( /|^ cOS-kyj:), 

- Nn 

'Tipi /6.4/ 


or 


r 


p = 
b 




1 

2 


(52 


f>0 


n =1 


Kn(A2h Aatj 




/6.5/ 
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the total force along /4, 1/ , /4.2/ and/6,5/ is 


P=Pu^Pb = - 


2 jx £W 
? 02 


oo 




KnW^n'Aarlffiq Cc \) i 

y * 

Vn Nn ic 


/ 6. 6/ 


At t > T: 


P=Pb= 


4^x 

9 


oo 


rw 

6:^ 



KnlAan^ Asni) 





/6.7/ 


7» Brief Data Relating to Experimental Investigations On a Model 
of an Inclined Ship Elevator, 


In the years 1957 - 1958 in the Moscow Kouybishev Institute of Civil 
Engineering the author was in charge of experimental investigations 
on a model of an inclined ship elevator, A self-propelled boat- 
carrying tank model running along a section of the railroad was 
built to a scale 1 : 50, Here the data of the prototype were assumed 

to be: 1 ^ = 110 m, bj^^^lSm, h=3,65m, tan©=0,05, the speed of 

uniform motion ranging fr^m 0, 8 to 1,2 m/ sec, , acceleration in the 
range of 0, 2 to 4 cm/ sec. • Instantaneous (practically, very rapid) 
chamber stoppages were also effected on the model. The boat was 
represented by model of a barge of 3, 330 't.displacement, and 1^= 85m, 

b^=14m.and s=3,2 m. Modelling was conducted in accordance with the 
laws of gravitational similarity. 


In the course of the tests synchronous records were made on an oscillograph - 
of the speed of chamber motion (with the aid of a tachometer), of water surface 
oscillations (with the aid of an electrolytic resistance wave meters), of 
variation of the longitudinal force acting on the ship (with the aid of a inductive 
measurer). 


The test data indicates that the theoretical calculations of the water surface 
oscillations in the tank and of the longitudinal force acting upon the boat 
are practically exact. 
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DENOMINATIONS 


k - chamber length, 

- chamber width (or approach channel width), 
along the free water surface, 
h - depth of water in the chamber (or channel), 

U)k - area of the cross-section of the fluid in the free 
chamber (or channel), 

“ length of ship, 

^c - width of ship along the water line, 
s “ ship settlement. 


tJc - cross-sectional area of the submerged portion 
of the ship, 

W - Voluminal ship displacement, 

p and p^ - pressure in the fluid and at its free surface 


- acceleration of gravity, 

and iS ~ mass density and specific weight of water, 

- time, 

volume rate of flow entering the lock chamber (or channel) 
volume rate of flow entering the lock chamber (or channel) 
through the distributional system (the longitudinal culvert) 
per unit of length, 

“ incline angle of the shipways, 

- running speed of the boat-carrying tank (chamber), 

- acceleration of tank motion, 

► 2, ^ horizontal and vertical components of the 

acceleration of tank motion, 
ordinate of the free fluid surface, 
cross-sectional area of fluid, 

-wave force acting upon the ship, 

« force of inertia acting upon the ship (in the ship-carrying tank), 
j C 

■'k and c - celerity of wave diffusion in the free sections 
of the chamber and in the section occupied by 
a boat respectively, 

ic " t ; ' 


g 

? 

t 

Q 

q 

e • 

V 

j • 

jx. 

V 

ti- 

p 

b 

p 

u 


- X * r — 9 


Ic 


Ic 


1 . 




dimensionless parameters. 
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A STUDY OF ECONOMICAL NAVIGATION LOCKS 


25X1 


FOR HIGH LIFTS 


* * >l< 


WALENTY JAROCKI 

s.'c 

SUMMARY 


This paper describes the investigations on pneumatic chamber 
locks, requiring small amounts of water for sluicing and utilizing 
compressed air to raise the tank. This lock is designed for use 
on those reaches of navigable canals where a water deficiency 
appears, as well as in cases of considerable lift heights. 

In order to learn the principles and the functioning of pneumatic 
chamber locks , computations were made for and laboratory 
investigations carried out on a model of a mobile tank in a chamber 
using compressed air for raising the tank. 

The analytical study of the working conditions in a pneumatic 
chamber lock was based on the dynamic equilibrium condition 
resulting from d*Alembert*s principle. According to this prin^- 
ciple, a dynamic equilibrium equation between engaged forces 
was established. Having solved this equation, the relations were 
found which permitted the computation of air pressure in the 
chamber required to raise the tank, etc. 

For laboratory investigations a model chamber and a mobile tank 
with all the connecting devices was built. Particular emphasis 
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was paid to the problem of providing a suitable tightening in order 
to reduce compressed air losses and to prevent the seepage of 
water into the chamber containing air. This would occur near 
walls where the air is in contact with water. 

On the basis of laboratory model tests, diagrams were plotted 
showing the relation between time, tank motion, velocity and 
the pressure of the compressed air in the chamber for various 
water depths in the tank and for different values of tightening 
pressure. 


SOMMAIRE 


Ce m6moire d6crit l^etude faite sur des 4cluses pneumatiques , 
Celles-ci ne n4cessitent qu’une petite quantite d*eau pour I’^clusage, 
utilisant de Pair comprime pour elever le reservoir, Ces Refuses 
sont particulibrement utiles dans les biefs oh Peau n*est pas 
abondante ou aux endroits oh la hauteur de chute est considerable, 

Afin de mieux saisir les principes et le fonctionnement des ecluses 
pneumatiques, des calculs ainsi que des experiences furent 
effectues sur un modble reduit de reservoir mobile et un sas 
contenant de Pair comprime necessaire h elever ce reservoir. 

Les calculs d^analyse de conditions de fonctionnement des ecluses 
pneumatiques ont 6te effectu6s h Pappui de la condition d*6quilibre 
dynamique s’ensuivant du principe de d*Alembert. Selon ce principe, 
une Equation d’4quilibre dynamique entre les forces agissantes a 
6t6 6tablie, Cette equation 6tant resolue, on a trouv6 les relations 
permettant de calculer la pression d*air necessaire dans le sas pour 
eiever le reservoir, etc. 

Pour effectuer les essais de laboratoire, un modhle du sas et du 
reservoir mobile fut construit. On a prSte une attention toute 
particulihre h Petanchement qui doit prevenir les fuites d*air et 
egalement la penetration de l^eau dans le sas. Ceci se produirait 
prfes des murs oh l*air est en contact avec Peau, 

D’aprfes les experiences de laboratoire faites sur le module, des 
graphiques furent traces montrant la relation entre le temps, la 
vitessede inouvement du reservoir et la force de pression de Pair 
comprime dans le sas, pour differents niveaux de remplissage du 
reservoir et pour differentes pressions d^etancheite. 
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1. INTRODUCTION 


On some navigation channels, where water is lacking, or where the lift 
is great, the construction of navigation locks requiring small water 
amounts for sluicing is desirable. Included with this type are the pneu- 
matic chamber locks, which seem to have the most advantageous con- 
struction for the above mentioned conditions. 

Pneumatic chamber locks may be divided into 2 types: first - those 
that work on the pneumatic clock principle, and second - actual pneu- 
matic locks (fig. 1) . The first type of pneumatic chamber lock (Dulton, 
Morawek) is rather expensive - a strong foundation as deep as the fall 
height being required. It seems, however, that the second type, one 
designed by Schancer (fig. lb) and another by Tillinger (fig. l,c,d) 
is more suitable. 

The latter type consists of two chambers filled with compressed air, 

A moving tank filled with water is placed in each chamber. When any 
vessel enters the tank a certain amount of water is added and the valve 
is opened. As the air passes from one chamber to the other, one of 
the two tanks rises, while the other falls. 

While the tank is in its upper position, the chamber walls are under 
relatively small air pressure, approximately equal to the water 
pressure of 3,5 m. of depth, and does not increase with the increase 
in fall height. In normal chamber locks, however, pressure increase 
is proportional to depth. That means that pressure in pneumatic 
chamber locks is less than in normal chamber locks which require 
a greater wall thickness. 

It was reckoned that in pneumatic chamber locks the sluicing pro- 
cedure requires about 0.3 m. depth of water, independent of the height. 

So with a 10 m. fall height, and using economical basins, the amount 
of water required by a pneumatic chamber lock will be about 6% of 
that used by a normal chamber lock. With a fall as high as 12 m. it 
decreases to 5%, and to 3% in case of a 20 m, fall height. 

The Schancer pneumatic chamber lock has a moving water-filled tank, 
while in Tillinger*s system this tank is turned upside down (fig, 1) . 

In the second case the walls may be high, and tightening between tank 
and chamber walls ought to be placed as low as possible, Tillinger 
affirms that the pressure exerted by a column of water on the tighten- 
ing must exceed the compressed air pressure in the chamber, other- 
wise air can escape from the chamber » 
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2. ANALYSIS OF THE WORKING OF THE PNEUMATIC CHAMBER 
LOCK 

In order to obtain an analytical picture of the conditions under which 
pneumatic chamber locks work, the movements of the tank in a chamber 
filled with compressed air were observed. 

As a starting point of this discussion we will consider equilibrium as 
the condition not only the equilibrium of weight of the tank filled with 
water remaining under compressed air pressure, but also the setting 
of this bed in motion. 

The condition of dynamic equilibrium resulting from d'Alambert's 
principle is a differential equation of the second degree, motion 
acceleration being the second derivative of distance against time. 

The equation of dynamic equilibrium between acting forces will be as 
follows : 

- Q - FX'H - fd + kF = Q -F F y H d^Y 

g dt^ 

where: 

Q « weight of tank, 

H = depth of water column, 
i « unit weight of water, 

F ^ a cross section through the chamber on the level 
of the tank bottom, 
d = pressure in tightening, 
f « friction coefficient, 
k = compressed air pressure, 
g a acceleration due to gravity, 
y « stroke or travel of the tank, 
t a* time of tank motion. 

In this equation the negative factors represent forces acting downwards, 
while the kF factor, the pressure force directed upwards, has the plus 
sign. 

The dynamic equilibrium equation after transformation may be written 
in the following form: 

d^y _ ( Q -F F y H -F ' id - kF) 

dt^ ” Q -F F y H 




-fd -F kF 
Q+Fy H 


) 


g 


l-B-4 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 



Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


Assuming 


A 


- 14 - 


-fd 4- kF 
Q + F^ H 


( 1 ) 


we shall obtain a differential formula of tank motion: 



Ag 


( 2 ) 


In this equation an unknown factor is the travel of the tank in function 
of time i.e, y = f(t) . After obtaining the solution of the differential 
equation, we determine the distance in function of time, according to 
boundary and initial conditions. It may be noted, that the initial 
conditions are relative to that moment when the movement begins , 
while the boundary conditions are the beginning and the end of the 
stroke. 


The integral of the differential equation of the motion, has the follow- 
ing form: 


y = 


Agt2 

2 


+ B t + C 


Thus, in order to fix two optional constants B and C, it is assumed 
that y (0)= O, and y* (O) — O, that means, that in the initial moment, 
distance and velocity are equal to zero. With such an assumption 
the tank movement will begin without any jerk at the initial moment 
(t = O) . 


Under such initial conditions the form of the integrate is 

v^ Agt2 
2 


(3) 


Considering that Tq is the time required for the rise of the tank to the 
entire height Lq* with the value of A corresponding to the known para- 
meters Q,'}!' , H, d, k, F. Substituting those values in the equation (3) 
we obtain: 


Lo- 


^ Ag 


= 


A 

T 


(- 


1 + 


k F - f d 

<a. + 


■) 


k F« f d 




(Q 4- F H) 


(4) 
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This final equation includes all parameters. Let us assume the value 
we try to determine is the pressure k required for the elevation of the 
tank under given conditions. The following formula thus results: 


k 



d 



1 + 


2 Lo \ 



(4a) 


As we are not interested in the analysis of the phenomenon of tank 
motion expressing the relation between stroke and time, therefore, 
developing the equation (4a) parameter y/gt^ was written in the Lg/g Tq^ 
form. An essential problem is the determination of the pressure of 
raising the tank in T^ period required for reaching the final height Lq. 


Dividing equation (4a) by (5. j-ltH) we shall obtain the relation: 

F 


K 


F 




)r H ^ 


7 

L O 


_Q_ 

F 


+ 


7 1 


(5) 


The above equation, if compared with equation (4) does not give any new 
solution, as it is synonymous. Equation (5) , however, is of a great 
value for investigators as it shows which dimensionless parameters are 
governing the tank motion phenomenon in the chamber, and which are 
the similarity parameters in model tests. 


The dimensionless parameter 


Q 

F 


_k 

-F 


represents the relation 


between air pressure K in the chamber and the pressure exerted by the 
rising chest. The parameter ^ is the quotient of tightening 

+ JTH 

pressure to the pressure of the same kind as in above mentioned para- 
meter, This parameter includes the numeral coefficient f:F. The third 
parameter is connected with the tank motion phenomenon. 


When the parameters ^d;=: constant and -I- H =r constant 

F F 

the formula (4a) can be expressed by the hyperbola equation 

Lq 

k- C1 + C2 — ^ 


This equation confirms that if the time required for the entire tank rise, 
Tq, is not constant, then with the increase of the duration phenomenon, 
the air pressure maintaining the motion will decrease. 
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It is evident (eq, 3) that the tank is moving with uniformly accelerated 
motion without initial motion, A rapid braking might be extremely 
undesired after the tank reaches the height Lq, Previous braking should 
be applied thus reducing the air pressure in the chamber. 

Let us assume the uniformly accelerated tank rising motion in 0,9 
Tq period, with v = y' — A g t velocity, is followed by braking procedure 
within 0,2 To period. This may be made by reducing the compressed 
air pressure in the chamber. 

Assuming that the second phase of motion (retarded motion) from 0,9 
Tq moment, the motion velocity will decrease linearly so that after 
0,2 Tq period it will be equal to zero (fig, 2) , Under such conditions 
the duration of the tank motion will be 1, 1 Tq of the period adopted. 

During 0,9 Tq period (fig, 3) the tank will travel the following distance: 
^1= g A (0,9 To)2 = gATQ2 (6) 

The remaining distance is: ' 

Lq - Li«_^i9_ Ag Tq2 (7) 

In order to obtain the zero velocity of the uniformly retarded motion 

at the end of ^ A g Tq^ distance, when the starting velocity was 
2 

0,9 Ag To» the question arises of adopting a suitable acceleration value 

y" = “AQg. 

The retarded motion equation has the following form: 

y = - Aq g t^+ Bt - 4 - C (8) 

Since the motion velocity and acceleration are: 

y' = - Aog t + B (9) 

y"- - Aog (10) 

The following values might be applied in order to determine the 3 constant 
parameters (Aq, B, C) in motion equation (8); 

a) The coordinates of point I (eq, 8) 

(t. 0.9 To, AgTo^). 
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b) The coordinates of point II 

(t = 1,1 T„. AgTo^) , 

c) The coordinates of point III 

(t=-0.9 To. 0.9 Ag T^) 

By substituting in equation (9) the coordinates bf the point III, we obtain 
the relation: 

- 0.9 Ao g To +B = 0,9 A g To (H) 

The coordinates of the points I and II might be substituted in the equation 
of stroke (8) : 

- Ao g To^ + 0,9 B To -I- C _ A g To^ (12a) 

2 2 

- Ao g T„2 1, 1 B To + C - Lo = A g 1^2 

Subtracting by sides the equations (12a) and (12b) we eliminate C and 
obtain the following relation: 

- (1.i2-0.9^) Aq g To^-hO. 2 BTp =r Lq - AgTp^. 

- 0,2 A^ g T^q- 0.2 B=_2^ AgTo. 

The above equation forms with the equation (11) the following system: 

- 0,2 Ap g T^^. 0.2 B , AgTo; 

2 

- 0.9 AogTo4-B= 0.9 AgTo 

Multiplying the upper equation by (-5), and adding it to the second one, 
we eliminate the unknown factor B: 


Since the value Ao might be computed: 

A^=.4,25A (13) 
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In order to evaluate the B constant, we substitute the value of Aq in 
equation (11) obtaining: 

3 = 0,9(4,25-1-1) A g To = 4,725 A g To (14) 

By substituting Aq and B in equation (12b) we obtain the C value; 

C=AgTo^ (2,5713-5,1975+0,5)= -2, 126 A g To^ (15) 

Subsequently the constants A^, B and C are substituted in the retarded 
motion equation (8) thus: 

y =r - ^ ^ — + 4,725 A g To t - 2, 126 A g 

The above obtained retarded motion equation is adopted in such a form, 
that in the moment t=? 0,9 Tq when the braking begins, the tank will be 
in point It will stop on the level when t = 1, 1 

Let us take, for verification, the derivative of the tank motion velocity: 
y'= - 4,25 A g t -+ 4,725 A g To 

This velocity value must be equal to zero, when t ^ 1, 1 T^* Substituting, 
we obtain: 

- 4,25 X 1,1 A g Tq 4-4,725 A g T^ = (4,725-4,675) A g T^ - O. 

3, THE PRELIMINARY STUDIES ON TIGHTENING AND TANK MOTION 
IN THE CHAMBER 


In order to carry on some laboratory experiments on pneumatic chamber 
locks, a model chamber was constructed* It has the form of a reinforced 
concrete tank 100 cm. high. The general appearance of the devices applied 
for tank motion tests in the chamber is shown on the Fig. 4. A3 mm, 
thick, metal mobile tank was stiffened with angle plates. The compressed 
air was conducted to the chamber. 

As it is difficult to obtain suitable tightness between water and compressed 
air, a special attempt was made to the problem of tightening between mobile 
tank walls and chamber walls. Several types of tightening with a rubber 
pipe filled with water were designed and tested. The water pressure in 
this pipe was varying with the variation of water-level in a water reservoir 
joined to the water-supply network and the rubby tightening pipe. This pipe 
was settled in a steel sheath welded to the tank, with an outlet along it 1/3 
circuit wide. 
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It was found that the most suitable tightening is a rubber pipe 2,6 cm. of 
diameter settled in sheath covered with leather tape retained by thin 
leather straps fastened to the borders of the sheath (fig. 5) . This tight- 
ening was applied in the preliminary investigations consisting of 45 tests. 

In general the whole device was working in a quite satisfying manner, 
except some cases when the tank moved headlong. 

It was confirmed, according to the preliminary investigations, that the 
projected tightening fulfilled the designed purposes, and the other elements 
worked well. Only the chamber depth was found to be insufficient for 
certain research proceedings. 


4. SOME DETAILED INVESTIGATIONS ON WATER- FILLED TANK 
MOTION IN AN AIR-CHAMBER 


The detailed investigations began after a supplementary chamber of 0,3 m. 
height was constructed, and after the rubber pipe filled with water was 
applied as the tightening. 

The investigations involved the following measurements: 

H - depth of water in the tank (cm.) , 
d - water pressure in the tightening pipe (m.H^O) , 
p - pressure before the Venturi meter (mm^Hg) , 
h - pressure on the Venturi meter (mm.Hg) , 
k - pressure in the chamber (cm. H^O) , 
z - air-pressure on the manometer (atm.) , 

L - differences of levels between the upper and the lower 
tank positions (cm.) , 

T^- tank rising time (min.) , 

T 2 - tank falling time (min.) , 

t - air temperature (C°) , 

b - barometric pressure (mm.Hg) , 

The main purpose of our investigations was to establish a relation between 
water depth in the tank, air pressure in the chamber, water pressure 
in the tightening pipe, and tank motion in an air chamber. 

The investigations were conducted with the temperature of the environ- 
ment varying from 16,5® to 24®C. , and barometric pressure being about 
1 atm. 

The height of the pressure applied in our investigations varied from 
0,5 to 4,0 m, of water column in the tightening pipe, while the depth 
of water in the channel was from 5 to 25 cm. (every 5 cm. subsequently) . 
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Pressure p on the conduit was changed within 5 to 80 cm. of Hg (every 
5 cm, approximately) , There was 40 cm, of difference between the 
upper and the lower tank positions. 

A negligible air flight was observed, and it was found that a certain 
minimal amount of water was leaking from the tank to the chamber, 

A total of 23 test series were made, which included 338 measurements. 
According to data thus obtained, the = f/p/ relation was shown by 
23 diagrams. These diagrams (which are not presented in this paper) 
were used to procure 5 synthetical diagrams Tj = f(p) for 5 measurement 
groups. Every diagram considered measurements made under the same 
water depth in the tank, and involved all pressure values applied in the 
tightening pipe under these conditions. On the basis of results above 
obtained, 5 synthetical diagrams v=f(p) were deduced. Fig. 6 and 7 
show diagrams of time and rising velocity for a tank filled with water 
to H 15 cm, and H ~ 20 cm, of depth. 

These diagrams prove that when the time of tank rise increases, the 
pressure in the chamber required for the existence of the motion 
phenomenon decrea-ses, and vice versa. Equation (4) leads to the 
same conclusion. 

Diagrams showing the variations of air pressure, K, in the chamber 
required for different pressure values in the tightening pipe are drawn 
on Fig. 6-7. The lines obtained on these diagrams are generally 
parallel, this solution being in conformity with equation (1) . These 
diagrams suggest that, as the pressure in the tightening pipe increases, 
a greater air pressure in the chamber might be necessary in order to 
raise the tank (eq,4) , 

Fig, 8 shows the relation curves for Tj = f (p) and v = f (p) plotted for 
different amounts of water in the tank. According to these curves, 
and to the equation (4) it is evident that a greater pressure is required 
for raising a heavier tank up within the same time. 

In table 1 are collected the values of dimensionless parameters, k:c, 
resulting from the many computations arising from the different 
pressure, d, in the tightening pipe and for different weights of tank 
filled with water -c. According to the above mentioned table, a diagram 
was plotted (Fig. 9). It shows that as water pressure in the tightening 
pipe increases, the value of the dimensionless parameter, k:c, also 
increases; thus a higher air pressure might be applied in the chamber. 

Fig. 9 shows an agreement between experimental records and an analytical 
solution expressed by equation (5) , According to Fig, 5, it is evident 
that if the pressure in the tightening pipe decreases to zero, the value 
of the k:c ratio will be greater than one; the equation (5) implies the 
same solution. If in this equation (which is the straight line equation) an 
assumption is made that d= 0, the k:c ratio will become equal to the 
two remaining components, the last being equal to one. 
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TABIE 1. 


Coiig)utation of the dimensionless parameters 
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Pressure 
d in the 
tighten- 
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nnmraniinraiiinimin 

Depth 
of water 
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mrammmmmmmmmi 

Pressure 
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Air pres- 
sure in 
the cham- 
ber 
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m 
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H cm 
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p 
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p 

kg/cm 
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1 1 

0,5 

10 

0,0257 

0,0312 

1,52 1 

1 ^ 

0,5 

15 

0,0287 

0,0583 

1,55 I 

1 ^ 

0,5 

20 

0,0357 

0,0469 

1.59 i 

I ^ 

0,5 

25 

0,0587 


1 

1 




The average 

1,55 1 

1 5 

1,0 

10 

0,0257 

0,0342 

1,44 1 

1 ^ 

1,0 

15 

0,0287 

0,0403 

1,40 1 

1 

1,0 

20 

0,0357 

0,0467 

1,40 1 

I 8 

1,0 

25 

0,0587 

_ 0.0555 

1.58 1 

1 



The average 

1,40 1 

1 ^ 

1,5 

10 

0,0237 

0,0375 

1.58 1 

1 

1,5 

15 

0,0287 

0,0470 

1.64 1 

1 

1,5 

20 

0,0357 

0,0509 

1.51 1 

1 

1,5 

25 

0,0387 
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1.48 1 

1 




The average 

1.55 1 
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10 

0,0257 

0,0399 

1,68 1 

1 

2,0 

15 

0,0287 

0,0460 

1,60 1 
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2,0 

20 

0,0557 

0,0549 

1,65 1 

1 

2,0 

25 
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1 




The average 

1,64 1 
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5,0 

10 
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1 
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5,0 

15 

0,0287 


1 

1 

5,0 

20 

i 0,0337 

0,0620 

1,85 I 

1 

5,0 

25 
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0.0780 

2.01 1 





The average 

1,92 1 

LmJ 

nmmminmmmmmti 

inmmmmmmmmm 
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5. CONCLUSIONS 

As a conclusion to the analytical computations and the laboratory exper- 
iments performed on a model device utilizing the principle of a pneumatic 
chamber lock, the following remarks may be stated: 

a) There is full agreement between the analytical solution and the labor- 
atory model experimental results. 

b) As the total tank rising time increases the compressed air pressure 
maintaining the movement decreases. 

c) Tank rising velocity rises with an increase of air pressure in the 
chamber. \ 

d) Increasing the water pressure in the tightening pipe, requires 
greater air pressure in the chamber in order to raise the tank at the 
same time. 

e) Tank motion in the chamber remains smooth (without shocks) , 
providing that the position of the mobile tank is absolutely horizontal. 

f) As the tank rising in the chamber proceeds with a uniformly 
accelerated motion adequate braking should be applied. 

g) Required tightness can be obtained at the confluence of water and 
compressed air. 

h) Established formulas afford the possibility of the calculation of the 
required factors such as friction coefficient, compressed air pressure 
in the chamber lock in the initial and final motion periods, etc. 

i) Because of the small model size, it was impossible to state the 
advantage of Tillinger*s method, which proposes the solution of the 
tighting problem on the lower borders of an upside down turned tank. 

For the same reason we were not able to determine the condensed air 
losses in relation to the tonnage, time of sluicing and temperature of 
the environment. 

j) It appears from what studies we have made, that the application of 
the pneumatic chamber locks is quite possible and economic. It is 
probable, however, that in practice, several difficulties will arise 

at the time when details of construction of such a sluice are being 
worked out. 
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Figure 1 - Pneumatic chamber locks schemes. 
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Figure 2 - The diagram of tank motion velocity. 
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Figure 4 - The cross-section through the chamber and the tank, 
and the general appearance of the devices. 
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Figure 9 - The relation between pressure d in the tightening and 
the dimension less parameter. 
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SUMMARY 


The paper presents a brief description of laboratory investi- 
gations of the total hydrodynamic pulsating load on bottom 
outlet sliding gates. The problem of scale-modelling of this 
phenomenon is also considered. 

The investigations showed that pulsations are unperiodical 
oscillations recorded as a number of positive and negative 
impulses different in form, value and duration. The oscil- 
lations are considerably dependent on head pressure and 
hydraulic conditions as well as on the vortex zone forming 
behind the gate. The value of maximum amplitude of the 
total hydrodynamic pulsating load pertaining to the average 
hydrodynamic loads can reach 55-60% and higher. For the 
investigated scheme of gate location the value of pulsation 
under study reached the highest magnitude with full gate 
opening (n = 1,0), When the lower edge of the gate is raised 
above the top of the conduit, i, e. when n > 1. 0 the magni- 
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tude of pulsation first sharply increases reaching the 
maximum at n = 1. 13 - 1.14, then gradually decreases 
and finally drops to zero at n = 1. 35 “ 1. 40. 

The paper presents the universal design diagrams for 
determination of the total hydrodynamic pulsating load. 

Investigations of scale -modelling of the total hydro- 
dynamic pulsating load upon sliding gates has been 
carried out on models built to different scales (the 

scale series). Th? investigations showed that the ^ 

basic principles of this theory relative to the auto- » 

model region completely proved out. Within the 
range of Reynolds numbers adopted in the experi- 
ments with the three scale models (2760 ^ 80000) 

in the case of transition from one model to another 
the absolute magnitude of the maximum amplitude 
of pulsation varies in proportion to the cube of linear 
scale coefficient while the dominant pulsation frequen- 
cies vary inversely proportional to the square root of 
the linear scale coefficient. 


SOMMAIRE 


Le present rapport comporte une breve description 
des resultats des recherches exp6rimentales sur la 
charge de pulsation sommaire hydrodynamique sur 
les vannes Sl glissiere de fond. Le probleme des 
etudes sur les modules est aussi traite, 

Les recherches ont demontre que les pulsations 
6tudi6es representent les oscillations desordonn^es, 
enregistr6es comme un nombre des impulsions 
positives et negatives, differentes en leurs formes, 
valeurs et duree, 

Ces oscillations dependent consid6rablement 
de la chute, ainsi que du regime hydraulique de 
comportement de l*ouvrage et de la position de 
zone tourbillonaire qui se forme derrifere la vanne, 
En rapport aux charges hydro dynamiques moyennes 
la valeur de l*amplitude maximum de la charge de 
pulsation hydrodynamique sommaire sur les vannes 
cL glissifere peut atteindre 35-60% et m^me da van- 
tage. Pour le scheme de la position de vanne 
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6tudi6 le valeur de pulsation atteint son maximum cL 
I'Ouverture totale de la vanne (n = 1.0). Quand le bord 
inferieur de la vanne est mont6 au-dessus de la paroi 
sup4rieure de la conduite d'eau, c. -k. -d. quand n > 1 .0, 
la valeur de pulsation au d6but s'accrott brusquement et 
aun = 1.13 - 1.14 atteint la grandeur maximale, aprfes 
quoi elle commence se reduire et au n = 1 . 35 - 1 . 40 
elle baisse jusqu'k z6ro. 

Dans le rapport sont pr6sent6s les graphiques universels 
de calcul pour la determination de pulsation de la charge 
sommaire hydrodynamique. 

Les etudes sur modules reduits de la charge sommaire 
hydrodynamique sur les vannes ^ glissibre de fond ont 
ete executees sur plusieurs modules avec differentes 
echelles (la serie d'echelles). 

Les etudes ont prouve que pour le phenombne etudie les 
principes fondamentaux de la theorie d' experimentation 
sur les modules concernant le domaine d'automodble, 
sont entiferement confirmes. Dans les limites des valeurs 
des nombres de Reynolds existant dans les trois modules 
de la serie d'echelles (2760 ^ Re 80000) en passant d'un 
module k I'autre, la valeur absolue de I'amplitude maxi- 
mum de pulsation etudiee varie proportionnellement au cube 
du coefficient lineaire d'echelle, tandis que les frequences 
de pulsation predominantes varient en proportion inverse de 
la racine carree du coefficient lineaire d'echelle. 


1. INTRODUCTION 

The experience gained from hydraulic structures now in service as 
well as the data obtained from model and field investigations show 
that gates are subject to a vibration which is due to the pulsation of 
hydrodynamic loads acting upon them . 

Under conditions of ever increasing head pressures and diameters 
of outlets closed by gates, the problem of vibration appears to be 
of considerable importance. 

However, up to the present moment, the state of development of 
theoretical and experimental research of gate vibration does not 
meet the requirements established by common practice of design 
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and operation of hydraulic structures. In this connection thorough and 
systematic studies of gate vibration and development of proper measures 
to prevent or eliminate this phenomenon are of great theoretical and 
practical interest. 


2, ESSENCE OF THE PROBLEM IN QUESTION 
AND ITS PRESENTATION 

As it has already been pointed out in our paper gate vibration in 
field conditions is an extremely complicated phenomenon which, 
other factors being equal, depends, as to acting pulsating loads, upon 
the structure and rigidity of an oscillating system, supports and sup- 
porting parts as well as upon external and internal forces of non-elastic 
resistance. 

Theoretical solution of the problem and reliable scale-modelling of the 
vibration phenomenon with observance of similarity of basic determin- 
ing parameters is extremely difficult and apparently cannot be obtained 
practically. 

Because of this, when studying the vibration phenomenon in laboratory 
conditions, the problem is divided in two stages: 

a) experimental, comprising investigations of amplitude and pulsation 
frequency of total hydrodynamic load (exciting force) for different 
structures and different service conditions of deep gates; 

b) theoretical, when main parameters characterizing gate vibration, 
based on the theory of induced oscillations, are established. Inves- 
tigations of the exciting force may be carried out on oscillating as 
well as on fixed gate models. 

The systematic investigations carried out at the All-Union Scientific 
Research Institute of Hydrotechnics since 1953 have been made on 
oscillating gate models suspended on strain gauges by using highly 
sensitive measuring apparatus. 

This investigational work threw some light on the pulsation of the total 
hydrodynamic load acting on bottom outlet sliding gates and helped in 
studying the problem of scale-modelling. 


A. S. Abelev, Basic Problems of Design and Research of Gate 
Vibration in Hydraulic Structures, Proc, of All-Union Scientific 
Research Inst, of Hydrotechnics named after B. E. Vedeneyev, vol, 
54, 1955. 
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3. INVESTIGATION OF THE TOTAL PULSATING HYDRODYNAMIC 
LOAD ACTING ON BOTTOM OUTLET SLIDING GATES 


The investigations were conducted on oscillating gate models by using 
the methods described in our papers l)» 2), The gate models were sus- 
pended on six wire resistance strain gauges (Fig. l), which enabled 
recording of the pulsation of the total hydrodynamic load acting on slide 
gates both in vertical (strain gauges I and II) and horizontal directions 
(strain gauges III, IV, V and VI). 


The tests carried out. showed that the pulsation of total hydrodynamic 
loads represented unperiodical oscillations as a number of positive and 
negative impulses different in form, value and duration (fig. 2). 

The investigation was intended to determine the value of maximum 
pulsat ion amplitude P rOOlX of the total hydrodynamic load acting on 
gates, and the dominant most frequently repeated oscillation frequency 
of this load for various service conditions of the structure. 


The value of the maximum pulsation amplitude PmClXwas determined 
by photographic recording of oscillations as well as through visual 
observations of deflection of a light indicator on an oscillograph screen^^* 


To determine the frequency range (spectrum ) of recorded oscillations 
and the dominant frequencies . the experimental data were subjected to 
statistic analysis. According to the data obtained.curves were plotted 
of oscillation distribution in periods — X ("f) or in frequencies 

Xm ' '• 



Since the investigations were intended to determine the frequency of 
maximum pulsation, in the analysis of the experimental data, the oscil- 
lations with very small amplitudes were not taken into account and 


See p, 4. 

2) A, S. Abelev, Investigation of Total Hydrodynamic Pulsating Loads 
for Design of Vibration of Bottom Outlet Sliding Gates, Proc. of the 
All-Union Scientific Research Inst, of Hydrotechnics named after 
B. E. Vedeneyev, vol. 58, 1958. 

3) While investigating pulsating loads, the readings of instruments 
corresponding to the average hydrodynamic loads were compen- 
sated by the electric method. 
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considered were only those oscillations, of which the amplitudes were 
not less than (0, 5-0, 3) P'maX. 

Such curves had a pronounced maximum and enabled to determine the 
magnitude of a dominant period and its corresponding dominant fre- 
quency 



It is obvious that the longer the analysed oscillogram is, the more accur- 
ate the results obtained by such analysis will be. As a rule, in oscilla- 
tion frequency tests, the length of each oscillogram was not less than 1500 
mm. which corresponds to a filming period of not less than 30 seconds. 

It should also be pointed out that the natural oscillation frequency of the 
system always exceeded considerably (several times) the maximum 
oscillation frequency and was quite sufficient for measurement of pul- 
sations without introducing an error. 

The investigations described below havte been carried out with a sliding 
gate installed in the inlet section of a conduit as shown in fig, 3. The 
peculiar feature of this installation is that the gate is placed in grooves 
without any intake beam in front of it. 


The gate model was provided with double sided skin plate and bottom seal. 
The latter was designed so that it did not prevent oscillations of the sys- 
tem^). 

Maximum head pressure upon the gate model above the bottom of the 
conduit was in the order of 150-155 cm. 


The research has been conducted at various openings of the gate 




a 


The investigations showed that the value of maximum amplitude and 
pulsation frequency of the total hydrodynamic load, depend, to a consid- 
erable extent, upon the hydraulic conditions in the conduit and location 
of the vortex zone forming behind the gate, A special series of experi- 
ments has been conducted to study this problem, the experiments having 
been carried out under the conditions of relatively small submerging of 
the top of the conduit at constant downstream water level (DWL), 


See p, 5, n. 2, 

2) The downstream water level (DV/L) in the present series of experi- 
ments has been determined at a considerable distance from the out- 
let section of the conduit behind the recovery area. 
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various gate openings, various upstreams water levels and various air 
discharges Q,^ behind the gate. The air discharge varied from CL^O 
to 0-2 O i. e. the value corresponding to the unlimited access of air 
into the conduit K The results of the investigation enabled to outline three 
characteristic regions of hydraulic conditions in an operating conduit: 

region of pressure conditions , with the hydraulic jump roller coming 
over the bottom edge of the gate, the region of transition condition , with 
the roller moving away from the gate but located inside the conduit, and 
the region of pressureless condition , with the roller coming out 
of the conduit, to which air can get freely from the downstream side^). 

To explain this, fig, 4 shows the relation between the maximum ampli- 
tude of the total pulsating hydrodynamic load acting on a gate in the hori- 
zontal direction, and the varying upstream level (UWL). The variation 
of UWL determines the different hydraulic conditions of the operating 
structure. In these experiments, the DWL remained constant. 

The experiments showed that, other factors being equal, the value of 
pulsating load under the pressure condition exceeds that of the pres sure - 
less condition, while the pulsation frequency increases as a rule in the 
course of transistion from the pressure condition to the pressureless 
one. In case of vacuum, air delivery behind the gate results in consid'er- 
able changes of hydraulic conditions and in reduction of the pulsation 
amplitude of the total hydrodynamic load for all those cases when the 
conduit works under the pressure condition; if the conduit operates under 
pressureless conditions, free access of air below the gate determines 
the tendency towards a certain increase in pulsating loads. The experi- 
ments show that with a given gate opening and constant downstream 
water level,the value of maximum amplitude of the total hydrodynamic 
pulsating load is the greatest either on the boundary line between pres- 
sure and transition conditions (at heads lower than maximum) or at 
maximum head pressure in the region of pressureless condition (fig. 4). 

In this connection, when estimating pulsating load and expected gate 
vibration, particular attention should be paid to the hydraulic con- 
dition of operating structure as well as to the value of acting head 
pressure; in certain cases an increase in effective head pressure 
results in changes in hydraulic condition, thus considerably reducing 


Since aeration of the flow is not scaled up, it is obvious that, in the 
experiments with air discharge behind the gate,it is possible on the 
basis of model investigations to evaluate only the quantitative aspect 

of the problem - the influence of air discharge on pulsating load and 
expected vibration. 

See p. 5, n, 2 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


10-A 


1 


1 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


the pulsating load and the possibility of development of dangerous vib 
ration. 


With respect to averaged hydrodynamic loads, the value of maximum 
pulsation amplitude of the total hydrodynamic load can reach 35-60% 
and higher. It should be kept in mind that for the design values of 
total pulsating loadstone should take account not of the amplitude value 
P*mCLX but of the corresponding deviation of hydrodynamic load from 
its total averaged value. However, since this deviation is not the same 
on either side and the location of the averaged load line is not known, 
the load, equal approximately to 65-0, 75) P*noclX may be assumed 

for a tentative design value. This, however, needs further verifica- 
tion. 

The results obtained from these investigations have refuted the general 
opinion that vibration is most dangerous at partial gate openings, that 
is, with 0 < n < 1, 0 and becomes negligible when the gate is lifted up 
to the full height of the conduit, i. e, at n = 1, O, For the installation 
that has been tested by the author (fig, 3) with the gate located in the 
conduit inlet section and with the bottom seal provided, the value of 
maximum pulsation amplitude of total hydrodynamic load reaches the 
highest magnitude at full gate opening (n = 1, O), 

When the bottom edge of the gate is raised above the conduit top, that 
is, with n > 1,0 the value of maximum pulsation amplitude of total 

hydrodynamic load first sharply goes up and with n = 1, 13-1, 14 inc- 
reases, compared to n r 1, 0, by 2, 5 - 4 times, then starts decreas- 
ing and drops to zero at n = 1, 35- 1, 40, A more detailed considera- 
tion of the problem was given in our paper 0. This paper also gives 
recommendations to reduce pulsating loads at n = 1. 0 and n > 1,0. 

The next extensive cycle of experiments conducted by the author was 
aimed at establishing the criteria for determination of PPnQX and 
Not under different operating conditions of deep sliding gates. 

As the above described experiments have shown, other factors 
as to the design and scheme of gate location within the conduit being 
equal, the unknown values "P rnOLX and N[j[ depend to a considerable 
extent upon the hydraulic conditions of operating structure. In this 
connection, we adopted the following approach to the problem in 
question. The search for design criteria has been determined for 
a certain fixed hydraulic condition of the operating structure; then 
were determined the correction factors (condition coefficients) which 
enabled to evaluate pulsation loads for other hydraulic conditions of 


1 ) 


See p. 5, n. 2, 
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the operating structure. 


The pressure condition, i.e. such a condition in which the longitudinal 
boundaries of the flow represent rigid walls, was assumed as a charac- 
teristic hydraulic condition, for which the design relations were deter- 
mined^'. 


In the design charts for the determination of P'mCLX shown in figs. 5 and 
6 the results obtained from investigations are presented in terms of di- 
mensionless coordinates A and n the dimensionless parameter A 
being taken as equal to the following characteristic values: 


(1 

where CO is the area to which the pulsating load is applied. It is evi- 
dent that CO = b^ for the pulsating load acting in horizontal directions 
(see fig* 5) and = for the pulsating load acting in vertical direc- 
tion (see fig. 6). 


It can easily be determined from equation (l) that dimensionless para- 
meter /A is the mean value of maximum strain of pulsating load 

f-h'mny) PmaX 

\r -^^^meCin CO expressed by the corresponding liquid column 

( p* max) mean 


( hpmax) 


ipmaxj mean — ^ 

and related head pressure H, that is 


A 


(hpmax) mean _ 

M 


(p'max)meQn p'max 


( 2 ) 


iTH 

The investigations showed that the relations (?' max)^ = flH) and 

(P'max)v = f (H) are the linear functions of M for all gate openings; 
in this connection, by introducing the dimensionless parameter A con- 
taining H in the denominator, it is possible to replace the design charts 
plotted for the determination of (P‘max)h and (P*max)y by the two uni- 
versal curves respectively 

A, _ (Pmaxlh _ ■r(n) nnri A _ (PmaxV _ r 

(figs. 5 and 6). 


In investigations of the bottom outlet gate position, knowing its dimensions 
and head pressure H we can determine, by using the above 


In this series of experiments the downstream water level elevations 
were measured by the piezometer established in the outlet section 
of the conduit. 
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curves, the values of and (p' lnnC3l>c)v for pressure condition of 

the operating structure at any gate openings 0 < n^ 1, 0, 

The experiments have shown that for pressure condition of the operating 
conduit the value of relative submergence of the conduit top from the 
downstream side does not influence the value of P'tnOlX and the nature of 
oscillations under study, 

A special series of experiments was aimed at studying the above- 
mentioned condition coefficients, that is, relative variation of the maxi- 
mum hydrodynamic load for the case of transition from pressure condi- 
tion , for which the design relations have been obtained (figs. 5 and 6), 
to other hydraulic conditions of operating structure. 

The investigations have shown that for the most unsatisfactory hydraulic 
state within the range comprising the pressureless and pressure condi- 
tions the value (Prnax)h can increase, in comparison with corresponding 
loads fixed for pressure condition, more than three times while (P rnOx)v 
increases more than twice. 

At the same time in the case of transition from pressure condition to the 
pressureless one the values of (P max)h can decrease almost four times 
and those of (Pmav)v more than thrice. Thus, from the above it is 
possible to draw a very important practical conclusion to the following 
effect; to reduce pulsation loads and, therefore, possible vibration of 
a gate, efforts should be made to provide a pressureless condition for 
an operating structure. If it is not possible, measures should be taken 
to ensure a steady pressure condition. It is not recommended to oper- 
ate structures under intermediate hydraulic conditions which cause a 
considerable increase in pulsation loads. 

To conclude the paragraph it should be mentioned that the extensive 
experimental data, relative to the frequency characteristic of the 
phenomenon under investigation, require more accurate definition 
and additional analysis prior to publication. 


4. SCALE-MODELLING OF THE TOTAL PULSATING HYDRO- 

DYNAMIC LOAD ACTING ON BOTTOM OUTLET SLIDING GATES 

When studying total pulsating hydrodynamic load on oscillating gate 
models, the problem of scale-modelling of this phenomenon should 
be investigated. This would make it possible to examine the scale 
effect, that is, the distorting effect of small dimensions of the model, 
and to develop reliable formula for applying scale laws to the proto- 
type. 
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In other words, the problem involves the examination of the theory of 
scale-modelling of this complicated phenomenon and the determination 
of the boundaries of the automodel region. This paper gives a descrip” 
tion of investigations carried out by the author in 1955-1956, with 
regard to the problem of scale-modelling of maximum amplitudes and 
dominant frequencies of the total hydrodynamic pulsating load. It should 
be mentioned that this problem, as far as we know, had never been 
studied before either for the conditions of flow around the gates or for 
any other boundary layer conditions. 

The papers published on the subject deal with the problem of scale - 
modelling of the pulsating pressure in separate points of a streamlined 
profile located on pressure gallery walls behind the gate^) and on the 
dam' s apron^). 

In the investigations described below the problem of scale-modelling 
of the total pulsating hydrodynamic load has been studied by carrying 
a special series of model experiments, the models being built to dif- 
ferent scales (scale series). 

The investigations were carried out on three geometrically similar 
models of slide gates made to scales 1, 1:1, 5 and 1:3 which are further 
referred to as and . The experiments were carried out 

for the pulsating load acting in horizontal direction with various open- 
ings of sliding gates installed in the inlet section of rectangular con- 
duits according to the scheme, being as shown in fig. 7. In addition 
to the similarity of the gate models, there was ensured the complete 
geometrical similarity of acting pressures in headrace, of grooves, 
conduits and tailrace. The basic dimensions of the experimental 
installation shown in fig, 7 refer to the largest model made to scale 
Ml. 

For the two remaining models the above mentioned basic dimensions 
of conduits, gates, heads etc, were reduced by 1,5 and 3 times res- 
pectively. 


^ D. I. Kumin, N. A. Preobrazhensky and G. A. Yuditsky, Scale- 
Modelling of the Pulsating Pressure on Pressure Gallery Walls, 
Proceedings of the All-Union Scientific Research Institute of 
Hydrotechnics named after B. E. Vedeneyev, vol. 52, 1954. 

M.S. Fomitchev, Studies of Scale-Modelling of the Pulsation 
and Hydrodynamic Pressure, USSR Academy of Science, 
Proceedings (isvestiya AN USSR, OTN) Technological Series, 
No. 11, 1956. 
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The frequency of natural oscillation of the system for the three models 
of the scale series considerably (several times) exceeded the maximum 
frequency of oscillations under study and was quite sufficient for pul- 
sation measurements without introducing an error. Investigations were 
conducted under pressure conditions without admitting air behind the 
gate. 

1° ♦ SCALE -MODELLING OF THE MAXIMUM AMPLITUDE OF THE 
TOTAL HYDRODYNAMIC LOAD. - To obtain reliable experimental 
data, testing of the three models was repeated for most of the gate open- 
ings and in some cases was carried out three to four times. 

To investigate the problem of scale-modelling of the maximum pulsation 
amplitudes of total hydrodynamic load Prncix the experimental points, 
obtained for the three scale models, were plotted in the diagram (fig. 

8) where the ordinates give n and the abscissas give the dimension- 
less numbers A, the values of the latter being computed by means of 
the formula (2) according to which for a load acting in horizontal direc- 
tions A equals to 

A p'max p'mox (3) 

IfcoH Yb,h,H 

As can be seen from this diagram the experimental points of the three 
models of the scale series lie within a single curve with an acceptable 
degree of accuracy. As to the scatter of experimental points, it is not 
the result of systematic deviations due to the scale effect and can be 
entirely attributed to the measurement error lying within the limits of 
± 15%. 

The results obtained show that for the three models of the scale series 

A = idem (4) 

Equation (4) shows that the phenomenon under study for the three models 
is taking place in the automodel region. 

It should be pointed out that since the gravity effect of liquid in pressure 
condition is negligible while it is nil in horizontal conduit, similarity 
of the phenomena in this case is determined, as is known, by the iden- 
tity clause of Reynolds criterion, conformance to which ensures iden- 
tity of the values of Euler criterion connecting pressures and velocities 
in similar flows. 

Since the turbulent disturbances in the flow under consideration develop- 
ing in the zone of gate location and causing hydrodynamic pulsating load 


0 
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are large in magnitude and intensity, the automodel region, that is, the 
zone of quadratic resistance, within which similarity is taking place 
despite nonconformance, the identity clauses of Reynolds criterion, is 
accomplished with relatively small numbers Re. 


To confirm the aforesaid and establish the boundaries of the automodel 
region, fig. 9 gives the relations (ReR,n). The diagram (fig. 9) 

is plotted by using the experimental data obtained on the three scale 
models for various H adopted as a parameter. The numbers 
have been computed by means of the formula 

Do =_9L_B 

The analysis of the data given in fig, 9 shows that within the range of 
Reynolds numbers in the experiments for the three models of the scale 
series (2760 ^ 80000) the value of the dimensionless parameter: 


A — P max 


practically does not depend upon Re, which cor- 


responds to the automodel region. In this connection the pulsation 
similarity determined for the three models with respect to their maxi- 
mum magnitudes and dominant frequencies (see parag* 4 , 2^), is found 
out with relatively small values of Re, This enables us to draw an im- 
portant practical conclusion that the investigations involving studies of 
the total pulsating hydrodynamic load acting on sliding gates can be con- 
ducted on models of relatively small scale. 


Proceeding from equation (4) and denoting all the components of equation 
(3) by index “f" for the larger model and by index '^m** for the smaller 
model we can write 


(Pmax)f _ (Prr>ax)m _ A _ .({em 
■ h.rH^ *■ 


Since density of liquid is considered to be constant and thus 
it is not difficult to obtain from equation (3) the following expression for 
the ratios of the maximum amplitudes of total pulsating hydrodynamic 
load in the case of transition from one scale model to the other 

I 


(Pmax)f _ / , _ Up hif 

(P;max)rp M, 

oC 


(P'max) 

therefore 
Then 


= c^% 


e 


( P max)^ - { P max)jpC( 


( 6 ) 

(7) 

( 8 ) 


The result obtained shows that with respect to the maximum amplitudes 
of the phenomenon under investigation, the general principles of the 
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scalc~inodelling theory relative to the autornodel region are completely 
proved . Within the i\a.uge of the Reynolds niin\bers used in tests on 

the three models of the scale scries (2760 4 80000) in the case of 

transition from one model to the other the absolute value of the maximum 
amplitude of total pulsating hydrodynamic load acting on deep sliding gates 
varies in proportion to the cube of the linear scale coefficient and, thus , 
the scale effect for relative values of the amplitude may be considered 
as being equal to zero. 


Let us compare now the results obtained with respect to the scale coef- 
ficient oCPrrtoLX^^^ maximum amplitudes of total pulsating hydrodynamic 
load, with the data of corresponding investigations intended to study the 
problem of scale-modelling of the pulsating hydrodynamic pressure in 
separate points of a streamlined profile^'* L From equation (?), we 
obtain ^ 

P max “ ^ e 

if we go now from Pmax to corresponding average stresses of pulsating 

load (Pnnax)fY,eah= oC(Pmclx)tT\eciinfor the 

above stresses will be apparently equal to: 

_ cCptnax _ cC e 


oC (pmax) mean 


OC 


= cCi 


(9) 


CO oLq 

which is in full agreement with the results obtained from the investiga- 
tions of scale-modelling of the hydrodynamic pulsating pressure in 
separate points of a streamlined profile placed in the flow region with 
intensive turbulent mxing^^’^). 


2° . SCALE -MODELLING OF THE TOTAL PULSATING HYDRODYNAMIC LOAD 
BY FREQUENCY. - In test of scale series was also investigated the problem of 
scale -modelling of the total pulsating hydrodynamic load by frequency. 

To this end the pulsations were recorded on film and the curves of oscil- 
lation distribution by pulsation impulses in periods or in frequencies 
have been plotted on the basis of statistical analysis of oscilligrams as 
described above (see 3). To obtain the reliable research data, pul- 
sation recording in each experiment has been repeated not less than 
three times. 


Fig, 10 presents some of the results obtained from the statistical anal- 
ysis of the oscillograms recorded for the three models by measuring 
cell V at n. a 0,2. The curves of oscillation distribution in periods 

See p. 11, n, 1, 

See p» 11, n. 2. 
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presented in this fig. show that the maxiinums of the 
above curves are shifted in the direction of larger periods with increas- 
ing scale of the model and thus the dominant frequencies of oscillations 


under study are reduc ed respectively. Such curves plotted for 

Tq 

all experiments enabled to determine regularity in variation of the 
characteristic frequencies depending on model scale. 


It should be pointed out that the condition of pulsation similarity accord- 
ing to frequency is determined, as is known, by the identity of the cri" 
terion of homochronism which is also called the Strouhal number. 
According to this, the further analysis of the experimental data con- 
sists in determining the Strouhal number for all scale series experi- 
ments* Comparison of different Strouhal numbers for models made to 
different scales enable to establish similarity of total pulsating hydro- 
dynamic load by frequency and determine the scale coefficient for the 
dominant frequencies. 


The Strouhal number has been determined by using the equation 

= (‘ 0 ) 

in which the height of the conduit *^h*^ was adopted for the character- 
istic linear dimension of pressure flow^)* 

The dominant frequencies and the corresponding Strouhal numbers for 
the total pulsating hydrodynamic load applied to the upper and lower 
parts of the gate were determined by strain gauges III (or IV) and V 
(or VI) respectively^'. 

The Strouhal numbers computed for pulsations are plotted in diagram 
shown in fig. 11 where the abscissas give the Strouhal numbers and the 
ordinates give “n”. 

As the diagram shows, all experimental points corresponding to pul- 
sations, recorded by strain gauges III (or IV) irrespective of the model 
scale, lie within a single curve with an acceptable degree of accuracy 


It is known that in kinematically similar flows the scale of turbul- 
ence which influences the frequency of pulsation under study is pro- 
portionate to the linear dimensions of the flow. 

The symmetrically located strain gauges III and IV gave absolutely 
identical pulsation recordings. The same can be said about sym- 
metrically located strain gauges V and VI, 


1 ) 

2 ) 
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while the points corresponding to pulsations recorded by strain gauges 
V , (or VI) lie within the other curve. As to some scatter of the experi- 
mental points, this, as the diagram (fig, ll) shows, resulted not from 
the systematic deviations due to the scale-effect but from the measur- 
ing error lying within the limits of t 15%, 

The result obtained shows that the Strouhal numbers for all three scale 
models may be assumed as being constant 

Sh =? idem { 1 1) 

On these grounds we could admit that the phenomenon in question for 
all three models is taking place in the automodel region . To confirm 
this and to determine the boundaries of the automodel region figs, 12 
and 13 give relations Sh * “F ( Rbr j ^ ) for pulsations recorded by strain 
gauges III (or IV) and strain gauges V (or Vl), 

The above relationships have been plotted according to the experimen- 
tal data obtained for the three scale models with various n taken as 
a parameter. 

It can be seen from the above diagrams (figs, 12 and 13) that within 
the range of Reynolds numbers adopted in the experiments for the 
three scale models (2760 ^ 80000)^' the Strouhal numbers are 

practically independent of which corresponds to the automodel 

region . 

Proceeding from equation (ll) and denoting the components of equation 
(lO) by the index **f*^ for the larger model and by the index those 

for the smaller model, we can write the following: 

h = idem 


Then it is not difficult to obtain the relation of the dominant frequen- 
cies of total hydrodynamic pulsating load in the case of transition from 
one model to the other: 

Nd.f . _/ -Vtf h m (13) 

Ndm hf 

Since the pressure conduits in scale models connect the two open water 


The values of have been computed as explained above, 

see page 13, 


lO-A-16 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 



1 


Approved For Release 2009/05/06 : CIA-RDP80T00246A008200300002-6 


levels, the elevations of which were predetermined as being propor- 
tionate to the linear dimensions of the conduits, and the discharge 
coefficients of the conduits, as it follows from the experimental data, 
practically did not depend upon RBr , then as is known, in the case 
of transition from one model to another we have 

^ q/v= oCe^ 

Further, if we take into consideration that 


^ (15) 

hf oCe 

and substitute equations (14) and (15) for equation (13), we obtain 



(16) 


In compliance with this we have 



(17) 


The result obtained shows that, as to the frequency characteristics 
of the phenomenon in question, the basic principles of the scale- 
modelling theory relative to the automodel region proved out with 
sufficient degree of accuracy* 


Within the range of Reynolds numbers adopted in the experiments for 
the three scale models (2760 ^ ^ 80000) in the case of tran- 

sition from one model to the other the dominant frequency of the total 
pulsating hydrodynamic load upon sliding flat gates varies inversely 
proportional to the square root of the linear scale coefficient and thus 
the scale effect for relative values of frequencies is equal to zero . 


It must be noted that a similar result has been obtained in the inves- 
tigations mentioned above when studying the scale laws for the 
characteristic frequencies of hydrodynamic pressure pulsation in 
separate points of a streamlined profile in flow regions with inten- 
sive turbulent mixing. 


1 ) 

See p. 11, n. 1. 
See p. 11, n. 2, 
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_ » Some additional investigations and conclusions^ The additional 
investigations showed that the above scaling-up formulas for the total 
hydrodynamic pulsating load obtained under pressure condition and 
expressing the scale laws of the phenomenon under study are true not 
only for pressure conditions but can apparently be used in the case of 
pressureless and intermediate conditions» 

It follows from the above data that scale’-modelling of the total hydro - 
dynamic pulsating load on bottom outlet sliding gates and hydrodynamic 

pulsating pressure in separate points of a streamlined profile placed in 
flow regions with intensive turbulent mixing are dependent on Froude*s 
law of gravitation similarity. The scale laws of pressure pulsation in 
Separate points which are in good agreement with model results were 
proved out by the corresponding field investigations^^# This suggests 
that the scale laws of pulsation of the total hydrodynamic load ) proved 
out by testing the scale model series^ could be applied to the proto- 
type with sufficient degree of accuracy* 


5* GENERAL CONCLUSIONS 

The problem of the total hydrodynamic pulsating load on bottom outlet 
sliding gates and scale-modelling of this phenomenon has been for the 
first time studied in the investigational work carried out by the writer# 

The results obtained enable us to draw a number of important conclu- 
sions and generalisations which are of interest both from the viewpoint 
of finding out some common regularities of the phenomenon under in“ 
vestigation and with respect to their practical application in the design 
practice to evaluate the dynamic loads and expected vibrations. The 
investigations did not expose any particular circumstances which would 
prevent the use of sliding gates for operation under high head pressure 
both at full or partial openings# Taking account of the acting averaged 
and pulsating loads as well as of the above remarks, the sliding gates 
can be designed for high head pressure and will offer reliable service, 
provided they are properly made. At the present moment similar in- 
vestigations with other schemes of location of sliding gates have been 
completed and investigations of tainter and throttle type gates are 
under way. 


1 ) 


See p. 11, n, 2. 
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LIST 

OF SYMBOLS 



A - is 

a dimensionless parameter 



Q - denotes height of gate lifting above the ^ 

bottom of the conduit 

• 


b - denotes width of the conduit 

• 


b, - 

** part of gate to which load is applied 



H - 

** — t 7DWL- acting head pressure 



h - 

height of the conduit 



- 

working height of the gate 



h'p - 

excessive piezometric height 



M - 

** scale of the model 



. 

” number of pulsating impulses of 

equal duration or equal period 



1 

Z 

E 

number of pulsating impulses of 
equal frequencies 



s:m- 

N'* - 

•• summation of impulses considered 

in frequency analysis of oscillograms 

pulsating frequency 


j 

Nd - 

** dominant pulsating frequency 

• 

n = 

a 

h 

relative gate opening 


In our preceding papers the pulsating frequency was denoted by 
symbol V which, however, appeared to be inconvenient since 
it is known that in hydraulics and hydromechanics this symbol 
denotes generally the kinematic viscosity. That is why from now 
and on we shall denote the pulsating frequency by symbol N * 
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Pmciy« PmClX“Pmin “ denotes maximum amplitude of total 


Pmax cind Pmin- ” 

(Pmttx)iii+iv <ind(PmQx)iv+vi - »* 

(Pmax)h Qh(l(Pmax)y- •» 

I 

(PmQy)mean “ " *' 

Q. - » 

0-2 ^ “ 

R - 

Re^ = — - “ 

Sh - » 

T - 

\ - •• 

iX _ 

- •* 

oLq - 


pulsating hydrodynamic loads 

averaged total hydrodynamic loads» 
acting on the gate in different periods 
of time 

maximum pulsation amplitudes of 
total hydrodynamic load acting in 
horizontal direction upon the upper 
lower parts of the gate 

maximum amplitudes of total pul- 
sating hydrodynamic load acting on 
the gate in horizontal and vertical 
directions respectively 

mean value of maximum strain of 
pulsating load 

flow discharge 

air discharge behind the gate 
is hydraulic radius 
Reynolds number 
Strouhal number 

duration or period of pulsating im- 
pulse 

dominant pulsating period 

mean flow velocity in conduit behind 
the point where the jet is widening 

mean flow velocity in section I-I 
under the gate, jet compression 
being ignored (see fig, ll) 

linear scale efficiency 
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denotes scale efficiency of dominant pul- 
sating frequencies of total hydraulic 
load 




scale efficiency of maximum ampli- 
tudes of total pulsating hydrodynamic 
load 


^(Pmax) 


mean 


oC V - 

o6gj - 

t - 
6 - 
V - 

GJ - 


scale efficiency of average values of 
strain of pulsating load 

scale efficiency of velocity 

scale efficiency of area 

volume weight of water 

thickness of the gate 

kinematic viscosity 

area 
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Schema of an oscillating model of the sliding gate. 1 - gate; 

2 - rods; 3 - strings; 4 - cantilever beams for the wire resis- 
tance strain gauges; 5 - tension springs; and Pv hydro- 
dynamic loads acting upon the gate in horizontal and vertical 
directions . 


^ I, Q)sec 

c:«^i ^ I > 




Figure No. 2 

A typical recording of total hydrodynamic pulsating load. 
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Figure No. 3 


Plan and longitudinal section of the experimental installation. 
1 - gate; 2 - conduit; 3 - aeration tube; 4 - diaphragm for 
measurement of air discharge. 



Relations (Pmax) jj, = f ^ “ j" ^ 

and (PnnClx)y^=|(U.)at downstream water level 22.0 cm; 

n " 0.2; Q. 2 5 Uupstream water level; (a) - region of 

pressure condition; (B region of transitional regime; 

© - region of pressureless regime. 
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Figure No. 7 

Plan and longitudinal section of a scale series experimental 
installation. Model made to scale . 1 - pressure tank; 

2 - gauge glass; 3 - gate; 4 - conduit; 5 - piezometer; 

6 - tail race flume; 7 - shutter; 8 - discharge. 
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Figure No. 8 

Relations A - f (n) plotted according to the experimental data. 
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Figure No . 9 

Relations A«^(Re^»n) plotted according to the experimental 
data. 
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Figure No. 10 


Curves of oscillation distribution in periods for pulsations 
recorded by strain gauges V. 



Figure No. 1 1 

Relations Sh = f (n) for pulsations recorded by strain gauges III 
(or IV); and V or (VI). 1 - pulsations recorded by strain gauges 

III (or IV); 2 - pulsations recorded by strain gauges V (or VI). 
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